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Low-Drag and Suction Airfoils 


The Eleventh Wright Brothers Lecture 


SYDNEY GOLDSTEIN* 
University of Manchester, England 


(1) INTRODUCTION 


nation LIKE FIRST TO SAY how much I appreciate the 
honor of being invited to deliver the Wright Brothers 
Lecture. To anyone whose work is associated in any 
way with the aeronautical sciences, it must be a source 
of pride and gratification to be invited to be a chief 
participant in one of the greatest occasions in the world 
of aeronautical research, an occasion designed in honor 
of the two great pioneers, Wilbur and Orville Wright. 
In my own case these feelings are shared with a feeling 
of humility and of my own unworthiness for the task. 

The story I am to tell in this lecture is still incom- 
plete. It is largely the story of the researches carried 
out in Great Britain arising from an idea due to Dr. 
A. A. Griffith. This idea is at once very general and 
very simple. When, at a sufficiently high Reynolds 
number, a fluid flows along an immersed solid surface, 
a boundary layer is formed; if the flow in the boundary 
layer cannot be controlled, harmful results may follow; 
the boundary layer may be removed by suction in an 
attempt to avoid these harmful results. Griffith’s idea 
was that if boundary-layer suction is to be applied, the 
shape of the solid surface should be specially designed. 
It is the application of this idea to airfoil design which 
forms the main subject of this lecture and it may be as 
well to stress that this does not include a discussion of 
the effects of suction on conventional airfoil shapes, or 
on any shapes not specially designed for the purpose; 
such a discussion should really be separate, since the 
amount of common ground is not large. 

Part of the story has been excellently told already by 
Mr. Relf in his Wilbur Wright Memorial Lecture to the 
Royal Aeronautical Society in 1946: I shall try to fill 
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in some of the details, as far as they are known at pres- 
ent. 

In assessing the probable performance of a suction 
wing, the state of the boundary layer (turbulent or 
laminar) at the slot through which part of it is to be 
sucked away is of great importance. It is therefore 
advisable to spend a short time on a discussion of the 
principles of the ordinary laminar-flow or low-drag 
wing, and of the present position with regard to obtain- 
ing laminar boundary layers in flight. Some mention 
of the principles is moreover necessary in order really 
to understand why a suction wing gives a theoretical 
prospect of better performance than the best low-drag 
wing without suction. 

The work described in this lecture was all carried out 
in the Aerodynamics Division of the National Physical 
Laboratory or at the Royal Aircraft Establishment, and 
I am indebted to the Directors of these establishments 
and to the relevant government departments (the 
Ministry of Supply and the Department of Scientific 
and Industrial Research) for permission to use the re- 
sults. I should also like to thank Mr. W. A. Mair for 
assistance with the preparation of the diagrams. 


(2) Low-Drac WINGS 


The principle underlying the design of low-drag wings 
is simply to design to keep as much of the boundary 
layer laminar as possible without boundary-layer sepa- 
ration before the trailing edge. The transition to tur- 
bulence should be as far aft as possible, without any 
risk of the boundary layer separating; apart from stall- 
ing incidences, this risk is greatest, at the Reynolds 
numbers of flight, at moderately large lifts with the 
boundary layer made accidentally turbulent at the 
leading edge by some form of accidental roughness. 
The whole matter is therefore bound up with transition 
to turbulence in boundary layers and with turbulent 
boundary-layer separation, both matters which, to this 
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day, are very imperfectly understood. Although, 
thanks largely to the excellent wind-tunnel researches 
by the N.A.C.A. in this country, much is known of the 
performance and behavior of low-drag wings, the basis 
for a complete physical understanding of the processes 


involved has not yet been constructed. Mathematical . 


researches in Germany and experimental researches at 
the National Bureau of Standards, corroborated by re- 
searches at the California Institute of Technology, have 
certainly laid a basis for understanding the instability 


, of laminar boundary layers, but it is abundantly clear 


that this instability is not itself the same thing as, nor 
need it necessarily lead to, transition to turbulence. 
Apart from the deeper physical reasons, this is clear 
from the fact that low drags, indicating a large extent 
of laminar boundary-layer flow, have been obtained in 
wind tunnels in this country up to Reynolds numbers of 


30 to 50 millions, and also from the flight experiments 


at the Royal Aircraft Establishment, Farnborough, 
mentioned later. The physical processes leading to 
transition are still largely unknown, and the same is 
certainly true of turbulent boundary-layer separation. 
We are therefore at present forced to rely on certain 
overall experimental results, which fortunately give 
us much of what is needed for practical purposes. 

The first such result is that on a smooth, true surface 
a boundary layer may be kept laminar if the velocity 
outside the boundary layer is increasing, and the pres- 
sure is decreasing, in the direction of flow, but that 
transition to turbulence takes place very soon when 
there is an adverse pressure gradient. Thus low-drag 
wings should be designed so that, at the design lift 
coefficient, and over as wide a range of lift coefficients 
as possible, the velocity outside the boundary layer 
rises from the front stagnation point to a position as far 
back along the wing as possible. A second overall ex- 
perimental result, however, must now be taken into 
account at once—namely, that the larger the unfavor- 
able pressure gradient the greater the risk of boundary- 
layer separation. The further back the position of 
maximum velocity, and therefore of maximum suction, 
the smaller is the fraction of the chord gver which the 
pressure recovery, which must occur before the trailing 
edge is reached, can take place; the greater therefore, 
for a given maximum suction (which is related to the 
thickness of the wing), is the unfavorable pressure 
gradient over the rear of the wing; and so the greater 
is the danger of boundary-layer separation. For any 
normal wing thickness there is therefore a limit imposed 
by these considerations on the position of maximum suc- 
tion. From a consideration of the American experi- 
mental results, M1. Richards! has deduced rough rules 
for the maximum rate of velocity decrease that should 
be allowed over the rear of a wing if the danger of tur- 
bulent boundary-layer separation at moderately high 
lifts, with the wing accidentally roughened and the 
boundary layer made turbulent at the leading edge, is 
to be avoided; if U is the velocity of the wing relative 
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to the undisturbed air, c the chord of the wing, and q 
the relative velocity of the air just outside the boundary 
layer and the wing at a distance x along the chord from 
the leading edge, then, for example, a rough criterion 
is that at a lift coefficient of 0.95 the nondimensiona] 
adverse velocity gradient —(c/U)(dg/dx) at seven- 
tenths of the chord from the leading edge should not 
exceed 1.27, at any rate when the maximum suction js 
at half chord. 

In practice a further limitation is set, not only to the 
length of laminar boundary layer which can be achieved, 
but also to the length that should be aimed at, by rough- 
ness and waviness of the surface. Beyond certain quite 
small limits it is known from experiment that even with 
a favorable pressure gradient roughnesses, or waviness 
due to manufacturing limitations, will induce, a transi- 
tion to turbulence in the boundary-layer flow. If 
transition occurs forward of the position of maximum 
velocity, the drag of the wing will be greater than if the 
maximum velocity had been designed to be at the posi- 
tion of transition, since there will be a greater area with 
a large velocity outside the boundary layer in the turbu- 
lent region, and this produces increased skin-friction 
drag. The aerodynamic design of the wing and the 
probable state of the surface should therefore be 
matched so that the transition occurs at, or just aft of, 
the position of maximum suction. 

The properties of low-drag wings that have been 
mentioned are all televant to a consideration of suction 
wings; other properties, such as maximum lift and the 
behavior of control surfaces, need not detain us, though 
some remarks on designing for maximum lift will be 


made later. 


(3) Tests on Low-Drac WINGS 


Two series of flight tests have been carried out at the 
Royal Aircraft Establishment. One? was on a King- 
cobra (P 63), the wing having an N.A.C.A. 662x-116 
section at the root and 662x-216 section at the tip, so 
the maximum suction at the design lift coefficient is at 
60 per cent of the chord from the leading edge and the 
maximum thickness-chord ratio is nearly 16 per cent. 
The profile drag was measured at a section 137 in. from 
the aircraft centerline. The results are summarized 
in Fig. 1. With the surface in the state in which it was 
received a drag coefficient of 0.0075 to 0.0080 was 
found; the reduction of roughness by smoothing yielded 
a drag coefficient of 0.005 to 0.006 near C, = 0.2, 
with sharp rises on both sides of that lift coefficient. 
A test portion about 3 ft. in span was then very care- 
fully prepared in order to decrease the waviness of the 
surface. The test portion was cleaned down to the 
metal surface, and was given two coats of primer paint, 
and then a paint-type filler was sprayed on; local areas 
found to be low were specially filled, and the wing was 
rubbed down with special curved rubbing blocks with 
strips of carborundum paper (grade 280C) held over the 
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As received at RAE. 
w Roughness reduced 
0-006 
| Waviness and 
punt Roughness raduced 
» 
Lowest drag obtained 
5 0-002 
3 Approximate 
14 Reynolds NO xiO4 
02 04 0-6 
Aircraft lift coefficient. 
Fic. 1. Measured profile-drag coefficient of a section of a wing 


of a Kingcobra aircraft. 


curved surfaces. The drag was now greatly reduced, 
as shown in Fig. 1; a drag coefficient of 0.0028 was ob- 
tained at C, = 0.1, and the lowest drag recorded (with 
the surface polished) was 0.00245. The surface wavi- 
ness was measured with an Ames dial, which could be 
read to 1/500 in., and which was mounted on a base 
supported at points 2 in. apart. The measurements on 
the test section as received and in the final condition 
are shown in Fig. 2. The measurements were static 
measurements made on the ground; deflections in 
fight were not measured. It appears that on a wing of 
this type, at the Reynolds numbers of the tests (up to 
18 millions), a very low drag is obtainable if the fluctua- 
tions on the front portion of the surface can be reduced 
to 1/1 oo0 in. measured on a 2-in. gage; larger fluctuations 
may be allowed further back. 

It may be noted that the use of paint and filler, by 
itself, does not necessarily reduce waviness, since they 
may be rubbed down in such a way that the final sur- 
face, while smooth, follows the waviness of the metal 
surface. 

In order to obtain the results shown in Fig. 1 it was 
necessary not to pick up flies and other insects during 
the take-off and climb. If encountered, flies would 
stick to the surface, generally over the first 10 per cent 
of the chord; in the words of the original report, ‘‘the 
fies on hitting the surface disintegrate and are found 
in the form of small parts of the flies’ bodies stuck to the 
surface with a particularly potent glue.’’ Transition 
to turbulence is, of course, produced by the resulting 
“roughness,’’ and the drag may thereby be raised by 50 
to 100 per cent. In the earlier months of the year 


(March and April) flies could be avoided by flying 
before 9 a.m. G.M.T., but with the approach of warmer 
weather some other method had to be found. A large 
sheet of paper, 3 ft. wide, was placed over the leading 
edge and covered the test portion on the top and bot- 
tom surfaces to 30 per cent of the chord from the leading 
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edge; the paper was kept in position by adhesive 
tape at the edges; a loop of string round the paper at 
the leading edge had its free end in the cockpit. No 
flies were encountered above 5,000 ft., and when this 
height was reached the pilot pulled the string and tore 
the paper along the leading edge; the wind then blew 
off the two halves of the paper from the top and bottom 
surfaces of the wing, leaving a clean surface free from 
parts of flies’ bodies. 

Apart from flies no other trouble was encountered. 
No trouble resulted from the effects of dust or water. 
Water may be used to wash the wings if they are pol- 
ished, and rain has no subsequent effect on the value of 
the drag. Polishing the surface has very little effect 
on the drag, but the use of wax polish is advisable if 
the wings are to be washed. 

Once the surface had been prepared it did not dete- 
riorate (except at the skin-joint) over a period of six 
months, during which 50 flights, occupying about 40 
hours, were made and the aircraft was exposed to 
weather for about 200 hours. During the remainder 
of the time it was kept in a hangar, but no special pre- 
cautions were taken against weathering.* The dif- 
ficulty at the skin-joint, where a gap between the two 
skias of 5 to 10 thousandths of an inch was filled with 
filler, which had to be rubbed down frequently, was 
probably caused by variations in the gap under load, 
and could be avoided by redesigning the joint. 

Experience on another aircraft, which in fact was left 
in the open for over 2 years, suggested that quite long 
exposure to weather would not by itself cause serious 
deterioration of the surface. 

There are, of course, other aspects of the problem of 
the maintenance of low-drag wings besides the effect of 
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Fic. 2. Surface waviness measurements with a 2-in. gauge on 
the test section of the wing of the Kingcobra. ------ Original 
condition. Final condition. 
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Fic. 3. Measured profile-drag coefficient of a section of a 
moderately low-drag wing fitted to a Hurricane aircraft. First 
series, as received; second series, with waviness reduced. 
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Fic. 4. Surface waviness measurements with a 2-in. gauge on 
the test section of the wing of the Hurricane. ---- Original 
condition. ~——Final condition. 


weathering. It is worth recording, therefore, that the 
R.A.E. Staff report® that efforts to impress the ground 
crews with the importance of avoiding damage during 
handling and servicing were completely successful and 
that ‘‘the ground crews became in fact enthusiastic 
custodians of the surface finish.” 

The second series of flight tests‘ was on a Hurricane 
with specially fitted low-drag wings. These wings 
were designed, by methods described in the next section, 
to be of moderate low-drag type, with maximum suction 
at mid-chord, and with a maximum thickness-chord 
ratio of nearly 18 per cent at the root and nearly 15 per 
cent at the tip;’ the test section was just over 17 per cent 
thick. The lowest drag to be expected on this wing is, 
of course, greater than on the Kingcobra. Fig. 3 shows 
the drag measured on the wing as first received from 


JOURNAL OF THE AERONAUTICAL SCIENCES—APRIL, 


1948 


Messrs. Armstrong Whitworth Ltd. and after reduction 
of waviness by the firm by a technique similar to that 
described above. The results of measurements of the 
waviness by a 2-in. gage are shown in Fig.4. With this 
moderately low-drag wing fair results are obtained, even 
with considerable waviness, over a moderate range of 
lift coefficients and up to Reynolds numbers of about 
14 millions. When the deflections are reduced to about 
+1/, 9 in. over the front portion of the surface the 
drag coefficient is about 0.0044 over the whole range of 
lift coefficients from 0.1 to 0.5, with a slight increase at 
the lowest lift coefficients and highest Reynolds num- 
bers (which may have been due to the fact that the de. 
flections exceeded 1/1009 in. at some places on the sur- 
face). Otherwise—with regard to both trouble with 
flies and satisfactory maintenance—the experiences 
with the Kingcobra were repeated with the Hurti- 
cane. 


It is now possible to find the position of transition to 
turbulence in a boundary layer along the whole span of 
a wing in flight by a chemical technique due to Mr, 
W. E. Gray.’ A photograph illustrating the results ob- 
tained in this way on the Hurricane is shown in Fig. 5; 
on the lighter parts of the wing the boundary layer is 
laminar. Fig. 5 shows part of the upper surface of the 
Hurticane after a flight at 205 m.p.h. at C, = 0.22 (R= 
14 X 10°); the general transition to turbulence is at 


Fic. 5. Transition to turbulence in the boundary layer of a 
wing of the Hurricane in flight. The boundary layer is laminar 
on the lighter parts of the wing. The general transition to tur- 
bulence is at 55 per cent of the chord from the leading edge, with 
“‘wedges”’ of turbulent flow further forward. (Cz = 0.22, R = 
14 X 108.) ; 
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55 per cent of the chord from the leading edge, with 
“wedges’’ of turbulent flow further forward. 


The conclusions from the flight experiments may be 
summarized as follows. Low drag may be obtained in 
fight if an accuracy of the surface over the relevant 
parts of *'/1,00 in. on a 2-in. gage can be achieved. 
Although the best method of obtaining the required ac- 
curacy of surface is not yet clear, such an accuracy 
could be both achieved and maintained without exces- 
sive difficulty. Particular attention to joints and rivets 
js still needed when maintenance problems are con- 
sidered. The problem presented by the adherence of 
flies and other insects on the surface near the leading 
edge is serious, and no practical solution has yet been 
found. 


(4) MATHEMATICAL AIRFOIL THEORIES 


If the boundary layer does not’separate, the velocity 
just outside the boundary layer on an airfoil is, except 
at and near the trailing edge, very nearly equal to the 
velocity calculated from the potential flow at the lift 
experienced by the airfoil. If, on the basis of experi- 
mental knowledge, we are prepared to guess at the 
position of transition to turbulence and the occurrence 
or nonoccurreuce of turbulent separation and at the 
reduction of the lift from its theoretical value on the 
Kutta-Joukowski hypothesis, we can forecast many of 
the aerodynamic properties of an airfoil section when we 
have calculated values, according to potential theory, of 
the velocity on the surface. ‘It is also an advantage to 
be able to design sections of airfoils which will have spec- 
ified velocity distributions, or velocity distributions of 
a specified type, on the surface. In this way we can 
design and choose promising models for experimental 
tests in wind tunnels and, eventually, for flight tests if 
the wind-tunnel results are satisfactory. It should be 
pointed out at once, however, that our ability to guess 
at all correctly about the occurrence of turbulent separa- 
tion is very limited, so that, for example, maximum lift 
cannot be predicted with any confidence. Neverthe- 
less, such theoretical calculations as we can make, and 
particularly the solution of the inverse problem of de- 
signing a section with a specified velocity or pressure 
distribution, provide a valuable tool in investigating 
low-drag wings according to the principles discussed in 
Section 2. 


A solution of the first problem, that of computing 
the velocity distribution on the surface of a given airfoil 
section, was given a good many years ago by Theo- 
dorsen,* and by Theodorsen and Garrick.’ This is a 
full and complete solution, but careful computation is 
required, especially for cambered airfoils, if the com- 
puted points are to be smooth and accurate; I have 
seen computed points that looked more like experimen- 
tal points with a fair amount of scatter. In its original 
accurate form the theory is not easy to use for the solu- 
tion of the inverse problem—of designing a wing section 


with a velocity distribution on the surface of specified 
type—or for forming rapidly judgments on the probable 
relations between the geometric and aerodynamic prop- 
erties of the wing. It is therefore well worth while to 
approximate Theodorsen and Garrick’s results. 

For the first, and crudest, approximation we linearize 
all the formulas completely, retaining only quantities 
of the first order in the thickness, the camber, and the 
incidence of the wing; for a symmetrical airfoil at 
zero incidence this leads to the formula 


g/U=1+8 
where 
= — (1/m)P (E)dé/(E — x) 


U is the relative velocity of the wing and the undis- 
turbed air, qg is the relative velocity at a distance x from 
the leading edge measured along the chord in fractions 
of the chord length, and y,(x) the half-thickness at 
that position also measured in fractions of the chord 
length. The integral is an improper integral, and P de- 
notes that the principal value is to be taken. The 
prime denotes a derivative. 

This formula, obtained by linearizing Theodorsen 
and Garrick’s accurate solution, may also be obtained 
by physical reasoning, since we may suppose that we 
could allow approximately for the variation in breadth 
of the wing by supposing sources, of strength 2Uy,’ (x) 
per unit length, placed along the chord; the same 
formula would be obtained in this way, and was so ob- 
tained by Squire. Now the leading edge is a stagnation 
point, and so is the trailing edge if it is not a cusp; and 
it is clear that a distribution of sources along the chord 
that would accurately produce the disturbance flow 
cannot extend right up to the leading edge or the trail- 
ing edge, whereas the approximate distribution does so. 
For this reason this first approximate expression must 
be wrong at and near the leading and trailing edges, and 
in particular will not give stagnation points at those 
places. 

When we find the first, completely linearized, approxi- 
mation for a cambered airfoil at any lift, the additional 
terms (found as before by approximating to Theodorsen 
and Garrick’s accurate solution) may be shown to be 
exactly those found from the vortex-sheet theory of in- 
finitely thin airfoils, as set out, for example, in Glau- 
ert’s Aerofoil and Airscrew Theory. Let a be the inci- 
dence and —§ the theoretical no-lift angle, and let the 
lift coefficient be taken as 


Cr = a sin (a + £) 
let 
x = 1/,(1 — cos @) = sin? 1/0 


so that @ is 0 at the leading edge and 7 at the trailing 
edge; let y.(x) be the ordinate of the centerline, 
measured in fractions of the chord length, and let 


| | 
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dy,/dx = F(#) = >> A, cos (0<@< 7) 
n=0 
so that 
= (1/m) fo" F(t)dt, 
and 
A, = (2/r) f'*F(t) cos nt dt 

forn 2 1. Also let Cz, be the lift coefficient de- 
fined by 


[(1/ao) + (1/2m)]Crom, = At 


Crop. May be shown to be the only lift coefficient 
for which, on a linear theory, there is no infinity at the 
leading edge in the velocity distribution for an infin- 
itely thin airfoil of the shape of the given centerline; 
it is also approximately equal to the lift coefficient at 
which the forward stagnation point is at the leading 
edge. 


Then the first approximation is 
1/1 1 
g/U=1+ 8, + Cr — Crop. ) X 


2 Na 2)” 2 


where g, was given above and 


a= - F(t) cot 5 (t — 0)dt = 


>A, sin n 6 
n=1 
When Cz, = Ci», and a = 2, g/Uis 1+ g, gy 
simply. The positive signs refer to the upper surface 
and the negative signs to the lower surface. 


The no-lift angle is given by 8 = 1/24; — Ap, and the 
moment coefficient at zero lift by (Cir)o = m(Az — Ai) + 
4. 


We have assumed that the theoretical no-lift angle is 
sufficiently accurate; the formulas may be adjusted 
to allow for any discrepancy. 


The first approximation will, as we have seen, be 
wrong near the nose and perhaps also near the tail; 
the second approximation consists merely of mserting a 
factor, obvious from the accurate solution, to remedy 
this. Approximation II is 
q + 1/5Cy2)|sin 


w+ sin? 0)” | 


1 + 1/2Co Cr (2 
(y? + sin? + 2\do x 


(1 + cos 0) | 


where 
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Fic. 6. Computed velocities on the surface of an NACA 0012 
airfoil section at zero incidence. ———Accurate computation. 
© Approximation III. ----- Approximation I. 


= 2y,/sin 0 
= 2y,/sin 0 


For the third approximation we attempt to come as 
near as we can to the accurate solution while still’ keep- 
ing the computations simple. With the upper signs 
referring to the upper surface, and the lower signs to the 
lower surface, Approximation ITI is 


= |AF sin (0 + = 8) + BF cos (6 + 8) * 
CF 


where 
F = &(1 + e')/(Y* + sin? 6)” 


e = €,(0) + €,(0) 
e’ = = 


the prime denoting a derivative, 


€,(0) = —(1/2x)P cot 1/2(t — 0)dt 


= cot 1/o(t — 0)dt 
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Fic. 7. Computed velocities on the surface of a symmetrical airfoil 24 per cent thick, EQH 2450, at zero incidence. 
Approximation II. 


computation. © Approximation II. x: 


note that y,(@) is to be taken as an even function and 
¥.(@) as an odd function of 6, and 


A = cos (a + B) = [1 — (C,?/a0?)]'" 
B = sin (a + 8B) = Cx/a 
C = Cy/(2ne™) 


Thus on Approximation III the calculation of all the 
subsidiary quantities, y, €, Co, is still linear, and the 
results of the computations of these quantities can all 
be scaled and added to allow for variations of thickness 
and camber. The final formula is not linear, but of the 
three coefficients A, B, C, only A is not linear in the 
lift coefficient, and results at other lift coefficients are 
rapidly found from results at any two. We may note 
that slightly more accurate values for g/U are usually 
found if the factor e® in F (but not in C) is replaced by 
exp (Cy + [e’y]), where 


[ev] = + * 


and for a routine procedure the use of this factor is 
recommended. 

For ease of reference, the original nomenclature of 
Theodorsen and Garrick has been largely retained. 
Mr. W. Chester, a research student at Manchester, has 
however recently extended these approximations to the 


Approximation I. 


calculation of the streamlines and the velocity at points 
not on the surface, and we have had to change the 
notation, since y is required for the stream-function. 


It is of interest that for nearly all airfoils for which 
algebraic formulas have ever been suggested, the quan- 
tities g,, Zi, €s, € can be found in closed form in terms 
of simple functions. Examples are given in the orig- 
inal reports.® 


In Approximation III the restriction to small inci- 
dences has completely disappeared. It is still to be 
expected, however, that the results would not apply 
very well for large thicknesses and cambers. A num- 
ber of very careful exact computations have been car- 
ried out for comparison, and the results are rather sur- 
prisingly good at large thicknesses; the effects of large 
camber are somewhat more marked. Some comparisons 
are given in Figs. 6, 7, and 8. In Fig. 6 the results 
(other than Approximation II) are shown on a large 
scale for NACA 0012 at zero incidence; Fig. 7 shows 
results at zero incidence for a symmetrical 24 per cent 
thick airfoil with its maximum thickness at mid-chord; 
and Fig. 8 the results for the upper and lower surfaces, 
separately, for a 12 per cent thick airfoil with 4 per 
cent camber, maximum thickness and maximum cam- 
ber being both at mid-chord; C, = 4.8 sin (a + 8) 
has been assumed, and the results are shown for C, = 


"0 
| 
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Fic. 8. Computed velocities on the surface of an airfoil 12 per cent thick with a 4 per cent camber, EQH 1250/4050, at Cz = 
0.6. Cr/sin (a + 8) has been assumed equal to 4.8. ———Accurate computation. © Approximation III. x: Approximation II. 


Approximation I. 


0.6. In all cases, the comparisons at other lift coef- 
ficients are rather similar. 


Before we leave this branch of the subject, we may 
remark that, in terms of 0, g, sin 0, ¢,, and g,; have 
Fourier sine series that are conjugate to the cosine 
series for 2dy,/d0, y,, and dy,/dx, respectively, and ¢, 
has a cosine series conjugate to the sine series for y,, the 
integral results previously written being alternative ex- 
pressions of these facts; but the series should not be 
used for computation unless they terminate after a 
small number of terms. [The series conjugate to 


(A, cos nO + B, sin nO) is, by definition, 


n=0 
sin — B, cos 

When we turn to the inverse problem—of designing 
an airfoil section with a specified type of velocity or 
pressure distribution on its surface—it is found that the 
best procedure is to design the airfoil section, and at 
the same time to find Approximation III to the velocity, 
from an assuined knowledge of Approximation I. In 
fact a symmetrical airfoil, or the fairing for a cambered 
one, can easily be designed from a knowledge of the 
quantity g, defined above, and centerlines may be de- 


signed from a knowledge of the quantity g; With 
some practice, however, it is not too difficult to make a 
moderately good guess (which may in any case be im- 
proved later, if required) at rough values of the factor 
by which Approximation I is multiplied to yield Ap- 
proximation II, at any rate in those regions, particularly 
near the nose, where Approximation I is known to be 
considerably in error. If the given velocity distribu- 
tion is regarded as Approximation II, Approximation I 
may therefore be at any rate roughly found. Moreover 
C, is known from the pressure distribution; theoreti- 
cally Cz,»,, may be determined directly, if ao is assumed 
specified, from a knowledge of the velocity at the lead- 
ing edge; but in practice this approach is difficult, and 
it is better, having estimated the first approximation to 
q/U on the top and bottom surfaces, to add to find g, 
and subtract to find 
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(= +. =) Crop. = Ar = 2 


to find Cz,5;,, and hence g;. 

To design a symmetrical airfoil or fairing fii a 
knowledge of g, we may either use the knowledge, di- 
rectly, that if y, is.the half-thickness, —2dy,/d0 has a 


cosine Fourier series, without a constant term, con- 


jugate to the sine series for g, sin 0, so that 


dy, 1 1 
— J ~(t — 
2 g,(t) sin cot )dt 
and then integrate for y,, or we may integrate first, 


writing 


= sin = gs dx 


and then prove that —2y, has a Fourier sine series con- 
jugate to the Fourier cosine series for G(@), so that 


2y, = (1/2r)P cot 1/2(t — 0)dt 


G being taken as an even function of 0. 
If the radii of curvature at the leading and trailing 
edges are py and pr, respectively, then 


(2pz)'* = (1/m) (1 + cos A) dt 
= (1/7) (1 — cos 


There must be certain restrictions og g, so that the air- 
foil contour does not cross itself; the most important 
are that the above expressions for (2p,)'* and (2pr)'”? 
should be positive, with the proviso that for a sharp 
trailing edge pr = 0; it may be shown that if g, is 
finite at the trailing edge, then a sharp trailing edge 
must be a cusp. 

To calculate Approximation III to the velocity at the 
surface of the airfoil so designed we use the formula 


above; yw, is found from y, in the usual way, but C), 
é, €,’ are more easily found from g, and G. In fact, 
= 

and 


€, = G(6) cosec 6 — Cy tan 1/20 


If g, is finite at the leading and trailing edges, «, = 0 
there; ¢,’ is found by differentiation, and a limiting 
process shows that, at the leading and trailing edges, 


€,'(0) = */2[g.(0) — Co] 
= — Co] 


It is shown in the original report that y, may be 
found in closed form in terms of simple functions for a 
large class of assumed values of g,. A particularly 
simple case is that of the so-called “‘roof-top’”’ airfoils 
(of which an example is shown in Fig. 9), where g, in- 
creases linearly from a value a at x = 0 to a value bat 
* = X,, and then decreases linearly to a value cat x = 
1; then be 
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Fic.9. Example of ‘roof-top type’ of symmetrical low-drag air- 
foil, and the theoretical velocities on its surface at zero incidence. 
----- Approximation I. ———Approximation III. 


Vs = + + che 


where fo, fi, fe are fixed functions of @ for a given Xj, 
and may be tabulated once and for all for a set of values 
of X;. The radii of curvature at the leading and trailing 
edges may be found in terms of a, b, c, X; before any 
other work is done, and it may be verified that suitable 
values of a, b, c have been chosen from this point of view. 
Of course, the values of a, b, c, X; also tell us about the 
velocity gradients we have chosen in the accelerating 
and decelerating regions. Although it is possible to de- 
sign better airfoils, a preliminary study of these rather 
simple cases helped us to get our ideas clear. 


When we turn to problems of-designing centerlines 
from a knowledge of g,, and of calculating Approxima- 
tion III to the velocity at the surface and other proper- 
ties of a cambered airfoil designed in this way, we first 
note that the Fourier cosine series in 6 for dy,/dx — Ao, 
where y, is the ordinate of the centerline and A» the 
constant term in the Fourier cosine series in 6 for dy, + 
dx, is conjugate to the sine series for —g,; so if g,(x) is 
given as a function of x, and g;(sin? '/20) is taken as an 
odd function of 8, 


gi(€) 


— Ax may be found by integration, and Ay from the 
condition that y, = Oatx = 1. 


However, y, and A» may also be obtained explicitly, 
since it may be proved that if 

Gx) = 

and 

A, = (2/m) So" gs sin 6 d0 = (4/m)G,(1) 


then 


; 
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_ x(1 — 2) 4G,(1) 
iA E log, x + (1 — x) log, (1 — »| 


and 


_ 1 — G1), 


The formula for y, is a transformation of the integral 

- that expresses that ¥.(= 2y,/sin @) has a Fourier sine 

series conjugate to the cosine series in @ for [1/24,0 — 
2G,(x)] cosec 6. 

To calculate Approximation III to the velocity we 

find y, from y,, but ¢, from g; by use of the formula 


€.(0) = 2[Gi(x) — cosec — Ao 
The no-lift angle is given by 1 
B = €(m) = 1/241 — Ao 
and the moment coefficient at zero lift is 
(Cu)o = — So'(4x — 1)gi(x)dx 
also 
€(0) = —1/241 — Ay 


Ws, €s, and Cp are found from the design of the fairing, so 
all relevant quantities are known for the use of the 
formulas written previously. The formulas for the 
radii of curvature at the leading and trailing edges are 
the same for a cambered as for a symmetrical airfoil. 
It may be noted that g; must be zero at the leading 
and trailing edges if dy,/dx is to be free from singulari- 
ties and e,’(@) is not to be discontinuous at 6 = 0 and z. 
- In the calculation of y, formulas may be found in 
closed form in terms of simple functions for a wide 
variety of assumed forms for g;;_ the results for special 
cases have, in fact, been given by a number of authors. 
The remark made in considering the direct problem, 
that the Fourier series should themselves be used if 
they terminate after a few terms, holds for the inverse 
problem also. 
A large part of the work in either the direct or the 
inverse problem consists in finding the sum of the 
Fourier series conjugate to the series for a given func- 
tion. Although in many useful cases formulas may be 
found for the sums in finite terms, the best use cannot 
be made of the theory unless the given functions may 
be varied arbitrarily, and the conjugates must then be 
found numerically. The same kind of computation 
has to be made in the exact theory, and a computa- 
tional formula was given by Theodorsen and Garrick. 
All suggested methods have this in common, that they 
suppose the given function known numerically at a 
number of equaily spaced points in the interval (—7, 
m). Theodorsen merely took the function as constant 
in each subinterval; this is definitely too crude, and 
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we worked out formulas in which the given function 
was approximated by a quadratic expression (i.e., its 
graph was approximated by a parabola) in each sub. 
interval. The results were, of course, much more accu- 
rate, but the formulas have been superseded by others 
derived by a method due to Germain,’ who approxi- 
mates the given function by a trigonometrical poly- 
nomial, or finite Fourier series, coinciding with the 
given function at the points where its value is known 
numerically. Watson’ has shown that general formu- 
las may be found for the constants occurring in the re- 
sulting formulas; moreover, formulas may be found for 
computing the derivative and the integral of the conju- 
gate in terms of the values of the given function, with- 
out its derivative or integral having to be found, and 
also for the integral of the conjugate multiplied by sin 
6, i.e. for its integral with respect to x instead of 6." 
The results are not without mathematical interest, but 
from the aerodynamic point of view the main point of 
interest is that such useful formulas are now known. 

The formulas for the velocity at the surface of an air- 
foil, especially Approximation II for airfoils designed 
from a knowledge of g, and g;, may be used to estimate 
rapidly the positions and magnitudes of the velocity 
maxima on the surface at any lift coefficient, and also 
to find the theoretical C,-ranges for low-drag airfoils, 
defined as the range of lift coefficients for which, on 
both the upper and lower surfaces, the nearest velocity 
maxima to the leading edge still occur at their designed 
positions at the design lift coefficient. As an example, 
we have the result that, with a) = 27, the C,-range of 
any one of the ‘roof-top’ airfoils of the type previ- 
ously described is 9.27(2p,)'s'", where, as before, 
p, is the radius of curvature at the leading edge, and 
s is the slope of the nondimensional velocity (i.e., s = 
dg,/dx) over the front of the airfoil. (Divide by 
m/ao + '/2 for other values of ap.) For a given radius 
of curvature at the trailing edge, given position of 
maximum velocity, and given value of the maximum 
velocity or maximum thickness, there is a value of s 
which makes the C,-range a maximum; for example, 
for a cusped airfoil of given maximum thickness, with 
maximum suction at mid-chord, this value is s = 0.982, 
where ¢ is the maximum thickness (not the half-thick- 
ness) as a fraction of the chord. When the trailing 
edge becomes blunter (i.e., when pz increases) the 
most favorable value of the slope (from this point of 
view) decreases, and the maximum C,-range also de- 
creases slightly. For cusped airfoils the C,-range is 
nearly proportional to ¢”*. 

I have ventured to summarize some of the results of 
the new approximate airfoil theory for two reasons: 
no account of the theory has yet been published, and 
the formulas form a basis for judgments, and subsequent 
calculations, in designing, for test, airfoils of a new and 
unorthodox type. But the theory, like all mathemat- 
ical theories, must be used in conjunction with good 
engineering judnment, controlled experiment, and, 
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LOW-DRAG AND SUCTION AIRFOILS 


eventually, full-scale trials. The need for judgment 
makes itself felt at once. Knowledge of the position 
of maximum velocity on an airfoil, for example, may 
be considered to provide information on the probable 
position of transition to turbulence in the boundary 
layer, since transition takes place soon in an adverse 
gradient. But because of our lack of knowledge of the 
processes leading to transition, care is needed in using 
this information. Certainly other factors in addition 
to the velocity distribution—waviness and roughness 
of the surface, turbulence in the air, Reynolds number 
—affect the position of transition even on a flat or 
convex surface; and the effect of the velocity distribu- 
tion, and its interaction with other factors, do not de- 


pend solely on the position of the maximum velocity. . 


If the velocity falls off steadily and not too slowly after 
the maximum, then in practice, at high Reynolds 
numbers, transition will not be delayed to any appre- 
ciable extent beyond the maximum. But whether 
transition occurs before the maximum will depend on 
the velocity gradient and the other factors mentioned; 
and moreover transition may occur well after the 
velocity maximum in certain circumstances—for ex- 
ample, with a good surface and low turbulence in the 
air, if the maximum is followed by a small fall in veloc- 
ity and the velocity then begins to rise again, or even if 
the fall in velocity after the maximum is slow enough. 
The former state of affairs applies near to the nose of a 
good many airfoils, including some N.A.C.A. airfoils, 
at certain lift coefficients. So for a judgment of the 
position of transition, knowledge of the position of the 
velocity maximum can be only of restricted use for a 
preliminary ‘‘sorting out.”’ 

All this has, of course, a bearing on the definition of 
the C,-range for low-drag airfoils. A definition of the 
“theoretical” range has been given. A second possible 
definition is the range of lift coefficients over which the 
drag stays low and practically constant when measured 
on models made as carefully as possible and tested in a 
low-turbulence wind tunnel at a given Reynolds num- 
ber, say 30 X 10°; this ‘‘wind-tunnel” definition is 
useful to those to whom a suitable tunnel and model- 
making facilities are available, but we had in Great 
Britain—and still have—no such tunnel. A _ third 
definition—the ‘‘practical’’ definition—is the most 
probable range of lift coefficients for low and nearly 
constant drag under practical conditions of manufac- 
ture, flight, and maintenance; in the words of the origi- 
nal report, ‘‘values according to this def~ition are the 
values we should all like to be able to give but none of 
us can.”’ The information in Section 3, however, has 
taken us some way since the report was written. 

One other remark should be made in connection with 
transition. Experimental evidence of the effect of 
velocity gradient on the tolerance that can be allowed 
for waviness of the surface is still scanty and not at all 
systematic; such evidence as we have indicates that, 
once a fair value of the gradient has been reached, any 
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further increase needs to be very large to make any 
practical difference to the tolerance. Thus, once the 
nondimensional gradient s is 0.1, for example, it is 
doubtful whether it would make any practical differ- 
ence to the tolerance if s were increased to 0.2; but if 
it is made very much bigger, for example multiplied by 
10, so that it becomes 1.0 instead of 0.1, then it seems 
that the tolerance on waviness may be increased, so we 
expect to be able to tolerate a larger waviness very 
near an airfoil nose than further back. Also if s is much 
decreased the tolerance on waviness becomes less; 
but exact knowledge is still lacking. It is of some in- 
terest that the values of s for maximum theoretical 
C,-range are reasonable values from the point of view 
of tolerance on waviness. 

It is not only knowledge about the transition to tur- 
bulence in the boundary layer that we lack; the in- 
fluence of compressibility, when part of the flow field 
is supersonic, is also still imperfectly understood. For 
those transport aircraft for which it is still hoped to 
avoid compressibility effects, however, we still believe 
that it is a good thing to keep the maximum veloc- 
ity low enough, if possible, to ensure that the whole flow 
field is subsonic at top speed; if that is impossible, 
the best method of design is still a matter of doubt and 
controversy, and I do not propose to enter on it here, 
contenting mryself with the remark that once we have 
decided what we want, the formulas found help us to 
design to get it. One thing we often do want is thin 
wings with as good a maximum lift as possible, and 
there is some evidence, though not conclusive, that a 
large radius of curvature at the leading edge is beneficial. 

The formulas of approximate airfoil theory help us 
to design airfoil sections with required properties, once 
the requirements have been translated effectively into 
conditions on the velocity distribution at the surface 
and on the geometry of the section. We can design 
for a low maximum velocity for a given thickness, a 
large theoretical C,-range, a large leading-edge radius 
of curvature, a low moment coefficient at zero lift, 
and soon; we can see what is physically possible and 
are helped to make judgments about the inevitable 
compromises between different desirable characteristics; 
and we can design for a mixture of aerodynamic and 
geometrical properties, as structural and manufactur- 
ing questions and control characteristics may require. 
A few reports by Thwaites''~'* and Miss Douglas’ 
have been issued on these aspects from the Aerodynam- 
ics Division of the National Physical Laboratory. 

Finally, we remark that the formulas written are, 
after all, only approximate; they are flexible and gen- 
eral and very useful for helping us to get our ideas 
straight and to decide what is possible, and are accurate 
enough for most airfoils of small camber and not too 
great a thickness; but our mathematical apparatus is 
not complete unless we have, not only an exact method 
for determining the velocity on a given airfoil, but also 
some more exact theory for designing an airfoil with a 


Fic. 10. Streamlines in a boundary layer near a suction slot. 


----- Dividing streamline. S: Stagnation point. 


given velocity distribution. This gap has been filled 
by Lighthill.!° Let us suppose the airfoil section con- 
formally transformed into a circle, with the same undis- 
turbed relative air velocity, and all lengths represented 
as multiples of the radius of the circle. The airfoil plane 
is taken as the z-plane, with z = x + ‘zy, and the circle 
plane as the {-plane, the circle itself being ¢ = e”®, with 
6 = 0 corresponding to the trailing edge of the airfoil 
section (note that the meaning of @ is different from that 
in the approximate theory); the Kutta-Joukowski 
condition is assumed satisfied at the trailing edge. If 
Ya is the ratio, at incidence a measured from zero lift, of 
the relative air velocity at the point on the airfoil sur- 
face corresponding to the point 6 on the circle to the un- 
disturbed relative air velocity, and go the corresponding 
quantity at zero incidence, 


qa = go[cos (1/29 — a)/cos */26] 


so that go is known as a function of @ if g, is assumed. 
For symmetrical airfoils g, must be assumed on one 
surface only, since go is the same on both surfaces; for 
cambered airfoils, g, may be assumed at two different 
incidences for the two surfaces, subject to conditions 
of continuity on g. The resulting dirfoil may then be 
designed from the formulas 


x = S'(2/q) cos x d(cos 8) 
y = J'(2/q) sin x d(cos 6) 


for the coordinates, where x is the sum of the Fourier 
series in @ conjugate to the series for log go, so that 


x(0) =—(1/2x)P J, log go(t) cot 1/2(t — 


and may be found by methods previously described. 
Thus the airfoil is designed from a knowledge of qo(6), 
which is found if required from a knowledge of the veloc- 
ity distribution at incidence in terms of 6; since x 
varies rather like cos 0, we can make a good guess at 
the velocity distribution in terms of x. But in order 
to obtain the correct undisturbed relative air velocity, 
and to ensure that the airfoil.section shall be a closed 
curve, go(@) must be chosen so that 


q dé = 0 
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JS—,log go cos 6 dé = 0 
JS =z, jog sin 6 = 0 


The chord c, as a multiple of the radius of the circle, is 
rather less than 4, and with a little experience can be 
guessed fairly accurately, so the formula for the moment 
coefficient at zero lift, 


(Cu)o = —(4/c*) go sin 26 dé 


- may also be used to design from the beginning with the 


value of (Cyy)o in mind. 

This completes the present summary of the mathe- 
matical methods used. We use all our judgment, 
knowledge of experimental results, and experience to 
decide what we should like to design an airfoil section 

‘for in terms of its geometry and the velocity distribution 

on its surface; we use the approximate theory to sort 
out what is possible and what are the best compromises 
to make, and for trial designs; and we use the exact 
theory to finalize the design, particularly for airfoils 
that are thick or highly cambered. 


(5) THe DesIGN oF SUCTION AIRFOILS 


In considering the principles of the design of low- 
drag airfoils, we saw that we should like the velocity 
outside the boundary layer to rise to a position as far 
back along the wing as possible, but that we are hindered 
by the danger of turbulent separation if the rate of 
velocity decrease at the back of the airfoil is too great; 
the danger increases, roughly, when the thickness of 
the airfoil increases and when the position of maximum 
suction is moved further back. If, however the whole, 
of the pressure recovery, or velocity decrease, is con- 
centrated over a very narrow interval along the chord, 
over, which the boundary layer, or as much of it as 
necessary, is sucked away to stop separation, all danger 
of separation is avoided, and we can have, if we wish, 
a favorable velocity gradient over the whole of the rest 
of the chord. More specifically, if the upstream bound- 
ary layer, or part of it, is sucked in through a slot, there 
must be a streamline that divides the air crossing the 
slot from that entering it, and that streamline must 
meet the surface at a stagnation point (Fig. 10). With 
the stagnation point behind the slot the flow in the 
boundary layer along the surface is reversed in direc- 
tion between the stagnation point and the slot, and all 
danger of separation is avoided with a well-designed 
slot if the design of the airfoil is such that the pressure 
recovery takes place wholly between the upstream lip 
of the slot and the stagnation point. Up to now suction 
airfoils have actually been designed with a discontinu- 
ous drop in velocity; the velocity rises to a given posi- 
tion on the chord, drops discontinuously, and thereafter 
rises, stays constant, or at any rate does not decrease 
rapidly enough to produce any danger of separation, to 
the trailing edge; the velocity distribution has this 
character on both the top and bottom surfaces. 
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LOW-DRAG AND SUCTION AIRFOILS 


When Dr. Griffith suggested that airfoil surfaces 
should be designed specially for suction, the mathe- 
matical apparatus of the approximate theory was 
ready and was at once applied; the exact theory of 
design followed later. On the approximate theory the 
simplest thing to do is to take the first, linear, approxi- 
mation to g/U for a symmetrical airfoil to rise linearly 
fom 1+ aatx = 0tol + bat x = X, where it 
drops discontinuously to 1 + c, and thence proceeds 
linearly to 1 + dat x = 1 (Fig. 11); the ordinates of 
the airfoil are given by 


y = afo + + dfs + c)fs 


where fo, fi, and f. are the same functions as occurred 
in the design of ‘“‘roof-top”’ airfoils, and f; is easily cal- 
culated. 


The first suction airfoil, designed in this way and 
tested by Mr. Richards at the National Physical 
Laboratory, was about 16 per cent thick and had the 
velocity discontinuity and suction slot on each surface 
at seven-tenths of the chord from the leading edge. The 
shape of the airfoil section, and Approximations I and 
III to the velocity distribution at zero lift, are shown 
in Fig. 11. Transition to turbulence and boundary- 
layer separation were observed both by filaments of 
smoke introduced into the boundary layer and by fine 
silk threads, about '/, in. long, attached to the surface. 
The model as made and tested did not, in fact, repro- 
duce accurately the computed shape, but Mr. Richards 
soon found several interesting results, among which we 
may mention here that the method of design was satis- 
factory, that separation could be avoided by sucking 
away less than half the air in-the boundary layer, and 
that, in spite of the favorable velocity gradient behind 
the slot, the flow in the boundary layer became turbu- 
lent soon after the slot except at low Reynolds numbers. 
The airfoil surface is concave to the flow aft of the slot 
(Fig. 11), and there is no doubt that the transition to 
turbulence is caused by the instability of flow in a 
boundary layer over a concave surface; this instability 
had been studied theoretically by Gértler,'* and transi- 
tion on a concave surface of constant curvature was 
studied experimentally by Liepmann," of whose results 
we were informed just as the possibility of several other 
causes for the transition after the slot had been elimi- 
nated by special experiments. Liepmann’s results 
fitted ours as well as could be judged in view of the 
vatiable curvature of the surface of our airfoil. It is 
not possible to design an airfoil of this type without a 
concavity aft of the slot, but the concavity can be re- 
duced by increasing the velocity over the tail and the 
radius of curvature at the trailing edge; Richards and 
I soon convinced ourselves, however, that at Reynolds 
numbers of 20 millions or more it would be impossible 
to design an airfoil for which the gain due to increasing 
the length of the laminar portion of the boundary 
layer would not be more than counterbalanced by the 
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Fic. 11. Symmetrical suction airfoil, 16.2 per cent thick, and 


the theoretical velocities on its surface at zero incidence. 
Approximation I. Approximation III. 
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Fic. 12. Symmetrical suction airfoil, 30 per cent thick, and 
the theoretical velocities on its surface at zero incidence. 
Approximation I. Approximation III. 


increase of form drag. We therefore decided to accept 
turbulent flow in the boundary layer aft of the slot, 
and the next airfoil designed had the velocity discon- 
tinuity and the slot further back, at eight-tenths of the 
chord from the leading edge; and, since theoretically the 
idea should work equally well for very thick airfoils, the 
next airfoil, which was still symmetrical, was 30 per 
cent thick. It was still designed by the approximate 
theory, and the shape of the section, together with Ap- 
proximations I and III to the velocity distribution at 
zero lift, is shown in Fig. 12. It should be stated at once 
that Approximation III agrees quite well with experi- 
ment, so there is no reason to believe that it is insuf- 
ficiently accurate for a symmetrical airfoil even 30 per 
cent thick. A comparison of Approximations I and 
III in Fig. 12, however, shows that the favorable veloc- 
ity gradient of Approximation I has been replaced for 
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Fic. 13. Symmetrical suction airfoil, 40 per cent thick, with a Cz-range from —1.35 to +1.35, and the theoretical velocity distribu- 
tion on its upper surface at Cz = 1.35. 


Fic. 14. Cambered suction airfoil, 41 per cent thick, with a Cz- 
range from 0 to 2.85. 
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Fic. 15, Cambered suction airfoil, 311/2 per cent thick, with 
a velocity discontinuity on the upper surface only and a Cr-range 
from 0 to 2, and the theoretical velocities on its surface at C, = 0 
and at Cy = 2. 


some distance forward of the discontinuity by an ad- 
verse gradient according to the more accurate Approxi- 
mation ITI, and the design criteria have to this extent 
not been fulfilled; this adverse gradient could lead to 
laminar separation at low Reynolds numbers and pre- 
mature transition to turbulence at higher Reynolds 
numbers. This defect in the design could be removed 


by an alteration in the first approximation, 1 + g,, 
to g/U, but thick suction airfoils have, in fact, lately 
been designed by use of Lighthill’s exact theory. Some 
specimens are shown in Figs. 13-17, of which Fig. 15s 
due to Mr. M. B. Glauert,'* and the rest to Mr. Light- 
hill. 

The original 30 per cent thick airfoil had a C,-range 
from —0.6 to 0.6, and it is evident that much larger 
C,-ranges may be obtained. Fig. 13, for example, 
shows a 40 per cent thick symmetrical section with a 
C,-range from —1.35 to 1.35, the velocity distribution 
on the upper surface at C, =:1.35 being flat both in 
front of and behind the slot, asshown. In Figs. 14, 15, 
and 16 three cambered airfoil sections are exhibited; 
cambered sections have not previously been men- 
tioned. No cambered suction airfoils were, in fact, 
designed by the approximate theory, although for- 
mulas are known for the design of centerlines with dis- 
continuities in g,; it is clear from the approximate 
theory that not only may cambered suction airfoils be 
designed, but by introducing discontinuities in both 
g, and g; the discontinuities in the velocity on the top 
and bottom surfaces may be treated independently, 
and an airfoil may be designed with a discontinuity on 
one surface but not on the other. Fig. 14 shows a 41 
per cent thick cambered airfoil, with a velocity dis- 
continuity on each surface, and a C,-range from 0 to 
2.85; Fig. 15 a 31'/2 per cent thick cambered airfoil, 
with a velocity discontinuity on the upper surface only, 
and a C,-range from 0 to 2; in this case the theoretical 
velocity distributions at the ends of the C,-range are 
also shown. The theoretical value of (Cy,)o for this 
airfoil is zero. Figs. 16 and 17 illustrate, by a couple 
of rather extreme designs, some of the curious shapes 
that can theoretically be designed. Fig. 16 shows a 
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cambered 30 per cent thick section with a C,-range 
from 0 to 2.69; the velocity discontinuity on the lower 
surface is rather far forward, and there would presum- 
ably be a rather long turbulent boundary layer on the 
concave surface behind it; but the C,-range is remark- 
able for a 30 per cent thick section. Fig. 17 shows a 
theoretical symmetrical section 70 per cent thick, 
together with the velocity distribution on the upper 
surface at a lift coefficient of 3. 


Of all the sections shown in Figs. 13 to 17, only that 
in Fig. 15 has been made and tested in a wind tunnel 
at the National Physical Laboratory, and the tests so 
far carried out are preliminary. It is still impossible 
to say whether the design is satisfactory or not. Cer- 
tain troubles were encountered, but the shape of the 
slot entry was not entirely satisfactory since an attempt 
was made to design a slot suitable for blowing as well 
as for suction, and it now appears that that is not really 
possible, at any rate on this shape. When flow without 
separation was established, good agreement was ob- 
tained with the properties forecast in the theoretical 
design. Further tests will be carried out with a new 
and improved slot. An airfoil of this shape has also 
been tested in Australia, where fairly satisfactory results 
were obtained by using more than one slot and round- 
ing off the surface near the design position of the dis- 
continuity in velocity. Final judgment must await 
further experiment. 
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Fic. 16. Cambered suction airfoil, 30 per cent thick, with a 
*Cr-range from 0 to 2.69, and the theoretical velocities on its sur- 
face at C, = 2.69 and at Cr = 0. 
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Fic. 17. Symmetrical suction airfoil, 70 per cent thick, and the 
theoretical velocities on its upper surface at Cy = 3. 


(6) Some ExPERIMENTAL RESULTS ON SUCTION 
_AIRFOILS 


Although the shape of the model of the first suction 
airfoil designed, 16 per cent thick with the slot at seven- 
tenths of the chord from the leading edge, did not ac- 
curately agree with the computed shape, it is of interest 
to reproduce at least one set of experimental results— 
namely, the measured velocity distribution at zero lift 
with and without suction.” These are shown, and 
compared with the theoretical distribution, in Fig. 18. 
The discrepancies between theory and experiment are 
largely due to the errors in shape. The results without 
suction are of considerable interest, and were closer to 
the theoretical values than we had expected. The 
boundary layer, which separates from the surface, is 
shown by these results to rejoin the surface at about 
0.85 of the chord from the leading edge, and this ac- 
counts for the high effectiveness found for control 
without suction in wind-tunnel experiments, both in 
two-dimensional tests at the National Physical Labora- 
tory and in tests in a wind tunnel at the Royal Aircraft 
Establishment on a model of a complete aircraft fitted 
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Fic. 18. Calculated and measured velocities on a symm 
x: Values measured with suction. 


Calculated values (Approximation III). 
ber = 3.85 X 108. 


og x 


etrical suction airfoil, 16.2 per cent thick (Fig. 11), at zero lift.. —— 
© Values measured without suction. 
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Fic. 19. Calculated and measured velocities on a symmetrical suction airfoil, 30 per cent thick (Fig. 12), at C, = 0 and C, = 
, upper surface;* 


0.883 


= 0.883, lower surface; = 0. 


——Calculated values (Approximation III). 
Values measured without suction (Cr = 0). 


x, A.© Values measured with suction. x: 


CL = 
Reynolds number = 0.96 108. 
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with wings of the shape shown in Fig. 11. The maxi- 
mum lift and aileron effectiveness were, in fact, larger 
with these wings without suction than with high-speed 
wings of a more orthodox shape which tapered from 12 
to 10 per cent thick. It should be remarked that re- 
sults without suction are liable to a large, and, on the 
whole, probably favorable, scale effect. 


The ‘‘measured”’ velocities were obtained, from pres- 
sures measured at the surface, by assuming the pres- 
sure constant across any section of the boundary layer 
and using Bernoulli’s equation, so that normally the 
velocity should be that just outside the boundary layer. 
In the regions of separated boundary-layer flow, how- 
ever, the results should perhaps be regarded purely as a 
measure of the pressure on the surface. 

Some of the experimental velocities*® measured 
at the surface of the 30 per cent thick symmetrical air- 
foil of Fig. 12 are compared with theoretical results in 
Fig. 19, where the full lines are again the theoretical 
results and the points the experimental results at lift 
coefficients zero (actually the mean of results on the 
upper and’ lower surface at C, = —0.02) and 0.883; 
the dotted curve refers to the results without suction 
at zero lift. 

The theoretical slope of the curve of lift coefficient 
against incidence for this airfoil is 7.9; experimentally, 
with sufficient suction to prevent separation, a straight 
line was obtained up to the highest lift coefficient, 2.3, 
measured, and the slope was 7.3; there was no sign 
of a stall. Thus these thick suction airfoils have large 
lift-incidence slopes, rather near the theoretical values; 
if sufficient suction is available to prevent separation, 
the stall is at such a high lift and incidence that the 
maximum lift for these airfoils probably has no practi- 
cal significance. Moreover parallel straight lines 
were obtained for the (C,, a) curves with a 20 per cent 
plain flap deflected, the increment in lift coefficient 
being about 0.045 per degree of flap deflection; meas- 
urements were made with deflections up to 14°. 

With sufficient suction to prevent separation, the 
graphs of the moment coefficient and hinge-moment 
coefficient against lift coefficient are also all represented 
fairly well by straight lines, and are in reasonable agree- 
ment with theory. The large rate of change of hinge- 
moment coefficient with flap deflection may necessitate 
soine special system of control balance, even with power 
operation. 

The lift coefficients without suction, shown in Fig. 
20 for a Reynolds number of 2.88 millions, are of con- 
siderable interest. It should be remembered, however, 
that results without suction are all subject to consider- 
able scale effect, probably favorable as the Reynolds 
number increases. This is illustrated in the present 


case by ‘the dotted curve in Fig. 20, which represents 
results at zero flap deflection at a Reynolds number of 
0.96 millions; there is no reason to believe that scale 
effect will be absent above 3 millions, and tests at higher 


° 


6° 10° 16° 20° 
weadence 
Fic. 20.° Measured lift coefficients, without suction, of a 
symmetrical suction airfoil, 30 per cent thick (Fig. 12), with a 
20 per cent plain flap. ——Reynolds number 2.88 X 10°. 
----- Reynolds number 0.96 X 108, zero flap deflection. 


Reynolds numbers are required. The lift curves without 
suction also certainly depend on the design; for the air- 
foil of Fig. 15 a lift coefficient greater than 1 is ob- 
tained without any bends in the curve or traces of a par- 
tial stall. 

Without suction, at a Reynolds number of 2.88 
millions, the moment coefficient changes linearly and 
agrees reasonably well with theoretical values derived 
for completely potential flow up to an incidence between 
6° and 7°, and then changes violently; this change 
appears to be due to the replacement of a gentle laminar 
separation by a violent turbulent one, with which the 
dips in the curves in Fig. 23 are also associated. The 
measured values of the hinge moments at a Reynolds 
number of 2.88 millions indicate overbalance following 
the response to control movement. However, scale 
effect on results without suction makes it impossible 
to draw firm conclusions. 

One other rather curious result without suction should 
be mentioned: at a Reynolds number of 2.88 millions 
the drag coefficient at small incidences and small flap - 
deflections is only about 0.010, which is a very low value 
for such a thick airfoil with such an unorthodox shape. 

We may sum up our results so far by saying that 
with sufficient suction to prevent separation, the be- 
havior in the wind tunnel is in accordance with theoreti- 
cal ideas; without suction the results are better than 
might have been expected. I might perhaps add that 
as a result of wind-tunnel tests on the 31'/2 per cent 
thick airfoil with a single slot, shown in Fig. 15, the 
Australians believe that in the event of suction failure 
lateral control, although poor, will be tolerable and an 
emergency landing should be possible. 

In this summary of experimental results no mention 
has yet been made of the quant.ty of air that must be 
sucked to prevent separation or of the effective drag; 
some discussion is advisable before we start to look at 
experimental results on quantity and drag. 
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Fic. 21. Streamlines in a boundary layer near a suction slot. 


Dividing streamline. S: Stagnation point. 


(7) SucTION QUANTITIES 


We have already remarked that separation was 
avoided on the first suction airfoil tested by sucking less 
than half the air in the boundary layer just upstream 
of the slot; this result was explained, and a means sug- 
gested for computing the minimum quantity necessary, 
by a simple calculation due to Sir Geoffrey Taylor. 
Let us assume that the total head does not change ap- 
preciably along a streamline in the boundary layer as it 
crosses the slot from section 1 to section 2 in Fig. 21. 
If 1%, wu: are the velocities and /;, p2 the pressures at 
those sections on the same streamline in the boundary 
layer, then 


— = p2 — pi 


If we now further assume that the static pressure is 
constant across each of the sections 1 and 2 of the 
boundary layer, 


— Pi = — U;?) 


where U;, U2 are the velocities just outside the layer, 
so 


uy? — ug? = — U2? 


and if the fluid in a streamtube is to have sufficient total 
head to cross the pressure discontinuity, “,? must be 
greater than U,”? — U,”._ Hence the portion of the up- 
stream boundary layer extending from the surface to 
the streamline where 


must be sucked away. If we assume that U2/U,; has 
the theoretical value in a potential flow, and that the 
velocity distribution in the boundary layer upstream 
of the discontinuity can be found sufficiently accu- 
rately, the minimum quantity that must be sucked can 
be calculated. Moreover, by using the equation of 
continuity for a streamtube we may compute, down- 
stream of the discontinuity, the distance from the sur- 
face of any streamline that has crossed the slot, and 
since the velocity, 12, on it is given by 


U2? = uy? — (U? 
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we can find the distribution of velocity in the boundary 
layer downstream of the slot, the momentum thickness 
there, and any other property of the boundary layer 
that we 

Gregory’ has tested the assumptions of the theory by 
boundary-layer explorations. He finds that the change 
of total head along a streamline in the boundary layer 
crossing the slot is small compared with the changes 
normal to the streamline, so Taylor’s assumption is 
justified. Near a slot, however, the surface of a suc- 
tion airfoil and the streamlines in the boundary layer 
have a large curvature, which, for some of the stream- 
lines, is enhanced by the sink effect; this curvature 
produces a large pressure gradient across a section of 
the boundary layer. Consequently if the sections 1 
and 2 of Fig. 21 are taken just upstream and just down- 
stream of the slot, corrections must be applied to the 
simple theory above to allow for the effects of curvature. 
One such effect is that on any particular streamline » 
will be less than the value given above, and on this 
account the minimum suction quantity would be 
greater than calculated. Another effect is that U; is 
less and U; is greater than the values that would be 
calculated from Bernoulli’s equation and the measured 
pressures on the surface with the assumption that the 
pressure is constant across a section of the boundary 
layer; these latter values are the ‘experimental veloci- 
ties at the surface’ with which theoretical results from 
potential theory were in good agreement. Hence the 
second effect makes U; less and U2 greater than the 
theoretical values, and on this account the minimum 
suction quantity would be less than calculated. The 
two corrections nearly cancel, leaving the simple cal- 
culation nearly right—as we should expect since the 
large curvatures are localized near the slot. 

Note also that the curvature of the streamlines in the 
boundary layer explains why the ‘‘measured’”’ discon- 
tinuity, is so near the theoretical one, instead of being 
softer, as we might have expected. 

When Taylor’s theory is applied, calculations are 
necessary of the velocity distribution in the boundary 
layer; for laminar layers these are not difficult, but for 
turbulent layers are still estimates rather than satis- 
factory calculations, but are probably sufficiently nearly 
correct for our purposes. For laminar layers the calcu- 
lated minimum quantities agree with experiment, but 
for turbulent layers separation can be avoided with 
rather less suction than the theoretical minimum, 
probably because of turbulent mixing. A few results 
are shown in Fig. 22, where the nondimensional coef- 
ficient Cg, equal to the quantity sucked per unit time 
per unit span on one surface divided by the product of 
the chord and the undisturbed relative air velocity, is 
plotted against the position of transition to turbulence 
in the boundary layer. The results refer to the airfoil 
of Fig. 12 at zero incidence. The increase of Cg as 
transition moves forward is not the same for all airfoils; 
it depends particularly on the ratio, U2/U,, of the veloc- 
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ities outside the boundary layer at the discontinuity, 
being less rapid for larger values of this ratio, as, for 
example, on the 16 per cent thick airfoil of Fig. 11. 

The necessary CQ is found in experiment to vary very 
little with C, when the lift is obtained by deflection of a 
fap; when the incidence is altered, the necessary Cg 
with a largely laminar boundary layer does not alter 
within the low-drag C,-range, but at higher lifts the 
value for the upper surface increases because transition 
moves forward, and, in addition, increasingly adverse 
velocity gradients produce a thickening of the turbulent 
boundary layer. However, airfoils may be designed 
with very large C,-ranges; moreover in flight the 
speed falls as the lift coefficient increases, and the total 
quantity that must be sucked does not alter so rapidly, 
and is reasonable even at the largest incidencés that are 
likely to be required. 

The quantities considered so far are those needed to 
prevent boundary-layer separation; if separation has 
taken place, larger quantities—which it is difficult 
to specify exactly, but which may in certain circum- 
stances be as much as double—are required to make 
the boundary layer return to the surface. 


There is some evidence that, with slot widths greater 
than the boundary-layer thickness, the necessary Cg 
increases gradually with increase in width; the losses 
in, the duct entries decrease, however. 

Finally we remark that though it is certainly neces- 
sary to have as clear an idea as possible of the minimum 
quantity that must be removed to prevent separation, 
it does not follow that that quantity will be used in 
practice. As we shall see, that quantity may not even 
necessarily give the least effective drag; the quantity 
that does that will depend on the losses in the ducts. 
Even that, however, does not settle the matter. We 
must consider what we are going to do with the air 
sucked away from the boundary layer. We may either 
use a separate pump, somewhat in the nature of a ducted 
fan, to remove the boundary-layer air, or take the air 
through the wing slots to the intake of a turbojet engine; 
and the arrangement to be preferred must be decided 
separately for each airplane design. Certainly with 
the latter arrangement the amount of air to be with- 
drawn will be settled by the requirements of the engine, 
rather than by purely aerodynamic considerations, so 
long as the losses in the ducting with that amount of 
air are not prohibitively large compared with those in 
other intakes. 


(8) Dracs 


Drags of suction airfoils are measured in wind tunnels 
by the pitot-traverse method; when we consider what 
is measured in that way we see that it is nearly the 
drag of the tail of the airfoil aft of the slots. For essen- 


tially what we measure is the loss of momentum in the 
wake, and the parts of the boundary layers that have 
lost most momentum before the slots are drawn into 
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Fic. 22. Minimum suction quantities to prevent separation 
on a symmetrical suction airfoil, 30 per cent thick (Fig. 12), at 
zero lift: theoretical and experimental values at two Reynolds 
numbers. Cg, equal to the quantity sucked per unit time per 
unit span on one surface, divided by the product of the chord 
length and the undisturbed air velocity, is plotted against the 
position of transition to turbulence in the boundary layer, meas- 
ured along the chord. For the points ringed the values are high 
because of the tendency to laminar separation at the lower Rey- 
nolds number produced by the adverse velocity gradient in 
front of the slot (Figs. 12 and 19). 


the slots and do not appear in the wake with this loss 
of momentum. 

As we have seen, the momentum thickness of a 
boundary layer just aft of a slot may be estimated, and 
thence the value of the momentum thickness at the 
trailing edge and the profile drag may be computed in 
the usual way.”* Gregory finds that when the theoreti- 
cal minimum quantity is sucked, the momentum thick- 
ness of a laminar boundary layer is increased by only 
about a quarter across the slot; the increase across the 
slot of the momentum thickness of a turbulent boundary 
layer varies appreciably with the ratio, U2/U;, of the 
velocities outside the boundary layer on the two sides 
of the discontinuity; for the airfoil of Fig. 12 the mo- 
mentum thickness of a turbulent boundary layer may 
be multiplied by a factor varying from 4.5 to 6 for 
Reynolds numbers between 10° and 10%. The increase 
of momentum thickness falls rapidly as the quantity 
sucked is increased. 

In the calculations of profile drag the boundary layer 
was assumed to be turbulent along the whole of the sur- 
face of the airfoil behind the slots. In the experiments 
on the airfoil of Fig. 12 (at Reynolds numbers of 2.88 
X 10° and 0.96 X 10°) there were some lengths of 
laminar layer behind the slots before the transitions 
to turbulence. We must expect the computed values 
of the profile drag to be higher than the measured values 
on this account, The comparisons of theory and experi- 
ment for the profile-drag coefficients of this airfoil are 
shown in Fig. 23. 
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Fic, 23. Profile-drag coefficients of a symmetrical suction airfoil, 30 per cent thick (Fig. 12), at zero lift: theoretical and experi- 
mental values at two Reynolds numbers. 


Some theory is necessary for extrapolation to other 
Reynolds numbers; calculated values of the profile drag 
of this airfoil (for the theoretical minimum suction 
quantities) are plotted against the position of transition 
at Reynolds numbers of 1, 10, and 100 millions as the 


dotted curves in Fig. 24; increasing the quantity sucked — 


would produce a rapid decrease of the profile drag when 
transition is far forward. 


To obtain the “effective” drag of the wing we have 
to take into consideration the power input to the 
pump, but one remark may be made before considering 
the equivalent “pump” drag. With sufficient suction 
the profile drag is much less than that of a normal wing, 
and the total drag that has to be balanced by a jet 
reaction (profile drag, induced drag, and parasitic drag) 
will also be considerably less; the requisite jet momen- 
tum is correspondingly reduced, and a turbojet may be 
profitably used at lower speeds. 


For a calculation of the equivalent pump drag we may 
suppose that the air taken in is discharged with the 
undisturbed relative air velocity and static pressure; 
the pump must then restore the loss of total head suf- 
fered by the air removed from the boundary layer at 
the slot, and must also make good the losses of head in 
the slot entry and the ducting. These losses, in any 
practical case, are still unknown. The ideal pump 
drag is found by neglecting them. For,a rough calcula- 
tion we may be content to neglect the effect on the 
static pressure of the curvature of the streamlines in the 


boundary layer; then if the boundary layer to a distance 
Y from the surface is sucked into the slot, the loss of 
head per unit time per unit span which has to be made 
good by the pump is ideally 


1 Mm uy 2. 
“J 
2U; 1 U; dy 
and if m is the efficiency of an assumed propulsive unit, 
no the efficiency of the pump, U» the undisturbed air 


velocity, and c the chord,. the ideal equivalent 
‘“‘pump-drag”’ coefficient is 


(Co) = c/o U; U; dy 

m/n2 is unknown; it may, in practice, easily be greater 
than 1. If we assume that it is equal to 1, and add 
the profile drag to the ideal ‘“‘pump”’ drag, we obtain 
the “‘ideal’’ effective drag shown by the full-line curves 
in Fig. 24 and cross-plotted in Fig. 25. It should be 
emphasized that these curves refer to the theoretical 
minimum quantity to prevent separation, and not to the 
quantity to give minimum ideal effective drag, which 
would be rather meaningless in the absence of knowledge 
of duct losses. 


We can say nothing definite yet about-losses in prac- 
tical ducting systems; the losses will, however, depend 
on the quantity to be sucked and we do know some- 
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thing about that. We may therefore carry out hypo- 
thetical calculations with reasonable assumptions about 
quantity and various assumptions about the efficiency 
of the slot and ducting system, and arrive at the follow- 
ing rather general conclusions. With transition far 
forward and an inefficient ducting system the effective 
drags are too large for the scheme to be attractive; 
with transition far forward and an efficient ducting 
system, a suction airfoil 30 per cent thick may probably 
be designed whose drag is about the same as that of an 
18 per cent thick airfoil with transition in the same 
place. If laminar layers can be maintained to the slots, 
at 0.8c, the drag with an efficient ducting system is ex- 
ceedingly low, and even with an inefficient ducting sys- 
tem the quantity could be kept so low (and we should 
dearly keep the quantity low if we could.not have an 
eficient ducting system) that the effective drag would 
still be less than that of a thin low-drag wing with 
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Fic. 24. Calculated ideal ‘effective’ drag coefficients of a sym- 
metrical suction airfoil, 30 per cent thick (Fig. 12), at zero lift. 
----Profile-drag coefficients. Ideal ‘effective’ drag coef- 
ficients, equal to the sum of the profile-drag and ideal pump-drag 
coefficients, with losses in the slots and ducting not included. 
he coefficients are section coefficients, with no allowance for 
nduced drag; they are calculated for the theoretical minimum 
quantity to prevent separation, and not for the lowest drag, and 
or an assumed propulsive unit with an efficiency equal to that of 
he pump. 
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Fic. 25. Calculated ideal ‘effective’ drag coefficients of a sym- 


metrical airfoil, 30 per cent thick (Fig. 12), at zero lift. The 
results shown in Fig. 24 are cross-plotted against Reynolds num- 
ber. 


maximum suction at 0.6c, say. It is now apparent how 
important it is to keep as long a length of laminar layer 
as possible; the efficiency of the ducting system is im- 
portant if wé fail to maintain a sufficient length of 
laminar layer or if we wish to take a large amount of 
air for an engine through the wing slots. 

Two more remarks may be made about drag, one 
theoretical and the other practical. Although the 
main purpose of the scheme is to allow a maximum suc- 
tion far back without danger of separation, we may 
note that ideally there is always some gain, inde- 
pendently of the prevention of separation, if the loss of 
head in a boundary layer is restored by a pump instead 
of appearing as a loss of momentum in a wake. In the 
simplified case when the pump sucks the whole bound- 
ary layer from the trailing edge of a flat plate, it may 
be shown” that the ratio of the ideal drags, with equal 
efficiencies for the pump and propulsive system, is 
0/0, where # is the energy thickness of the boundary 
layer, 


/(u/U)[1 — (u/U)?]}dy 


and 6 the momentum thickness, 


— (u/U)}dy 


the integrals being taken across the boundary layer. 
Since u/U < 1,%< 6. Fora laminar boundary layer 
this saving is 21 per cent; for a turbulent one, at a 
Reynolds number of 10 millions, only about 8 per cent. 

We should notice also that even when we know the 
duct losses, and include them in the calculation of ef- 


e 
. 
4 


210 


fective drag, we still shall not have the whole story. 
For if the air is taken to the intake of a turbojet engine, 
we should consider the influence on the efficiency and 
specific fuel consumption of the engine. In fact, the 
full effect on performance can be judged only for a com- 
pletely designed aircraft. 


(9) S Lor anp DESIGN 


The story of the research with suction wings has been 
largely one of honest, painstaking endeavor, and it is 
right that it should be noted that just as the theory 
needed improvement, and is still capable of further 
refinement, so the experimental technique originally 
employed needed improvement, and may still be elabor- 
ated. 

One of the main difficulties was to design a ducting 
system that would lead to uniform suction along the 
span with a pump attached to one end. The design of 
the slots themselves, with particular reference to the 
angle between the surface and the slot axis and the 
ratio of the radii of curvature of the front and back lips 
to the width of the slot, had been studied by Fage and 
Sargent.** A glance at Fig. 12, however, will show that 
the direction of the slot relative to the surface is largely 
dictated by the airfoil shape and the need to accommo- 
date the ducting and collector without sharp bends in 
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the slots; in the airfoil of Fig. 12 the slot was, in 
fact, cut starting backwards towards the leading edge 
nearly parallel to the chord. Moreover, the slots used 
are at least as wide as the boundary layer, and prefer- 
ably wider, to keep down losses. 


Ducting design has been studied by Rawcliffe,” and 
two views (on different scales) of a recent ducting sys- 
tem are shown in Fig. 26. The air enters at AB, which 
is along the span, and is drawn off through EF. In the 
other view CD is along the boundary of a section of the 
airfoil through the chord. There should first be an 
entry carefully designed to avoid separation of the flow 
from the wall;* this is followed by a two-dimensionaj 
diffuser with straight walls at a slope of 1 in 12 until the 
arrow FP is reached, where the partitions start which are 
shown in the first view. The straight-walled diffuser 
is continued until an expansion ratio of 21/2 has been 
reached, and the walls are then curved to preserve a 
constant area. The trailing edges of the partitions 
should be sharp, cusped, and extended in a straight line 
a short way downstream—about the width of one bend 
if possible. Uniform suction is obtained with ducting 
of this type; the losses in it have been measured at 
present only in static tests, with still air drawn in, and 
were found to be just under one-fifth of the mean 
dynamic pressute at the slot throat. 
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Fic. 26. Experimental ducting design. 
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(10) Suction WINGs. CONCLUDING REMARKS 


In spite of the work that has been done, it is clear 
that much is yet unknown about suction wings. Most 
consideration has been given to thick wings, and though 
itis important to realize that suction wipgs may be de- 
signed of any thickness required, certainly including 
thicknesses less than 30 per cent, it seems likely that 
something like this thickness may be the first to be 
used, so as to take advantage of the greatly increased 
stowage that will be possible in the wing and the possi- 
bility of burying the power plant in it. The greatest 
gains will be obtained by applications to flying wings, 
and it becomes feasible to design a much smaller flying 
wing than would otherwise be the case. Moreover the 
possibility of designing for a very small or zero (Cy)o 
anda very large C,-range, witha maximum C, which may 
be regarded as practically unlimited in relation to pres- 
ent requirements, should ease the trim problem on 
flying wings, though suction failure still needs further 
thought. 

The large thickness must necessarily lead to a low 
critical Mach number, and the main application must 
be to transport aircraft; a straight 30 per cent thick 
wing at a lift coefficient of 0.5 would have a (theoretical) 
critical speed of 350 m.p.h. at stratospheric heights. 
But the wings would certainly be used with sweep- 
back, which might increase this critical speed to 450 
m.p.h.; it may be that suction will to some extent de- 
lay serious effects of compressibility; and somewhat 
thinner wings may be used on large craft. 

Our ignorance is now apparent. We have no experi- 
mental information at all about the application of the 
principle to finite wings with sweepback, or about the ef- 
fects of suction on such wings when shock waves are 
present. 

Our knowledge of the effects of compressibility on suc- 
tion wings is limited to the case when the whole field 
of flow is subsonic. The theoretical principles are then 
experimentally verified, but the suction quantity neces- 
sary to prevent separation, which remains less than the 
quantity of air in the boundary layer, may increase with 
speed. This effect may, however, have been due to the 
fact that on the airfoil tested there was a decreasing 
velocity for a short distance upstream of the slot simi- 
lar to, but less marked than, that shown in Fig. 12; 
the airfoil tested was 22 per cent thick, with the slot 
at three-quarters of the chord from the leading edge.” 
It is true that tests* on an airfoil of orthodox shape 
(NACA 0020) .had indicated that, provided efficient 
slots and ducting could be installed and extrapolation 
to higher Reynolds numbers was justified—in particu- 
lar, provided that boundary-layer separation was still 
associated with a shock wave at the Reynolds numbers 
of flight—then with suction the drag coefficient be- 
comes equal to its value without suction and at the 
critical Mach number, M), for a rise of drag at a Mach 
number 0.1 higher than M>, and this would be attained 


-on these airfoils at the slots. 
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with an economical use of power at such Reynolds num- 
bers. But this information is not sufficient for us to 
feel any certainty about the behavior of a specially 
designed suction wing at speeds above the critical. 

There will be many practical questions to be settled, 
and the field for speculation still appears boundless. It 
is an interesting reflection how every improvement in 
aerodynamic efficiency produces difficulties in other 
ways. In this case it is exemplified by the increase in 
severity of gust s‘resses, and the increased need for 
gust alleviation for passenger comfort, arising from the 
large increase in the lift-incidence slope. 

There are still unsolved problems associated with the 
aerodynamics of even a straight wing. The design of 
practical and efficient ducting and the correct design of 
wing tips and of a junction with any fuselage or nacelle 
that might be fitted are obvious examples. . Although 
some excellent systematic experimental work on slot 
design has been done,” further work is needed, es- 
pecially with the curvatures of the surface encountered 
Theoretically, it may be 
possible to extend the calculation of suitable shapes for 
slot walls—e.g. to design walls with constant pressure” 
—so as to take some account of the vorticity in the en- 
tering flow; and in the theoretical designs of the more 
extreme shapes it may be possible to take account not 
only of the sink effect (which is easy) but also of the 
shape of the slot entry (which is much harder). Above 
all, it must be remembered how much there i: to be 
gained if the boundary layer can be kept laminar to the 
slot; when no account is taken of duct losses, this is 
clear from Fig. 23; with duct losses included the values 
of the drag coefficient would fall even more sharply 
with backward movement of transition. The problem 
of maintaining laminar flow on these airfoils is much the 
same as on ordinary low-drag airfoils, which we have 
already discussed, except that we are now becoming 
even more ambitious and seeking to maintain laminar 
flow over greater stretches of the surface. What shall 
we do about the insect problem—about the turbulence 
caused by portions of flies’ bodies stuck on the surface 
near the leading edge? We must remember that the 
turbulence from any one protuberance spreads span- 
wise, and that one spanwise stretch of turbulent flow 
will require a largely increased suction quantity, not 
only along that stretch, but also (if we have taken pains 
to achieve uniform suction, as we must) all along the 
span; if we do not suck enough to prevent separation 
at the stretch of turbulent flow; the spanwise length of 
separated flow will spread enormously as it passes 
through the designed discontinuous increase of pres- 
sure. Shall we be able to devise chemical or other 
means of keeping the surface clean and free from flies, 
or must we be prepared to suck away a turbulent bound- 
ary layer through a slot not too far behind the leading 
edge? Shall we arrange to suck similarly behind access 
doors and the like? If it is too ambitious to seek to 
keep laminar flow right up to a single slot on each sur- 
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face, shall we not put in several siots along the chord? 
If we do that, could we not save even more by locating 
the rear slot nearer to the trailing edge than the 20 per 
cent of the chord shown in Fig. 12, for example? If we 
are toyise several slots on each surface, shall we use a 
surface designed with a single velocity discontinuity, 
as heretofore, or a more orthodox surface, or shall we 
not rather divide the velocity discontinuity into the 
same number of parts as there are to be slots, placing 
one at each? Perhaps in this way we can design an air- 
foil with large favorable velocity gradients everywhere 
except at the slots and at the extreme tail; but we 
have still to see whether we can avoid concavity of the 
surface upstream of the last slot, and to consider what 
the economy of the system will be with the divided 
velocity discontinuity; to these last questions we 
should soon have at least theoretical answers. 

Finally we note that the purpose of the original 
scheme, as shown in the designs of Figs. 11 and 12, 
may be considered achieved if we have a rising velocity 


gradient over most of the chord and some means of pre-. 


veriting separation over the rear portion of the airfoil, 
and that separation may be prevented by blowing air 
with sufficient momentum into the boundary layer 
through a slot in the surface. Some preliminary experi- 
ments?*!, *? with air blown into the boundary layer have 
been made with the airfoils of both Figs. 11 and 12, 
and indicate that the quantities required to prevent 
separation are larger with blowing than with suction, 
and that when separation is just prevented the form 
drag is larger with blowing; on the other hand, the 
form drag with blowing decréases continuously as the 
quantity increases, and with sufficiently large quantities 
becomes equal to and then less than that with suction. 
The actual recorded values, for what they are worth, 
for the ratios of the quantities required to prevent 
separation were 5:2 and 5:3 for the airfoils of Figs. 
11 and 12 respectively; the ratios of the quantities 
to give the same form drag were 7:1 and 3:1, respec- 
tively. The power required for blowing has not been 
estimated, but was probably greater than that re- 
quired for suction. The effect of blowing, however, 
depends not only on the quantity of air ejected, but also, 
on the velocity and pressure with which it leaves the 
slot; neither these nor the width of the slot have been 
systematically varied in experiments. Moreover, the 
experiments mentioned were made with ducting, slot 
shapes, and airfoil shapes, all specially designed for 
suction. There is nothing like the same background 
of theory and systematic experiment for blowing as 
there is for suction; and it seems unlikely, for instance, 
that the airfoil surface itself should have the same shape 
in the two cases, at any rate behind the slot. There- 
fore blowing cannot be condemned on the strength of 
these experiments alone, especially since, on the face of 
it, it looks an easier job to provide ducting for blowing 
than for suction. In view of all these considerations 
it does not seem worth while to consider further the 


1948 


recorded experimental results, except perhaps to men. 
tion that at. incidence the lifts were all just a little 
greater with blowing than with suction, as was to be 
expected. 


(11) Suction ArrForts DESIGNED FoR 
MAxXIMvuM LIFTs 


Another application of the same basic idea—to design 
the shape of the contour of an airfoil section specially 
for suction at a few points, where the design velocity 
on the surface decreases discontinuously, or at any 
rate very, very rapidly—is to attempt to design a thin 
airfoil, suitable for high speeds, with the principie ap- 
plied to increase the maximum lift. Generally speak- 
ing, a theoretical velocity distribution on the upper 
surface of ‘a thin airfoil at a high lift coefficient has a 
high peak near the leading edge, a sharp fall immedi- 
ately behind the peak, and a more gradual decline over 
the rest of the surface; the sharp fall behind the peak 
leads to separation and stalling. If this sharp fall is 
replaced by a discontinuity, with the boundary layer 
sucked away there, and if the adverse velocity gradient 
over the rest of the surface is reasonable, separation and 
stalling may perhaps be postponed to higher lift coef- 
ficients without any adverse effects on the other quali- 
ties of the airfoil section. 

It is not obvious, a priori, how successful the applica- 
tion of these ideas will be in practice. When separa- 
tion occurs at a high lift on a thin airfoil we expect 
only a partial stall, with the boundary layer rejoining 
the surface, so it is not clear how much gain is to be ex- 
pected. There is, of course, a subsidiary favorable ef- 
fect of suction, in that the sink effect itself produces a 
favorable alteratioh in the theoretical velocity dis- 
tribution; but we should not expect this effect to be 
large, because, especially in view of the thinness of the 
wings and the difficulty of fitting ducting, the quantity 
sucked must be kept small. 

Since the airfoil is thin, it may be thought at first that 
the approximate theory may be used for design. The 
discontinuity, however, must be located near the lead- 
ing edge, where the first, linear approximation fails 
completely, so Approximation II must be used. The 
most elementary theory of design must therefore be 
more difficult for these airfoils than for thicker airfoils 
with the discontinuity far back. Actually Approxima- 
tion II has not yet been used to explore the possibilities 
of design of these airfoils; it would be instructive if 
this were done. So far, several airfoils have been de- 
signed by Lighthill’® ** and one by Glauert*‘ by the 
use of Lighthill’s exact theory; Lighthill takes the 


sink effect into consideration and relies upon it to some 
extent to achieve the required results; Glauert does 
not. 

Small-scale models of Glauert’s design and of one of 
Lighthill’s designs have been tested** ® at the Na- 
tional Physical Laboratory at Reynolds numbers of 
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about half a million. Both designs were about the 
same thickness (8.65 and 8.6 per cent thick, respec- 
tively). Without suction both airfoils had a maximum 
normal-force coefficient Cy of 1.08. On Glauert’s 
design a normal-force coefficient of 1.25 is obtained 
with a small suction quantity, the suction coefficient 
Co being not greater than 0.003, and increase of Cg 
has only little effect on (Cy)max, On Lighthill’s design, 
(Cv)max. Seems to increase almost linearly with Co, 
up to the largest value used, 0.019 (much too large to 
be practical); interpolation in the results shows that a 
‘(Cy) max. Of 1.25 should be obtained for a Cg of 0.005. 

I have discussed these results only briefly because we 
do not know that the airfoils tested approach at all 
closely to the best that could be designed. No thorough 
attempt has yet been made to explore the theoretical 
possibilities and limits of design. 


(12) SUCTION 


No discussion of suction airfoils would be complete 


| without some mention of distributed suction through 


porous surfaces, though the experiments on this aspect 
are still in a rudimentary state. Theoretical consider- 
ations®® §” point to the need of very small quantities to 
prevent boundary-layer separation under adverse pres- 
‘sure gradients, and if the practical difficulties of the 
mechanical weakness of porous sheet, of shaping the 
sheet, and of the danger of clogging by rain and dirt, 
can be overcome, the applications seem very, promising, 
particularly to increasing the maximum lifts of thin 
airfoils by distributed suction through porous noses 
fitted to them.* 

The size and spacing of the pores in the surface must 
be small compared with the thickness of the boundary 
layer, and the distribution of the pores should prefer- 
ably be a random one. The aim is that the concept of 
a velocity through the surface and continuous all over 
it—not spotty—should be realized as closely as possible. 
This velocity, even if small, produces a marked de- 
crease in the rate of growth of the boundary layer in an 
adverse pressure gradient and marked changes in the 
distribution of the forward velocity across it; and it 
is this effect that delays separation. 

Here, too, the nose may be shaped to give the best 
results, a large radius of curvature at the leading edge 
being desirable, but the shapes are much more ortho- 
dox than those designed for suction at a single slot. 

The basic idea of distributed suction is not new, but 
it is only lately, since suitable materials*—sintered 
metals, e.g. sintered bronze—became available, that 
systematic calculations have been made, and experi- 
ments became possible. 

Theoretically, distributed suction ‘should also help 
considerably to stabilize the boundary Jayer and pre- 
serve laminar flow. (For the case of zero pressure 
gradient see reference 40.) Experimental confirma- 
tion is still lacking, and it is still, inf particular, an open 
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question whether a porous material can be produced 
with a surface sufficiently smooth for this purpose. 


In some ways the use of a large number of slots on an 
airfoil of orthodox section, as in Pfenninger’s experi- 
ments,?® may be considered to be a half-way stage be- 
tween the use of a few slots and the use of distributed 
suction. But the underlying principles are so differ- 
ent in the three cases that this view should not be 
pressed too far. It seems that distributed suction, 
which works in quite a different way, requires a smaller 
quantity than suction through separate slots. The 
final decision on the merits of the various schemes must 
be left to aeronautical scientists, designers, and time. 


(13) ConcLusIon 


This lecture is in the nature of an interim report on 
work still in progress, and, without any conscious direc- 
tion on my part, the descriptions became less and less 
detailed and satisfactory as we moved on. Looking 
back, I believe that is as it should be, since the sub- 
jects mentioned became progressively less and less 
explored, the ideas—or the possibility of their applica- 


‘tion—newer and newer, until we reached a stage when 


we were considering matters on which experiments 
have only just begun. 


I believe that the work on suction airfoils is as neat 
an example as we could wish to find of the amount of 
work to which a comparatively simple new idea can 
give rise, and of how the problems that arise should be 
tackled by all the methods—theoretical and experi- 
mental—at our command. Mathematical theories 
can be used to design both the models for test and the 
experiments themselves, and to help to interpret the 
results; the laboratory results must themselves be 
checked by full-scale trials. 


One final word: the shapes of suction airfoils doubt- 
less appear strange and novel to designers. I venture 
to prophesy that if future experiments show that con- 
siderable savings may be obtained with safety by their 
use, the time will come when they will ‘look right’’ to 
designers, whose intuition, based on a deep scientific 
background used without conscious effort, will then 
lead them to choose such shapes for the purposes for 
which they were produced. 
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Ira H. Abbott, Aeronautical Consultant, N.A.C.A.: 
It is a privilege to congratuate Dr. Goldstein on 
his excellent presentation of the research being 
done in Great Britain on low-drag and suction airfoils. 
The subject of today’s lecture appears particularly 
appropriate to me because today is the tenth anni- 
versary of the first Wright Brothers Lecture, de- 
livered by B. Melville Jones on the subject, ‘‘Flight 
Experiments on the Boundary Layer.’ The results 
presented 10 years ago were heartening to those of us 
in the N.A.C.A. who were then just getting started on 
the problem of extending the laminar layer on wings. 
Much has been accomplished in the past 10 years, 
but, as pointed out by Dr. Goldstein, much remains 
tobe done. The problem of improving the efficiency 
of wings offers a continuous challenge that cannot be 
met adequately at subsonic speeds until answers are 
obtained to the questions raised by Dr. Goldstein. 
Even at supersonic speeds, it may safely be pre- 
dicted that optimum wing efficiency will require that 
more attention be directed to the boundary layer. 

It is difficult to discuss this lecture because I am so 
thoroughly in agreement with the basic ideas presented. 
However, it appears appropriate to point out the paral- 
lel between the work in Great Britain and that of the 
N.A.C.A. in this country. In both instances, the re- 
search has been directed toward the broad objective 
of improving the efficiency of wings by increasing the 
extent of the laminar boundary layer. Although the 
emphasis has been different, this objective has been 
sought, in both instances, through the design of low- 
drag airfoils and the use of boundary-layer control. 
To a large extent, the results of the different approaches 
to the objective have supplemented each other and il- 
lustrate the advantages of separate attacks on an ob- 
jective by independent research teams. It is easy to 
overemphasize the differences in the two approaches, 
but at the risk of doing so, I wish to discuss some of 
them briefly. 

The most notable difference in the two approaches is 
the emphasis placed on theoretical and experimental 
work. Dr. Goldstein and his associates have tended 
to rely mainly on theoretical results with spot checks in 
flight or in the wind tunnel. When rigorous theory was 
lacking, they devised approximate or empirical rela- 
tions to fill the gap. - The success of this approach is 
well shown by today’s lecture. 

We, in this country, have been unwilling to rely on 
calculation when means could be devised for the ex- 
perimental determination of overall airfoil characteris- 
tics. Knowledge of the boundary layer has not pro- 


gressed sufficiently as yet to permit reliable calculation 


of all airfoil characteristics over the useful range of 


angle of attack. The N.A.C.A. two-dimensional low- 
turbulence wind tunnels were designed and built as 
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Discussions of the Lecture 


the result of an early decision that practical airfoil 
sections could be obtained in a reasonable length of 
time only if facilities were available to permit system- 
atic testing under conditions closely approximating 
those of flight. Initial exploratory investigations of a 
wide variety of airfoil sections served to outline the 
possibilities. The objective was then established of a 
series of sections that would exhibit low-drag charac- 
teristics with smooth surfaces without appreciable 
sacrifice of any of the good characteristics of the more 
conventional sections when compared under conditions 
of equal thickness ratios, camber, and surface condi- 
tion. The N.A.C.A. 6-series airfoils are the result of 
our attempt to reach this objective. 

The solution presented by Dr. Goldstein to the so- 
called inverse problem of designing airfoil sections to 
have desired types of pressure distribution differs 
markedly from the method used by the N.A.C.A. It 
should be noted that there are difficulties in an ade- 


‘quate mathematical formulation of the problem. The 


designer is faced, nevertheless, with the real problem of 
accurately associating profile shapes with desired types 
of pressure distribution. The difficulty of this problem 
presented a serious obstacle to the design of low-drag 
airfoils. The N.A.C.A. 1-series airfoils were obtained 
through empirical methods. The N.A.C.A. 2- through 
5-series airfoils were designed by various approximate 
methods. Our experiences with: these approximate 
methods were so unpromising that we abandoned this 
line of attack before a satisfactory method was de- 
veloped and concentrated on the development of an 
exact method based on the Theodorsen solution. 

The remarks of Dr. Goldstein with regard to the dif- 
ficulties involved in accurate computation of the Theo- 
dorsen solution are well justified. Even though the 
method is theoretically exact, an accurate numerical 
solution by usual methods may require’ the use of 80 
points and not less than two approximations. It was 
found, however, that the epsilon-curve could be modi- 
fied to produce changes of the desired type in the pres- 
sure distribution and that the resulting profile shapes 
and pressure distributions could be accurately related 
without undue difficulty in the computations. Con- 
tributions to this method were made by Theodorsen, 
Garrick, Jacobs, von Doenhoff, Pinkerton, and Naiman. 
This method has been used successfully in a routine 
manner for the design of the N.A.C.A. 6-series air- 
foils. 

It would be presumptious to compare the Goldstein 
and N.A.C.A. methods without experience in the de- 
tailed use of both of them. However, some general 
comments may be offered. The N.A.C.A. method, in 
practice, has been applied to the design of symmetrical 
sections or thickness distributions, although it is not 
limited to this application. These thickness distribu- 
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tions are then combined with mean lines derived by 
thin airfoil theory. The pressure distributions of the 
final sections may be obtained by Theodorsen’s solu- 
tion or, approximately, by superposition of velocities 
using tabulated characteristics. This method has many 
advantages because the pressure distributions of a 
large variety of sections may be obtained quickly from 
a modest number of basic characteristics tabulated by 
the N.A.C.A. 

Although the N.A.C.A. method is sufficiently accur- 
ate for most engineering purposes, no great claim for ac- 
curacy can be made for the first approximation pres- 
sure distributions of the cambered sections. This situ- 
ation is satisfactory for sections of small and moderate 
thickness ratios, but we are finding it necessary to de- 
sign thick sections directly as cambered airfoils. The 
Goldstein method has been applied directly to the cam- 
bered sections and contains within it an approximation 
of high order for the effect of camber changes. It is 
thus apparent that both methods have certain advan- 
tages. It is hoped that a way may be found to com- 
bine the advantages of both methods without intro- 
ducing undue complication. 

The research reported by Dr. Goldstein on specially 
designed suction wings is of considerable fundamental 
importance. Such wings avoid the difficulty of pro- 
ducing long favorable pressure gradients without exces- 
sive pressure recoveries in the limited region between 
the point of minimum pressure and the trailing edge. 
It is well known that the highest values of the lift- 
drag ratio are obtained with the largest extent of lam- 
inar flow and the largest amounts of camber consistent 
with prevention of turbulent separation. With ordi- 
nary airfoils, turbulent separation imposes limits on 
these design variables and on the thickness ratio. 
With suction airfoils, these limitations are avoided, 
making it possible to design thick sections that should 
have extremely high lift-drag ratios if extensive laminar 
layers can be realized. 

As indicated by Dr. Goldstein, the advantages of the 
suction airfoils are dependent to a great extent on ac- 
tually achieving the design extent of the laminar flow. 
This problem is not solved, as is apparent from the cited 
difficulties with insects in flight. Large numbers of 
airplanes are using low-drag wings. Special efforts 
have been made on some of these airplanes to obtain 
sufficiently smooth and fair surfaces to permit reason- 
ably extensive laminar flow. So far as I know, none 
of these airplanes are obtaining extensive laminar flow 
under actual conditions of field maintenance and service 
operation. It must be concluded that the important 
problem is to stabilize the laminar boundary layer to 
prevent transition resulting from the types of surface 
roughness and irregularities that cannot be avoided 
without unreasonable maintenance or impractical re- 
striction of operating conditions. 

The N.A.C.A. is now concentrating its efforts in this 
field on the problem of obtaining a stable laminar 


boundary layer- by means of boundary-layer control, 
Research conducted just before the war indicated that 
the use of suction slots would permit a greater extent 
of laminar layer to be realized. The results did not 
indicate that the use of slots stabilized the layer in the 
sense of preventing transition resulting from roughness, 
This research was discontinued during the war in favor 
of the development of airfoil sections of immediate 
practical value. 

In resuming this research we have been aided by the 
result of German and Swiss investigations. It ap- 
pears that the use of porous surfaces is of particular in- 
terest. Preliminary N.A.C.A. results indicate qualita- 
tively that suction through porous surfaces can stabilize 
the boundary layer, to some extent at least, to the ef- 
fects of roughness. This result is believed to be signifi- 
cant, although the studies have not progressed as yet 
to any practical surface or to flight values of the Rey- 
nolds Number. It thus appears premature to discuss 
the practical problems of using sections of the type pre- 
sented today. If it proves possible to stabilize the 
laminar boundary layer, we shall be forced to revise 
our concepts of what is and is not practical. 


Major General L. C. Craigie, U.S.A.F., Director of 
Research and Development: Professor Goldstein is 
to be congratulated on his masterly presentation of 
this important subject. The paper is timely, since 
it comes during a period when increased interest and 
activity in boundary-layer control are being evidenced 
with a view toward increasing the performance of mili- 
tary aircraft. 

Boundary-layer control has always shown great 
potential promise but has not, in general, come up to 
expectations because of the fact that it has been largely 
used as a so-called ‘‘aerodynamic crutch.’”’ In other 
words, it has been used as an atid in overcoming bound- 
ary-layer difficulties on aerodynamic surfaces that were 
designed primarily for operation without boundary- 
layer control. 

Professor Goldstein has presented a constructive new 
approach in designing the basic airfoil for operation with 
boundary-layer suction. 

_It is believed that this new concept will have promis- 
ing applications in’ military aviation. Continued re- 
search along the lines suggested by Professor Goldstein 
should lead to wing designs which would considerably 
increase the range of high-speed jet aircraft, since it 
should decrease the wing drag and make use of the 
Jarge amount of surplus air required by thermal jet 
propulsive systems. It should be emphasized that the 
advent of jet propulsion provides a built-in pumping 
system that eliminates the previous handicap of having 
to provide an extra pump for the purpose. 

Some of the unexplored applications that Professor 
Goldstein has mentioned appear especially promising 
for the future; for instance, the possibility of using, 
for high-speed flight, thin small wings that still have 
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adequate lift for take-off and landing under conditions 
that might be encountered in service. 

It is believed that boundary-layer suction will help 
diminish losses in the supersonic régime. In super- 
sonic flow, the shock waves occurring around the sur- 
face of an airfoil are associated with a rapid rise of 
pressure through the shock. At the surface of the 
wing this strong adverse pressure gradient causes an 
essential increase in the boundary-layer thickness with 
strong tendencies toward separation. It is» believed 
that systematic exploitation of Professor Goldstein’s 
methods will lead to considerable aerodynamic improve- 
ments in both the external and internal flow systems 
of supersonic aircraft and missiles. 

In conclusion, Professor Goldstein’s paper shows 
clearly the great number of theoretical and experimen- 
tal investigations which must be carried on and inte- 
grated to make progress on a single new concept. This 
emphasizes the complexity of a research and develop- 
ment program leading to improvements in air materiel. 

As a representative of the United States Air Force, I 
wish to thank Professor Goldstein for his stimulating 
paper. I,am happy to see this continued cooperation 
between Great Britain and the United States, which, I 
am sure we all agree, is to our mutual advantage. 


Captain Walter S. Diehl, U.S.N., R. & D. Group, 
Bureau of Aeronautics, Navy Department: Dr. 
Goldstein is an outstanding authority on mathematics 
in the fields of fluid mechanics and theoretical aero- 
dynamics. The Institute is most fortunate in having 


- such a distinguished guest speaker for the Eleventh 


Wright Brothers Lecture. Anything that Dr. Gold- 
stein might say regarding the work that he has been 
doing is certain to be of great interest to our members. 
It is particularly pleasing to find that he has chosen 
to give us a complete explanation of one of the most 
interesting applications of mathematical analysis to a 
specific problem. 

One of the most striking of the broad aspects of Dr. 
Goldstein's paper is the revelation of the extent to which 
modern mathematical analysis can enter into the detail 
design of airplane wings. This is not a matter of acci- 
dent, but rather it is the result of many years of thought 
and hard work on the part of those concerned. Our 
N.A.C.A. can point with pride to its efforts in this field. 
The development of highly satisfactory approximations 
for the exact theory as outlined in the lecture appears 
to have been a major achievement. For one thing it 
shows that the mathemati¢ian is a normal human being 
and that he does not always insist on doing a thing the 
hard way. It is significant to note the close coordina- 
tion of the theoretical studies with the experimental 
work. This, as many of us firmly believe, holds the 
key to rapid progress. In it, the theoretical analysis 
is always being crowded by the experimental data— 
which may not always agree with the predicted values. 


One well-known professor is reported to have told his 
class that, ‘‘When theory and facts do not agree, it is 
just too bad for the facts.” That was many years 
ago. The modern mathematician is not satisfied until 
it has been fully demonstrated that the theory accu- 
rately predicts the facts. 

The lecture really consists of two parts. The first 
part is concerned with Dr. Goldstein’s research work 
on laminar-flow wings, covering ground generally simi- 
Jar to previous N.A.C.A. research. Most of our mem- 
bers will find this an interesting account, and especially 
so in the sections devoted to flight evaluation of the 
practical difficulties encountered in the attempts to 
obtain and then to maintain laminar flow. 

The second section is concerned with suction wings 
and boundary-layer control. This section presents 
new material that requires some detailed technical 
comment. The comments that follow represent my 
personal appraisal and do not necessarily reflect the 
views of the U.S. Navy. 

Suction wings and boundary-layer control schemes 
have been with us for many years. The Bureau of 
Aeronautics sponsored an N.A.C.A. flight investiga- 
tion of the Katzmayr wing more than 20 years ago. 
There have been numerous attempts made _ subse- 
quently, with varying results, but boundary-layer con- 
trol must still be classed by the designer as an incom- 
pletely demonstrated device. . One finds three broad 
points of view among those who are familiar with the 
subject: 

(1) It is the most promising method we have for 
getting real improvement in airplane performance. 

(2) It shows some promise in special applications, 
but, in general, it is ruled out by considerations of the 
overall effects on the airplane. 

(3) It is just another “boot-strap”’ proposition. 

Most engineers hold to the middle ground point of 
view and probably will continue to do so until the mat- 
ter is settled for them in a convincing manner. 

Although Dr. Goldstein predicts cautiously that the 
time may come when suction airfoils will ‘look right” 
to designers, his paper is to be considered as a report 
on what has been done theoretically and experimentally 
rather than as a proposal to jump immediately into a 
practical application. Before this can be done there is 
another job for mathematical analysis—and it is an 
important job. An airplane consists of wings, power 
plant, and controls. Any practicable combination 
must have reasonably aeceptable stability, control, and 
performance characteristics. There can be no doubt 
that these characteristics will be affected by the use of 
suction wings. The mathematician must now deter- 
mine the boundaries for the design parameters giving at 
least some promise for a practical application of the 
suction wings. 

I should like to compliment Dr. Goldstein on his 
masterly presentation of a highly technical paper. By 
all ordinary, standards the complex technical ground 
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covered would seem to limit rather severely the size of 
audience that could appreciate the details involved. 
Dr. Goldstein has been able to’ surmount the barriers 
with ease and give us, not a technical broadside, but 
a most entertaining lecture. I am sure that all of those 


.who were fortunate enough. to hear the lecture join 


with me in thanking Dr. Goldstein for his contributions 
to our knowledge of boundary-layer flow. 


Professor William R. Sears, Director, Graduate 
School of Aeronautical Engineering, Cornell Uni- 
versity: First I must express my share of apprecia- 
tion to Professor Goldstein for this masterful lecture. 
It is particularly fortunate when the Wright Brothers 
Lecture, as has been the case today, can be devoted 
to the unveiling of an engineering invention. Dr. 
Griffith’s ingenious idea, based on designing an airfoil 
contour for a discontinuous velocity distribution, 
could hardly have been conceived by anyone unfamil- 
iar with the theory of aerodynamics. We have here a 
refutation of the popular fallacy that ingenuity and 
scientific learning are for some reason incompatible 
and that, therefore, inventions can be made only by 
untrained geniuses—diamonds in the rough—whose 
imaginations are not stifled by familiarity with the 
laws of physics or their mathematical expressions! 
In this respect today’s lecture is especially fitting to 
honor the Wright Brothers; they were certainly 
inventors of a more scientific character than is per- 
missible according to the popular legend! 

I am sure that the suction wings will be studied with 
great interest by the designers of flying-wing airplanes, 
as Dr. Goldstein has predicted. The combination of 
low drag, high lift, small pitching moment, and in- 
creased volume will be extremely attractive. On the 
other hand, in recent years we have seen a greater de- 
mand for thick profiles for conventional airplanes than 
for flying wings. This has been because of the require- 
ments of spar depth. Perhaps we shall see this trend 
carried on to produce conventional airplanes of very 
great aspect ratio and very thick profiles. If sweep- 
back up to about 40° is used—as one of Dr. Goldstein's 
figures indicates—we shall be thankful for the large 
values of dC,/da. 

One of the most startling results we have seen this 
afternoon is the reasonable behavior of the suction 
wings without suction. It seems possible that they 
can be designed to be safely usable without suction. 
This would certainly increase the attractiveness of the 
sections. But here again we would be faced with our 
ignorance regarding transition and turbulent separa- 
tion. 

Certainly Professor Goldstein’s remarks on this sub- 
ject and his whole lecture, for that matter, should in- 
spire us to further research on transition and the turbu- 
lent boundary layer. And since it is doubtful whether 


we shall ever understand the turbulent boundary layer 
before we fully understand simpler turbulent flows, 
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we may take this as an incentive also to attack the 
turbulence problem in general. . 

In this connection, I am interested in the speaker’s 
rather deflationary remarks on the subject of transi- 
tion. I think that in this country we have felt a bit 
more optimistic, believing that the amplified laminar- 
boundary-layer oscillations do represent the essential 
phenomenon leading to transition. But if Dr. Gdld- 
stein means that we cannot, at present, predict the 
transition’ point on the surface of an airfoil, and could 
not even if we knew the Reynolds Number of laminar 
instability, then we cannot refute him. Admittedly, 
we still have a lot to learn! 


One of the other discussers has already mentioned. 


certain parallel American work, especially on low-drag 
wings. On one essential point there seems to be agree- 
ment between Teddington and Langley Field—that is, 
the question, how smooth is a smooth surface? East- 
man Jacobs, the father of low-drag airfoils, whose in- 
fluence is certainly present in these low-drag and suc- 
tion airfoils, used to say that a wing was smooth when 
its finish resembled that of a new automobile. This 
put the criterion in terms of a product mass-produced 
and sold for a profit! The British experience seems 
to have been similar, and Dr. Goldstein’s remarks about 
the effects of weather and service are also encouraging. 
One might question his opinion that larger waviness is 
permissible near the leading edge because of the favor- 
able gradient there; probably this means “‘larger’’ 
only in relation to the boundary-layer thickness there 
and, therefore, very small permissible waviness indeed. 


The Lighthill theory for calculation of an airfoil to 
fit a prescribed pressure distribution also has its counter- 
parts in certain researches in this country. I think 
there is general agreement that this “‘inverse’”’ problem 
has no solution in general. Lighthill’s theory might 
appear to contradict this, until we notice that Lighthill 
specified the velocity distribution in the transformed 
“circle” plane, not on the airfoil. This can always be 
done, but it does not strictly constitute a solution of 
the inverse problem, since the exact variation of pres- 
sure along the airfoil remains unknown until the profile 
is calculated. 


I am not sure that in this country we have had any 
analog of Sir Geoffrey Taylor’s ingenious approximate 
calculation of the suction quantity. It is inevitable 
that Professor Taylor’s name should appear in a lecture 
like this. They say he has a whole desk drawer full of 
unpublished calculations similar to the one I refer to, 
and he gets them out when baffled engineers come to 
Cambridge for help, as they apparently often do. 


Finally, I cannot resist adding a word about this 
surprising aerodynamic-entomological problem. As I 
understand it, Professor Goldstein has told us this: that 
low-drag and suction wings are not different from most 
aeronautical inventions. They show great promise, 
but they still have a few bugs! 
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The discussions that follow were prepared for the 
West Coast presentation of the Wright Brothers Lecture. 


Philip A. Colman, Chief Aerodynamics Engineer, 
Lockheed Aircraft Corporation: One of the significant 
points brought out in Dr. Goldstein’s paper is the 
importance of airfoil smoothness. Many people be- 
lieve that, with boundary-layer control, the drag of a 
wing could be materially reduced, while, at the same 
time, the degree of accuracy of maintaining airfoil 
contour and the surface smoothness requirements 
could be relaxed. This is just the opposite to what is 
required according to this paper. If an extended 
region of laminar boundary-layer flow is not obtained, 
the drag reduction is small, and the ducting problems 
and internal power requirements are excessive. 

Lockheed has at various times experimented with 
boundary-layer control. The last investigation in- 
cluded wind-tunnel tests of a laminar-flow-type airfoil 
section of the kind used on the Constitution airplane. 
In the wind tunnel, slots were developed to obtain 
optimum shape for uniform suction and low intake 
losses. These tests showed a drag reduction, and 
therefore a flight-test program was initiated. A P-38 
airplane was equipped with a “‘glove,”’ which had 
sufficient width to obtain representative flow condi- 
tions. This ‘‘glove’’ had the same airfoil section 
tested in the wind tunnel and included the optimum 
slot. The source of suction was from the turbosuper- 
charger. These tests showed first the importance 
and the difficulty in obtaining true aerodynamic 
smoothness. With a smooth airfoil and boundary- 
layer control, the drag was reduced, but the net gain 
was not of sufficient magnitude to encourage further 
test work. Since this occurred at the time that the 
laminar-flow-type airfoil was first coming into use, the 
decision was that the development of this new-type 
airfoil reduced materially the amount to be gained 
from boundary-layer control. 

Now Dr. Goldstein shows us that to gain the ulti- 
mate in drag reduction, the airfoil shape should be 
designed to fit with the boundary-layer removal 
principle. This now appears to be obviously correct. 
I hope to see rapid development under the principle 
he has presented. This is a tedious job but should be 
attacked vigorously today so that the benefits can be 
realized in actual airplane design in the near future. 


A. M. O. Smith, El Segundo Plant, Douglas Air- 
craft Company, Inc.: For many years, boundary-layer 
control, especially stall control, in an indefinite sort 
of way has appeared to promise improvement in aero- 
dynamic design. However, until recently, no one 
had indicated a method of full utilization of the ef- 
fects that can be obtained by suction. Thanks to the 


recent progress in laminar-boundary-layer stability 
theory, the German porous plate theoretical studies, 


and Dr. Griffith’s idea of a velocity discontinuity at 
a slot, we are beginning to obtain promising methods 
and theory for use in further attacks on the drag 
problem. Dr. Goldstein is to be congratulated for 
his work on this problem and for his excellent report 
on the British efforts. 

During the past few years great strides in aircraft 
high speed have been made, but range capabilities have 
hardly progressed. Jet and gas-turbine aircraft at 
the present have less range than the best reciprocating 
engine types of recent years. Hence, a paper stimulat- 
ing thought toward improving the intrinsic efficiency 
of aircraft seems timely. 

I am pleased also to learn about the recent develop- 
ments in airfoil theory. Since the methods were de- 
veloped to be timesavers, we intend to study them care- 
fully. 

The thought of the very thick wings considered by 
Dr. Goldstein, especially ahead of the control surfaces, 
is disturbing to the industrial aerodynamicist, whose 
most difficult problem is to produce satisfactory flying 
qualities for the various critical conditions of flight. 
For a suction wing it is control suction power off. 
However, no evaluation can be made until more test 
data are available. 

Some additional information on the overall problem 
of the boundary-layer control aircraft may be of inter- 
est at this time. The significance of drag reductions 
by suction can only be determined by study of the over- 
all performance of such an aircraft, including allowance 
for the deleterious effects on engine performance. 
When the work of Dr. Pfenninger' on control of the 
laminar layer was received, such a study was made of 
the performance of a large, tailless, six-engined jet 
airplane. Three cases were studied, (a) operating 
conventionally on ramming inlet jet engines; (b) using 
the jet engines to suck off the boundary layer and apply- 
ing the data of reference 1 to all surface area to estimate 
the drag coefficient C,, [the induced drag was un- 
changed from case (a)]; (c) the same as case (b) ex- 
cept that the drag data of reference 1 were corrected 
to full scale by extrapolation along the 10°Cg = 0.2 
line of Fig. 24? (full-scalé cruising Reynolds Number 
varied from 60 X 10® to 10® depending on speed and 
altitude). 

An additional 15 per cent parasite drag increase over 
case (a) was added for cases (b) and (c) to account for 
inefficiencies of boundary-layer control on a complete 
airplane. Based on reference 1, N.A.C.A. tests, and 
some of our own test data, the total head into the en- 
gines was assumed to be static pressure minus 0.4q; 
0.25q corresponds to the peak negative pressure on the 
wing at cruising and the other 0.15g represents ducting 
losses. It appears that this value is somewhat conserv- 
ative. For case (a) 100 per cent ram was used. The 
engines used were late production engines. The study 
was made possible because the engine data were in a 
form suitable for calculation of thrust and fuel flow at 
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TABLE I 


Range Comparison 


Condition a b c 
: 100% ram Per reference 1 Per reference'1 extrapolated 
Ram = Static —0.4q¢ to full scale Reynolds Num- 
ber Ram = Static —0.4g 
Absolute range, miles 5,290 5,500 13,600 
Average cruise speed, m.p.h. 541 522 591 
Initial speed and altitude 541 at 35,000 ft. 589 at 25,000 ft. 600 at 15,000 ft. 
Final cruise altitude, ft. 45,000 35,000 35,000 
Power-plant operation All engines operatirig 6 engines at start, 4 halfway Engines Engines 
out, 3 at finish on Suc- on 
Normal r.p.m. never exceeded tion Ram 
Start 1.7 0.11 
Halfway 1.6 .14 
Finish 1.5 .15 


Normal r.p.m. never exceeded 


negative rams. The engine data made it possible to 
adjust the tailpipe area to obtain optimum perform- 
ance. Consequently, except for doubt about the total 
head at the engine, the power-plant calculations are 
straightforward, using manufacturers’ data, and should 
be considered as valid as for any ramming system. 


In Table I are the results of this study. The same 
fraction fuel was used in all cases; hence, cases (b) 
and (c) are optimistic by whatever added structural 
weight becomes necessary to accommodate slots. 


Calculation of performance for this type of operation 
is tremendously difficult because drag loses its inde- 
pendence of power. In case (a) the altitude was ad- 
justed at all times so that all the engines were operating 
at good efficiency. In cases (b) and (c) partial engine 
operation was permitted if better economy resulted. 
In case (b) no engines operated on ramming air. In 
case (c) operation was at a specified Cg of 0.0002 
based on wetted area. Based on airflow considerations 
only, the minimum number of engines required is then 
determined by the engine airflow at maximum cruising 
t.p.m. As thus determined, slightly more thrust was 
needed than could be given by the engines when they 
consumed this amount of air. The needed extra'thrust 
is indicated by the fractional number of ramming en- 
gines. Since so little extra thrust over that given by 
the suction engines is needed, it is apparent that the 
aerodynamic and power-plant airflow requirements are 
nearly equal. Of course, it should not be forgotten 
that the air requirements for case (c) basically come 
from porous plate calculations, although the airflow 
and drag relations check rather well with the Swiss 
test using discrete slots. 


Two means for obtaining the necessary extensive 
- laminar layers appear possible—porous surfaces or 
single or multiple slots. Examination of Dr. Gold- 
stein’s data shows that his values are little different 
from those assumed in the above study and so it can be 


considered as a rough guide to the possible overall ef- 
fects of boundary-layer control with any of the suction 
systems. From the study the following tentative con- 
clusions can be made: : 

(1) It may be possible to increase the range of jet 
aircraft two and one-half to three times the present-day 
values. 

(2) Cruising speeds may be somewhat higher than 
for the ramming-jet airplane because optimum opera- 
tion is at lower lift coefficients, and, if there are Mach 
Number limitations, optimum altitude becomes lower 
where the velocity of sound is higher. 

(3) Because the optimum altitude is lower, fears of 
combustion difficulties due to abnormally low: pressure 
in the burners seem unfounded. 

(4) There is a satisfactory correspondence between 
boundary-layer control and engine air requirements. 

(5) In the aircraft studied, there was sufficient 
space for the necessary ducting. 

(6) Because of the change from ramming to pump- 
ing operation, very large drag reductions must be made 
before a boundary-layer control system begins to be an 
improvement in overall economy. 

As previously mentioned, the above study appears 
applicable to the British work whether for jet or gas- 
turbine operation. The single-slot system appears to 
be the most desirable from a practical point of view. 
Consequently, since it shows possibilities of such very 
large range increases, it is hoped that the above infor- 
mation will serve as further encouragement to Dr. 
Goldstein in his work. 
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Warping and Shear Lag in Closed. 
Cylindrical Shells 


Me LEON BESKIN* 


= ABSTRACT 

gines 

on The state of stress in thin-walled closed cylindrical shells is de- 
am termined hereafter in the neighborhood of sections with concen- 
ll trated loads or of sections corresponding to discontinuities of sec- 
14 tion properties, such as fixed ends. A method of successive ap- 
15 proximations is developed, using the conventional beam formulas 
-eded as the first approximation. The second approximation, con- 


sidered as a correction to the first, is determined and examined 
in some detail in the cases of torsion and bending. This correc- 
tion may be regarded as a type of stress distribution having the 
-all ef- property of producing no section warping to a greater degree of 
accuracy than the conventional solution. Superposition of the 


oe two solutions gives the possibility of writing arbitrary warping 
e con- conditions at given sections. For instance, at a fully restrained 
: section (built-in cantilever), the nonwarping solution replaces the 
of jet conventional one, but, along the beam, the nonwarping solu- 
it-day tion damps out, following an exponential law, and is replaced by 
the conventional solution. In the case of a section with a con- 
, cefitrated load, it is possible to write that the mean warping due 
r than to the conventional distribution is the same on the two sides of 
opera- the loaded section, thus defining the portions of the loads taken 
Mach by each distribution on each side. A similar method applies to 
lower the case of a change of section properties in shear, if the rigidities 
in shear change in the same ratio at all points of the section. 
The validity of the method, in the case of structures stiffened by 
ars of frames, is limited by the fact that the frame spacing is often of 
essure the order of, or even larger than, the characteristic length along 
which the nonwarping solution damps out, while it is assumed 
that the frames are infinitely closely spaced. However, the 
tween fundamental aspects of the phenomenon remain unchanged, and 
nts. a more complete analysis could be made along the same lines. 
ficient 
PRINCIPAL NOTATIONS 
yumMp- do = total area of a section 
made ay, a2 = flange areas 
ie da = elementary area (= hds) 
an A, A = polar areas from twist centers 
Ao = total area inside section contour 
pears B = bimoment 
r gas- b,c = dimensions of a rectangular section 
irs to Cz, Cy, Czy = constants [Eq. (67)] 
view E, = moduli of elasticity and rigidity 
P h = thickness resisting axial loads 
very jm & = principal moments of inertia of a section 
infor- To = torsion bending constant 
» Dr. P i = coefficient of rigidity in twist 
Ji = secondary coefficient of rigidity in twist [Eq. 
(54)] 
K = constant [Eq. (47)] 
= constants [Eqs. (69)] 
riction L,1l;,I. = characteristic perturbation lengths [Eqs. (50) 
uction, and (70)] 
on in Received October 1, 1947. 
LCA: * Design Specialist. Now with Design and Research Division, 
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Consolidated Vultee Aircraft Corporation 


M = applied torque (about twist center) 

M = torque distributed according to conventional 
formulas 

AM = complementary torque (= M — M), perturba- 
tion 


M., M,, M, = components of the applied moment (bending 
moments and torque) 

M., M,, Mz = components of the applied moment distributed 
according to conventional formulas 


AM,, AM, = perturbations to the bending moments 

p = normal stress flow (= ho) 

P = constant [Eq. (32)] 

q = shear flow (= tr) 

Qs, Oy = static moments (= /xdaand /y da) 

r = distance between the tangent to a contour and 
the center of twist 

s = arc along the section contour 

t, ty, te = thicknesses resisting to shears 

u, = functions defining secondary stress distribu- 
tions [Eqs. (28) and (38)] 

uo = constant [Eqs. (39) and (45)] 

w, Aw = warping and secondary warping 

Awe, Aw = unit secondary warpings [Eqs. (38)] 

x,y = principal coordinates of a section 

Xo, Yo = principal coordinates of a specific point 

xT, YT = principal coordinates of the twist center 

xr’, yr’ = principal coordinates of the secondary twist 
center 

= applied shearing forces 

be f = shearing forces distributed according to con- 
ventional formulas 

AX, AY = complementary shearing forces (perturbations) 

a = polar areas from centroid 

Aa, Aa’ = additional constants to polar-areas 

Y = shear strain 

de = ds/t 

€, € = function [Eq. (24)], and additional constant 

. [Eq. (18)] 

6 = angle of twist 

Ar, Ae = constants [Eqs. (69)j 

B = constant [Eq. (78)] 

&, Ey = functions [Eqs. (24)] 

o = normal stress 

T = shear stress 


FUNDAMENTAL ASSUMPTIONS AND METHODS 


ie CONVENTIONAL BEAM THEORY is valid under 
restrictive assumptions, which, essentially, are the 
following: 

(1) The structure is a prism or a cylinder with con- 
stant section properties.. 

(2) The sections are free to deform in their planes. 

(3) The sections are free to warp out of their 
planes. 
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(4) The loads are applied only at the two ends of the 
structure or are uniformly distributed over the length 
(Surface and volume forces). 

(5) Applied loads are distributed in a section accord- 
ing to definite formulas, which result from the assump- 
tions of the beam theory. 

‘Under these conditions, the beam theory can be ap- 
plied to closed cylindrical shells—i.e., structures that 
have sufficient local stiffness to be stable under the ap- 
plied loads as opposed to membranes that have no de- 
termined shape and can coincide with all the surfaces 
which are applicable on one another. The stabilizing 
structures required in the case of membranes prevent 
condition (2) to be fulfilled and, to some degree, condi- 
tion (3). 

However, the limit case of membranes with closely 
spaced rigid frames does not differ essentially from the 
case of shells if the following simplifications are intro- 
duced: 

(1) The influence of lateral contraction (Poisson’s 
ratio) is disregarded, since this contraction is locally 
prevented by frames. 

(2) The material constituting the shell is considered 
as concentrated along its mean surface. 

(3) All plane sections are considered as indeformable 
in their planes. 

Two methods of analysis can be developed from these 
assumptions; both have been presented by the author 
in a condensed form :! 

(1) Representation of tke stresses by trigonometric 
series or trigonometric integrals of the coordinate 2, 
parallel to the generators of the beam. This method 
has also been developed by von Karman? in a somewhat 
different manner. The author wishes to acknowledge 
that a previous work of von Karman on shear lag 
theory’ is at the origin of his research work in this field.* 

(2) Representation of the stresses by successive ap- 
proximations, using the conventional beam theory as 
the first approximation. This second method is the 
object of this paper. 


SECTION I—THEORY 


(1) General Equations - 

When the cross sections of a cylindrical membrane 
are rigid in their planes but can warp out of plane, a de- 
formed state of the membrane is defined by two quan- 
tities: the angle of rotation @ between a section de- 
fined by its coordinate z along the generators of the 
cylinders, and the section z = 0, and the warping (dis- 
placement) w along the generators, function of z and 
of the coordinate s along a directrix (cross section) of 
the membrane. The center C of the infinitesimal rota- 
tion d@ between two sections may have a constant loca- 
tion or may be variable. It is defined by its coordin- 
ates (xo, Yo) with respect to the principal axes (x, y) in 


each section. 


The shearing strain being designated by vy, the follow- 
ing well-known relation exists between w, y, and 6, 


Ow/Os = y — 2(d0/dz)(dA/ds) (1) 


where 2dA/ds = r is the distance from the center of 
rotation C to the tangent to the directrix of the section 
(Fig. 1), while A is the corresponding polar area, taken 
with an arbitrary origin radius. The shearing strain 
is a function of the shear flow gq, 


y = (2) 


where ¢ is the effective thickness of the membrane re- 
sisting to shear and G is the modulus of rigidity. 

The warping is a function of the normal stress flow p 
along the generators, 


Ow/0z = o/E = p/Eh (3) 


where o is the normal stress along the generators and 
h is the effective thickness of the membrane resisting 
to normal stress. Thicknesses h and ¢ are often unequal 
because longitudinal members (stringers) are included 
only in hk and because buckling may affect differently 
the rigidities under axial loads and under shear. 
Eliminating y and w between Eqs. (1) to (3), it is 


found 
while # and q are related by the equation of equilibrium 
(Op/dz) + (0g/0s) = 0 (5) 


Relations (4) and (5) constitute a system of differen- 
tial equations in two unknowns, and their solution, 
with given conditions at the edges, gives the general ef- 
fect of the interaction between shear and warping. 

It should be noted that Eq. (4) is incomplete, since 
the term representing the displacement of the center 
of rotation is missing in its right-hand member; this 
term is 


—E—— = —E—— cos 8 (6) 


where s; is the projected arc defining the motion of C 
and 8 is the angle between s; and r (Fig. 1). This 
term is omitted because solutions considered hereafter 
correspond to a fixed axis of rotation or are obtained by 
superposition of such solutions. 

The conventional beam theory corresponds to the 
case when the shear flow is independent of z, so that 
Eq. (4) is simplified, and, after integration, it is found 


o=— = —2AE— = 2A— (7) 


with a proper choice of the origin radius of the areas A. 

Since g is independent of z, Eq. (5) shows that the 
nortnal stress o is a linear function of z, and Eq. (7) 
shows the normal stress in a section is proportional to 
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follow- <3 sumptions made here (see “Introduction’’). Correct 
8, methods, using trigonometric series or integrals may be 

warranted if more accurate physical assumptions are 

(1) made. This has been done by the writer for the prob- 

iter of lem of shear lag in I-beams,‘ developing a method first 

ection presented by von Karman.® 

taken In Section II, a comparison between the correct and 

strain the approximate method shows that in many cases, the 
approximate method gives results that are extremely 

@) close to those of the correct method. 

sae (2) General Relations for Closed Beams 
The conventional beam formulas, with conventional 
flow p notations (see also ‘‘Notations’’), are 
= —(M,/I,)x + (11) 
(3) 0 M, (12) 
s and the area A. This is the general case of torsion-bending 4 “> 
isting of open sections, where B is the bimoment about the with 
equal pole and J is the torsion-bending constant.* In the - 
luded case of conventional bending, the center of rotation is at GM = ©, = (19) 
ently infinity. A» is the area inside the contour of the beam section; 
The beam theory gives a first approximation of the the torque /, is taken about a stress center defined by 
it is solution of the system of differential Eqs. (4) and (5). its coordinates about the principal axes 
In a second approximation, the value of q calculated 
‘ from Eq. (5) is introduced into Eq. (4), giving, for p Xr = —(1/Iz) $20, da 
= (1/1,)$2Q,d (14) 
(4) and 6, corrections Ap and Aé related by r = (1/I,) $20, da 
OAc _ A r) E dq dA dA the polar areas a being measured from the centroid. 
rium h Gldz (8) It is known that if the static moments satisfy the con 
(5) ditions 
Eq. (5) then defines a correction Ag, £O0Ads/t) = 0, £0,(ds/t) = 0 (15) 
eren- 
tion, (OAp/dz) + (0Ag/ds) = 0 (9) Eqs. (14) define the coordinates of the conventional 
l ef- The process can be repeated with the system of equa- twist center, sometimes called shear center. With 
| tiene this choice of the reference point for torques M,, the 
hl expression of the angle of twist, given by the general 
this s\h Gt ds dz} (10) d0/dz = (1/2A) (q ds/Gt) (16) 
(2A%q/ds) + = 0 becomes 
(6 This method is examined here in the case of one u wiloaiiaciins jual (17) 
) term, and appears as a first order correction to the with ‘ai 
conventional beam theory. It is generally sufficient | on 0° _ ds/t 18 
for practical problems, at least with the physical as- 
fter 
lb ‘ 

y Expressions (14) of the coordinates of the twist center are generally derived by applying a dummy torque to a 
th section loaded by given transversal forces that can have an arbitrary line of support. The condition of zero rota- 
h ; tion of the dummy torque defines the twist center as the point through which must pass all the lines of support. 
. eo This method is, however, incomplete, since it does not show the necessity of the existence of a fixed end to define 
- a twist center; if there is no fixed end, any lateral displacement can be superimposed to the twist, which displaces 

the center of twist. Hereafter, the end restraint is explicitly introduced, and the resulting secondary stresses are 
(7) determined. Then the moments and shearing forces that appear in Eqs. (11) and (12) must be considered as the 
portions of the total moments and shearing forces that are transmitted by conventional beam distributions, so that 
A. it can be written 
the ‘ 
(7) M, = M, M, +> AMy, (19) 
R=X+AX, Y=V+ Ay 
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where the quantities with a bar represent the total applied loads and the quantities preceded by a A represent the 


portions of the loads transmitted as perturbations. 


The secondary normal stress Ac is defined by integration of Eq. (8), in which the shear flow q is given by Eq. 


(12), 


ger dz I, 


G (20) 


“dz. t dz? 


where the areas A are measured from the’ twist center (x7, yz) and are related to the areas a measured from the 


centroid by 


2A = 2a + — yxr + (21) 


The stress Ao must be a univalued function of the arc s,a condition that is satisfied by the first and second terms of 
the parentheses in Eq. (20), since Eqs. (16) are satisfied. Defining d?A@/dz? by the condition 


2A = 


(€)/2A) (dM,/dz) (22) 


the stress Ao becomes entirely a univalued function of s and can be written in the form 


( 
Ao = —| —~—-(-- 


GL dz 2Ao \e 


where 


e= Si, (ds/t) = Si, de, 
& = See Q, de 


A\_ 
4) dz Ij! dz I, (23) 


The unknown constants 5», s;, and s, are defined by the condition that the secondary stresses must have no re- 


sultant force, which yields three conditions 


F — (A/Ao)]h ds = 0 (25) 


ds = 0, 


These conditions can be written in the form 


= 0 


da = + $e da) = da Sede \ (28) 
aS, *Q,de= daf,'Q, de wf”Q,de = da de! 
where the total area resisting to normal stresses is do, 
a = § da, with da = hds (27) 


It can be verified that the normal stresses due to the torque term dM/,/dz have no resulting bending couple, 


which corresponds to the conditions 


= 0, fydu = 0, with du = [(€/@) — (A/Apo) (28) 
which can also be written ‘ 
Xp = —(2A0/Iz)F [(€/e0) — (a/Ao)]y (29) 
vr = (2Ao/I,)F [(€/e0) — (a/Ao)]x da 


These equations define the twist center, since they 
can be identified with Eqs. (14). For instance, it can 
be shown that 


F$Q,da = (Ao/e)feyda — Say da 


if the second condition (15) is satisfied. Since dQ, = 
y da, the relation above can be transformed, introduc- 
ing condition (15) and integrating by parts, 


which gives the identity 


Q,Ao 


Thus, the two definitions of the twist center are 
equivalent. The amount of numerical computations 
required by the two formulas is approximately the 
same. Egs. (29) define the twist center by the property 
that it is the point around which the torsion occurs 
under the action of a pure torque, without introducing 
bending stresses; in these equations, the origin for 
which ¢ and a are simultaneously equal to zero can be 
arbitrarily chosen on the contour of the shell. The 
definition (14) of the twist center corresponds to the 
property that, if a transversal load is applied at that 
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point, it produces no twist. By application of the 
principle of reciprocity (Betti- Maxwell), the two defini- 
tions are shown to be equivalent, which proves directly 
the identity of Eqs. (14) and (29). 

The use of the twist center as origin for torques M, 
due to shearing forces, makes twisting and bending def- 
ormations independent of one another, and they can 
be examined separately. 


(3) Secondary Stresses Under Pure Torque 


(a) Solution of the General Equations.—Eq. (23) is 
reduced to [see Eq. (28)] 
EdM e& du 


G dz 2Ayda (30) 


where the pure torque is represented by / instead of 
M,. 

On this equation appears the great difference be- 
tween axial stresses due to torque in open and closed 
membranes. For open membranes, Eq. (7) shows that 
these stresses are proportional to the integral of the 
torque (bimoment); for closed membranes, axial 
stresses are proportional to the derivative of the torque 
and, very often, they can be disregarded or considered 
as secondary stresses. 

Stresses o are equal to zero if there exists a pole, origin 


of areas, such that dA/de = rt = const. The constant 
is equal to Ay/€ so that 7 is defined by 
r = Ao/te (31) 


It is easy to verify for a given section if lines parallel to 
the tangents to the contour drawn at distances r, Eq. 
(31), from these tangents pass by a common point. If 
such point exists, there is no warping under torque. 
This is the case for a circular cylinder of constant thick- 
ness, for a regular polygonal prism of constant thick- 
ness, and for a rectangular prism in which the thickness 
of each side is proportional to its width. In these cases, 
the center of twist coincides with the geometric center, 
a property that is hidden when Eqs. (14) are used. In 
addition, the following property of the twist center 
results from Eqs. (28): In structures that do not warp 
under torque (du/da = 0), the twist center is independ- 
ent of the flange material (areas da), since Eqs. (28) 
are identically satisfied; in structures that warp under 
torque, the twist center location depends on the dis- 
tribution of the flange material. 

The quantity du/da [Eq. (28)] represents the warp- 
ing distribution; both « and A contain arbitrary con- 
stants bnt from the first Eq. (26) it results that du/da is 
entirely defined. Since the normal stress is equal to 
zero for at least four points of a convex section, one of 
these neutral points may be used as a common origin 
for A and «. 

The secondary shear flow is obtained by introducing 
expression (30) for o in Eq. (9), defining 


E 2 


be = Pu, 
G dz? 2A, 


In this expression, u contains an arbitrary constant 
of integration, which will be defined by the condition 
that the secondary warping due to the shear flow be 
minimum. To write this condition, the quadratic 
mean is used. The warping should be equal to zero, 
but this condition is not generally satisfied, except for 
nonwarping sections for which no correction is 
required. The secondary warping Aw is determined 
by an equation similar to Eq. (1), and it can be writ- 
ten 


Aw = (1/G) Aqde — 2A'(d A°0/dz) (33) 
with, 
2A’ = — yr'x — (34) 


The quantity dA°0/dz is determined by the condition 
that Aw is univalued, 


G 2A, u de (35) 


dz 2AG 


so that 
Aw = (P/G)[fude — (A'/Av) de] (36) 


The secondary center of twist (x7’, yr’), pole of the 
areas A’, and the constant ao’ are determined by the 
condition that the normal stresses due to secondary 
warping have zero resultant and zero resulting couple. 
These stresses being proportional to the derivative of 
the warping, the corresponding conditions can be 
written as: 


Aw da = 0, Awxda = 0, f Awyda = 0 


(37) 


Special solutions # and A@ are first defined for u and 
Aw, which satisfy the relations 


fide = 0, Ah = £Amda = 0 


(38) 
and the expressions for u and Aw become 
u=u+ 


Aw = = AD) + (39) 
where u% is the constant defined by the condition of 
minimum secondary warping, equivalent to a condi- 
tion of minimum strain energy due to secondary normal 
stresses, since these stresses in any given section are 
proportional to the warping in that section. Thus, the 
condition can be written ‘ 


5(f$ Ad’ da) = 0 (40) 


Substituting the value of Aw given by Eq. (39) in 
Egs. (37) and noting that the first of these equations 
is equivalent to the first Eq. (25), because of the last - 
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condition (38), it appears that 
= (41) 
The coordinates of the secondary twist center differ 


from those of the principal twist center, and, by com- 
parison with Eqs. (29), it is found 


The developed expression of Aw is thus 


AD = + evto(du/da) — 


(42) 
Yr 
(x/I,) Atiox da — (y/I,) da (43) 


and can be introduced into Eq. (40), thus defining a quadratic function of uo, the minimum of which corresponds to 


— (%/I,) F da — (y/Iz) y da](du/da)da 


(44 


(du/da)*da 


Eqs. (28) show that, after integration, the last two 
terms in the brackets are equal to zero; the first term, 
after introducing the definition (38) of A® and inte- 


grating by parts, yields the final expressiofi for 1 
uy = _ de (45) 
(du/da)* da 


The secondary torque AM, corresponding to the 
secondary shear flow Aq, is, from Wie (32), 


AM =2 Agqda = 


expression which can be written in the form 
AM = —(E/G)(d?M/dz?)K (46) 
with 
K= dA = i dA 
u + —— a u 


The second term can be rewritten as 


de du 
du\? 
of da 
because the contribution of Su de is zero [Eq. (38)] 
and because the integration by parts which leads to the 
last expression is also giving a zero term, since @ is a 


univalued function. Thus, replacing m by its value 
(45), the constant K has the final form 


(du/da)*da 


(6) Determination of Stress Distributions —From 
Eqs. (19) and (46), it results that the total torque M is 


M = M+ AM = M — 


When M = M—the conventional solution—the warp- 
ing is total. When = AM, there is no warping, to the 
accuracy of the present solution. If M7 is constant, the 
solution of Eq. (48)—the origin of the axis z being at 


(47) 


the fixed end of a cantilever—is 
M:= —e7*), AM = (49) 
where 
= VEK/G (50) 
At the fixed end, the normal stresses and shear flows 
are 
oo = (E/G)(M/L)(¢/2Ao) (du/da) (51) 
Meu 
240K £2udA 


Under shear flows Ag that produce no warping, there 
is a twist AQ, and the rate of twist, defined by Eq. (35), 
is 


(52) 


qo = Aqo = 


dA%/dz = AM/GJ, (53) 
with 
Ji = ff AqdA/Aqde = JK/euy (54) 


or 


= 14 {e[f(du/da)? da}?/ fu? de} (55) 


The shear flow at ipl is 


AMeu 
+559 
and the total rate of twist at any section is 
M,AM M AM 1 


Because of the effect of end fixity, the total angle of 
twist is decreased by 


ML $ (du/da)* da]? 
GJ fi’ de 
If M is variable (concentrated or distributed applied 


torques), Eq. (48) defines the distribution of M in its 
two components M and AM. On the two sides of a 
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section where a concentrated torque is applied, the 
distribution is represented on Fig. 2. 

(c) Strain Energy.—Methods of strain energy can be 
applied to problems in which secondary stresses due to 
torque are considered. The strain energy due to nor- 
mal stresses is, per unit length of beam, 


= (1/2)(EC/J*G*) (dM/dz)? (59) 


or 


with 


C = 4A? $ (du/da)* da (60) 


The strain energy due to shear is, per unit length of 
beam, 
W, = (1/2)$(q + Ag)*de/G 


or 


W, = (1/2GJ)[M? + (¢ouo/K)(M? — M*)] (61) 


(4) Secondary Stresses Under Pure Shearing Forces 
Eq. (23) is reduced to 


(Zé be 
G\dz 
and the secondary shear flow Ag is 


1 42 1 f aa) 63) 

G\ dz I, Pd; ( 
From Eqs. (26), it results that Ag is univalued. The 

origins of the integrals in Eq. (63) are defined by the 


condition that the resulting moment with respect to 
the twist center (x7, yr) is zero, 


SIASi,da=0, dASi,da = (64) 


The stresses Ao defined by Eq. (62) generally have a 
resulting moment to which corresponds a shearing 
force defined by its components 


d d 
ax = Acx da AY = — Aoy da 
ds oy (65) 


Using expression (62) for Ao and —- by parts, 
it is found 


(62) 


C= 


ay 
dz I, 


AX = = 
2 2 
G\I, dz I, dz (66) 
Ay 
with the notations 
= £0,’ de, = £0, de 
Cry = £ 2:0, de (67) 


Using expressions (19) for AX and AY, a system of 
simultaneous differential equations in two unknowns 
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3 


z 


Fic. 2. 


2 2 

E E L dz I, dz (68) 
G Cy@X C,@Y 


the solution of which gives the second approximation 
for the stress distribution in closed membranes near a 
section where the loads or the section properties are dis- 


continuous. 
The solutions of Eqs. (68) are of the type 


X = x ta + (69) 
VY = Y — + Kye) | 


in which the constants Z; and L» are defined by the 
equation 


[C, — L*(G/E)I,)[C, — L*(G/E)I,] — Cu? = 0 (70) 


The constants K can be chosen as follows: 
Ku = Li(G/E) — (C,/I,) 
Ky = L2?(G/E) — (C;/Iz) (71) 
Kyl, => Kyl, = Coy 


and constants \; and ), are defined by end conditions. 
For instance, if the warping is equal to zero (fully re- 
strained end), both X and Y must be equal to zero at 
that end. Then 


Xx = + 
Y= + 


Constants \; and ), being determined, the stress 
distribution is defined similarly as for torsion: The 
bending moments are divided into two components, 
one corresponding to the conventional solution and the 
other to the corrective terms: 

Conventional distribution, 


(72) 


2 
M; = M, ™ 
i=1 


(73) 
M, = mM, 
t=1 
Perturbation terms, 
2 
AM, = 
(74) 


2 
AM,. 


i=1 


: 
= X and Y is found eo. 
|| = 
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The normal stresses due to MW, and M, are defined a, t(h) a, 
by the conventional beam formulas; the normal 1 - oe 
stresses due to AM, and AM, are defined by Eq. (62), 
which can be written in the form c 
Ao = — =| — K,e7 
alt 
] 


Fic. 3. 

SECTION I]—EXAMPLES 
problem of a tube with a restrained end section under 
the action of a torque applied at the other end is inves- 

The section is defined by its dimensions 6 and c, the tigated (Fig. 3). It is assumed that the distance be- 
reduced thickness of the webs under shear, ¢, the re- tween ends is large as,compared with the characteristic 
duced thickness h of the webs and stiffeners under nor-_ length L, so that only decreasing exponentials should 
mal stresses, and the corner flanges areas a;. The be considered. ; 


(1) Symmetrical Rectangular Section 


The secondary normal stresses at the flanges at the fixed end, calculated by trigonometric integrals,' are defined 
by 


_ 4(6 f dee (76) 
T ra GhJo w(uo + coth Bw + coth yw) + 4 — (1/By) 
in which the shear stress r is the conventional stress for a tube under torque, 
= M/2bcet (77) 


and in which the other wobatiods are 
B = b/(b+ 0), y=c/b+c), w= 2a/h(b +c) (78) 


while w is a variable of integration. 

The solution given by Eq. (76) is considered as an ‘‘exact’’ solution, since it satisfies the physical assumptions. 
The approximate solution developed in this paper can be compared to it, in order to find the order of magnitude 
of the errors resulting from the consideration of a first order correction. 

The constants and functions that appear in the approximate solution have the following values, with the origin 
of s being in the middle of each face: 


e= s/t, € = 2(b — c)/t, Ay = be 


Expression In Face In Face c 
dA/ds c/4 b/4 
s(b — c) s(b — c) 
wh 4b(b + c) 4c(b + c) 
+ c) Sh (b + 2c)h _ ae 4 (c + 2b)h ab 
b 12 b+c 12 b+e 


(du/da)*da = [(b — c)/48t(b + + h(b + 0)] 
Side = — ~ *) + (b2 + c2 + 


32t\b + ¢ b+e 180 
b+e 


The maximum secondary normal stress at the flanges in which the shear stress 7 is given by Eq. (77). 
of the fixed end is defined by Table 1 gives some values of the secondary stresses, 
calculated by the ‘‘exact’”’ Eq. (76) and the approximate 

o/t =VE/KG[(b — c)/4] (80) Eq. (80), in the case where h = t. 
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TABLE 1 |! 
Secondary Stresses in a Rectangular Tube Under Torque, Re- ao 
strained at One End ve H 
+ ah b/c (exact) (approx.) L/(b + ©) 
1 2 0.49 0.475 0.46 ajet, 
1 4 0.83 0.835 0.47 . | > Cc x 
Comparing the exact and the approximate values, it a, . 3 - 1 
appears that the error is small when the flange area is { a, 
relatively important and becomes appreciable for small = ib >| 


flange areas. The values with zero flange area ob- 
tained by the approximate method are meaningless be- 
cause the theoretical value is infinite. 

When hk = O0—i.e., when the webs do not transmit 
normal stresses—the approximate and the exact 
methods give the same result, 


o/t = (b — c)VEt/2G(b + ca (81) 
and the length Z is given by the relation 


b+c Watt 48 (b + c)\1 + 3u/\60 


(82) 


in which yu is the ratio of the areas concentrated at the 
flanges, to the distributed areas [see Eq. (78)]. If 
h = 0, L/b + = If w = 1 and 
h = t, L/(b + ©) does not vary practically with b/c and 
is equal to while it reaches for small values of 
For a nearly circular tube, assuming that 


du/da = k sin (2xs/C;) (83) 


where & is a small quantity and C, is the half circum- 
ference, 


L/C, = (1/2%)VEh/Gt (84) 


andifh = t, L = C,/4. 

These values show that the length L is a fairly con- 
stant fraction of the perimeter of the tube and that the 
distance between frames may be smaller than L, so 
that the use of this method, which assumes a continuous 
distribution of frames, may lead to incorrect results. 

Remarks on Some Previous Solutions.—Many authors 
have previously analyzed the special problem of the rec- 
tangular section, examined here as an application of 
the general theory. The case of the tube without 
flanges and stiffeners has generally been, investigated, 
and many authors have found that the maximum nor- 
mal stress o is given by the simple relation 


o/t = V3E/G = 2.8 (85) 


In this case, Eq. (78) gives a good approximation of 
the correct value, while it is apparent that Eq. (85) is 
basically incorrect, since it indicates the presence of 
secondary normal stresses in a square tube prevented 
from warping and subjected to a torque, which is ob- 


Fic. 4. 


viously incorrect because a square tube does not warp 
under torque. 

This error is generally due to the fact that the second- 
ary warping is defined* by the condition g Ards = 0, 
instead of the condition of minimum warping used here. 
The error resides in the unjustified extension, to second- 
ary shear, of the relation £7 ds = 0, valid for the princi- 
pal shear. Asa result, the constant K, defined by Eq. 
(47), becomes 


K’ = (du/da)*da (86) 


This constant K’ tends toward zero for a square tube 
and, if used instead of K in Eq. (80), yields the incor- 
rect result represented by Eq. (85). 

In addition, it can be shown that the erroneous basic 
condition is equivalent to the assumption that the 
secondary shears produce no general shear deforma- 
tion. Since the secondary shears replace the primary 
shears at a fully restrained end, this assumption means 
that there is no shear deformation at that end, leading 
to too high values for the normal stresses. 


(2) Unsymmetrical Rectangular Section 


(a) The Webs Have Different Thicknesses and the 
Flanges Are Unequal.—The web does not transmit axial 
loads (h = 0) (Fig. 4). 

The centroid is located at an abscissa x, from the 
center of figure, 


The twist center is at an abscissa x7 from the center 
of figure, 


(88) 


(b) The Corner Flanges Are Replaced by Flanges on 
the y Axis.—Then x, = 0 and 


(c) Generalization of the Previous Formulas.—An 
approximately rectangular section such as used in two- 


. 
— 
‘ ; 
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spar wing structures can be considered as a combination 
of conditions (a) and (b): condition (a) represents spar 
flange areas and the skin and stringers near the flanges; 
condition (b) represents the skin and stringers toward 
the center of the section. 

Introducing the ‘‘web centroid,” which is the cen- 
troid of the spar web gages, 


X» = (b/2)[(2 — + (90) 
this point coincides with the twist center if 
b b ct 


The first relation is satisfied because the definition of 
the web centroid has been chosen for that purpose; 
the second relation is satisfied if the constant 9g, - 


= + te) (92) 


is equal to unity. In this case, there is no warping, 
and the twist cénter is independent of the flange ma- 
terial. In the general case, the location of the twist 
center is a function of the flange material, and it can 
be shown that in both problems (a) and (b) the follow- 
ing formula holds: 


Xp = Xp — = [2/(¢ + 1)](x» — 


In general, for a conventional wing section, ¢ is 
greater than unity, so that the twist center is between 
the web centroid and the flange centroid, and it is 
generally closer to the web centroid (gy < 3). The use 
of the above formulas can be generalized, considering 
them as approximate expressions for nonrectangular 
sections. In this case, the quantity designated by x, 
is better represented by the abscissa of the centroid of 
the flange moments of inertia, rather than by the cen- 
troid of the flanges. 

The centroid of the flange moments of inertia is 


= Sxy’ da/ $y’ da (94) 


and this permits to write an approximate formula for 
the twist center: 


Xr + [2/(¢ + 1)) Xe’) 


(93) 


(95) 


GENERAL CONCLUSIONS 


For a structure without end restraint which warps 
under torque, there is no definite twist center under 
constant torque, since the direction of a warped end 
section is not defined. If a section is prevented from 
warping, it is possible to define a twist center [Eqs. 
(14) or (29)]. Because these solutions are approxi- 


mate, a new twist center i¢ found at each successive 
approximation, except for non-warping structures, 
such as a circular cylinder of constant thickness. In 
this case, s.mple geometric considerations define the 
twist center. 
The effect of an end restraint or of a change of load- 
ing or of section properties is to produce a redistribu- 


1948 


tion of shears and of normal stresses, and the deviation 
from the conventional distribution follows approxi- 
mately an exponential law, rapidly decreasing on 
both sides of the section at which the perturbation is 
caused. 
There exists a characteristic length, designated by L, 
along which the main perturbation damps out; since 
e—! = 0.37, 63 per cent of the perturbation is elimi- 
nated along this length L. It has been shown that this 
length is small as compared with the perimeter of the 
section and is frequently smaller than the frame spac- 
ing. In these cases, the method developed in this paper 
does not apply, and an analysis that introduces the 
frame spacing in an explicit manner may be required. 


SUMMARY OF FORMULAS 


fren of the Twist Center About the Principal 
es 


There are two equivalent sets of formulas: 


with the conditions of definition of Q, and Q,, 
£Q:(ds/t)=0, £0,(ds/t) = 


ao 


with an arbitrary common origin for the following in- 
tegrals, 
e= f,ds/t, = = {rds 


B—Secondary Stresses in Torsion 

The applied torque J is divided into a warping com- 
ponent M corresponding to the conventional distribu- 
tion of shear flows q, 


(b) 


and a nonwarping component AM, corresponding to the 
distribution of shear flows Ag and normal stresses a, 


uAM EdM eo du 


$2udA’ G dz 2Ayda 
in which the variable u is defined as follows, 


with, the areas A being taken from the twist center, 
d _ du A % 
u_eé f 


I(A/Ao)da’= $ (€/e)da 


and 


= de)/[eo$ (du/da)*da} 
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The relation between M and AM is 
AM = —(E/G)K(d?M/dz") 


with 


K = of + | 


The differential equation defining M is 
M = M -— (E/G)K(d?M/dz*) 
At a fixed end, M = 0Oand AM = M;; then, if M is 
constant, 


M = L=V(E/G)K 


C—Secondary Stresses in Bending 


The applied loads X and Y are divided into warping 
components X and Y, and nonwarping components AX 
and AY. The second components give normal stresses 
Ao and shear flows Ag which, in principal coordinates, 
are 


with the following definitions of £,, &, 


while the origins of the integrals in Ag are defined by 
fi,da=0, dA ft,da=0 


The warping components X and Y correspond to con- 
ventional stress distributions, and the relations be- 
tween the two sets of components are 


ET C,d2X 
dz? I, dz? 
4 
dz? I, dz*: 
with 
C, = FQ, de, C, FQ,’ de, de 


For the determination of X, Y and AX, AY, see 
Eqs. (68) to (75). 
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Letter to the Editor 
Dear Sir: . course, for large M, t would have to be very small and Mr. Lai- 


We must protest against the view of ‘Supersonic Axial Flow 
About a Bedy of Revolution” presented by E, V. Laitone in the 
November issue of the JOURNAL OF THE AERONAUTICAL SCIENCES. 
It is true that, if in this problem the linearized solution be sub- 
stituted in the exact equation of motion, the error is of the order 
#2, where ¢ is the thickness ratio of the body. However, in refer- 
ence 1 it was shown that this does not affect the validity of the 
terms in ¢? in the linearized solution. This is confirmed by the 
complete discussion in reference 2, where it is shown that, using 
the exact equations throughout, the pressure coefficient on any 
body of revolution in an axial supersonic stream can be expressed 
as a sum of (specified) terms of orders 2°, ¢? log ¢, ¢4, t* log t, and 
t log? t, if terms of order #* log* ¢ be neglected, and that the terms 
of orders of ¢? and ?? log ¢ are the linearized solution as given in 
reference 1. Thus, the linearized solution is correct if terms of 


order ¢* log? ¢ be neglected. 
We must stress that, for every Mach Number M, the linearized 
solution is a good approximation for sufficiently small ¢, Of 


tone’s “‘hypersonic’’ results are then valuable. We see no reason, 
however, why doubt should be cast on the linearized solution, 
independently of t, for M > v2, a number apparently chosen 
arbitrarily. 

The objections we raise are of practical importance, since the 
use of a formula not based on sound theory, even if chosen to fit 
known facts in one special case (here, flow past a cone), may in 
others give widely wrong results. 
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Low-Speed Flutter and Its Physical — 
Interpretation 


M. A. BIOT* ann LEE ARNOLD 


Brown University and Aeronautical Consultant 


ABSTRACT 


The phenomenon of low-speed flutter of an oscillating airfoil 
is investigated. The problem resolves itself into the study.of two 
phenomena: that of zero air-speed classic flutter and that of an 
oscillating airfoil with no vortex shed. It is seen that 
both cases involve oscillation of the airfoil about a nodal line 
located three-quarters of the chord aft of the leading edge. 
The study is based on the fundamental theory developed by 
Theodorsen. 

The relationship between the parameters required for zero- 
speed flutter for the general case is developed, and the physical 
interpretation of this phenomenon is discussed. 

For the latter case, the solution of the equations of motion, 
based on noncirculatory flow, and the equations resulting from 
the Kutta condition result in a relationship between the inertia 
and elastic parameters for which the airfoil will flutter (i.e., main- 
tain oscillations at constant amplitude) for all air speeds. The 
physical interpretation and the mechanism of this phenomenon 
are discussed. 

Practical implications, such as the importance of nodal-line 
location and the effect of the physical parameters of the airfoil 
on nodal line location, are discussed. 

The possibility of a new approach to the flutter problem is sug- 
gested—that of employing the ground vibration modes of the 
airplane to study the existence of a low flutter speed and as a de- 
sign guide to raise the flutter speed. 

Conclusions are also drawn regarding the influence of aspect 
ratio and the relative importance of theoretical aerodynamic 
coefficients against values of these quantities measured in the 
wind tunnel. 


NOTATION 

h = downward displacement of airfoil at elastic 
axis 

a = pitch of airfoil (nose upward) 

V = air speed 

w = flutter frequency (rad. per sec.) 

le = moment of inertia of airfoil about elastic axis 

M = mass of airfoil 

Se = static moment of airfoil about the elastic axis 

‘x = Iq + mpb4(!/; + a?) = moment of inertia of 
airfoil including apparent mass effect 

M = M + zpb? = mass of airfoil plus apparent 
mass 

Si = Sa — arpb® = static moment including appar- 
ent mass effect 

Wh = natural frequency in h degree of freedom 

Wa = natural frequency in a degree of freedom 

b semichord 

ab elastic axis location aft of mid-chord 


= wh/V = reduced frequency 
= air density 


Received December 11, 1947. 
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Raa, Tea, etc. = real and imaginary parts'of aerodynamic coef- 
ficients, not including apparent mass effects 
g = structural damping 


INTRODUCTION 


UCH HAS BEEN DONE to improve the methods of 
flutter analysis during the last few years, but no 
tools have been devised which can be applied with reli- 
ance during the early stages of design. Various “‘tricks 
of the trade’’ are resorted to (such as the set rule of 
moving the center of gravity forward), but the re- 
minder that there have been exceptions to these rules 
compels one to search further for more reliable guides. 
With regard to flutter, the only controls the designer 
has are those that pertain to mass and elasticity distri- 
bution. However, the two are interdependent and the 
influence of the change of one on the other and the 
joint effect on the flutter characteristics are too complex 
to allow for the determination of suitable design guides 
begotten from intuition alone. 

By examining the basic two-dimensional flutter 
theory of Theodorsen,! one becomes aware of a feature 
of the aircraft which incorporates both the inertia and 
elastic characteristics. This feature, the so-called 
nodal-line location, is shown to be the primary factor 
in the occurrence of low-speed flutter. In fact, it is 
shown that for an idealized airfoil the location of the 
nodal line at the three-quarter chord results in a zero 
flutter speed. 


LINE 


If an external vibratory load is applied to the wing of 
an airplane on the ground, at each of certain frequencies 
of excitation maximum amplitudes of wing response 
would be induced. It would be seen that at each of 
these frequencies every point of the system would move 
in phase with every other point and that the motion of 
any section, say B of Fig. 1, would be essentially the ro- 
tation of the section about a point that remains sta- 
tionary in space. The locus of these nodal points is 
the nodal line. 

Let us examine the behavior of one of these wing sec- 
tions (Fig. 2). We assume that the section is free to 
oscillate in only two degrees of freedom. The first 
resonant condition will involve oscillations about a 
nodal point a large distance forward of the airfoil; the 
second will involve oscillations about a nodal point on 
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Fic. 1. 


ZERO AIR-SPEED FLUTTER 


Let us consider the effect of an air stream flowing 
over the airfoil of Fig. 2. At any air speed below the 
criticai flutter speed, any small disturbance of the air- 
foil from its static equilibrium position will result in an 
oscillatory motion that will damp out with time. As 
the flutter speed is approached, the rate of decay di- 
minishes until at the critical flutter speed there is no de- 
cay and the airfoil maintains a constant amplitude os- 
cillation. At this critical speed there is no exchange 
of energy between the mechanical system and the 
stream. At any speed above that of flutter, the me- 
chanical system absorbs energy resulting in oscillations 
whose amplitude builds up with time. 

Consider now the airfoil of Fig. 2 in still air. If we 


Nodal 
Line 


Plan form of wing and nodal line. 


neglect the structural and air friction, any displace- 
ment from equilibrium will result in constant amplitude 
oscillations. However, this phenomenon will be a 


FIc. 


limiting case of flutter for only special distributions of 
2, Isolated section of the wing. mass and elasticity of the system. It will now be 
shown that the necessary and sufficient condition for 
zero air-speed flutter is that the mass and elasticity of 


or near the airfoil. It is the latter point that is im- the system be such as to result in a node at the three- 
portant in the occurrence of low-speed wing flutter. quarter chord aft of the leading edge. 


The equations of motion of an airfoil oscillating in a uniform stream are: 


— wa*)a — = + + (Ran + tLan)(h/b)] (1) 
—w*S.a — M(w? — w,")h = + + (Ren + (h/b)) (2) 


where w, and w, are the natural frequencies of the system in still air, the former for pure bending and the latter for 


pure torsion. 


I.., 5,, and M are the moment of inertia and static moment about the elastic axis, and mass, re- 


spectively, all including apparent mass contributions. In nondimensional form, Eqs. (1) and (2) become 


E apb* + (Raa + ile) + +|- pb? ta) | =0 (3) 
M 
E | +|- ( + (Ra + ita) | = 0 (4) 
The stability determinant becomes , 
Wa? &, 
Sa M 1 
where Iga = etc. 
Separating into real and imaginary parts, 
[— 75; (1-25) + [- + 
= 0 (6) 
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Tea’ | - Se + Ra | 


- pb? (1 + Re 


The factor 1/k of Eq. (7), put equal to zero combined 
with Eq. (6), is the condition for steady-state oscilla- 
tions at zero air speed but not the condition for flutter. 
It is the simultaneous vanishing of Eq. (6) and the sec- 
ond factor of Eq. (7) which gives the flutter solution. 
At zero air speed 


1/k = Raa = Ran = Rea = Ra = 0 
and 
Toa’ = ("/2—4)*, = —(1/2 + @) 
Tea’ = */,-—a, Ia’ =1 
Substituting, Eqs. (6) and (7) become, 


apb4 w? arpb 
= 0 (8) 
apb*® apb? w? 
and 
We? 
( (1/2 — a) 
5. 
S. 
= 0 (9) 


(3/2 — a) (1-5) 


Eq. (8) is the usual zero air-speed frequency equa- 
tion. The important equation is Eq. (9). Expanding 
it, we obtain 


w? w? 
3 
(2a — 1)bSa + =0 (10) 


which is the equation governing the oscillating motion 
of an airfoil elastically restrained at the elastic axis but 
constrained by a node at the three-quarter chord point 
as shown in Fig. 3. 

The total restoring moment for the system of Fig. 3 
about the three-quarter chord is 


Mr = [b?(!/2 — + (11) 


The moment of the inertia forces about the three- 
quarter chord is 


M, = (Ia + M(*/s — a)*b? + 2b(a — 1/2)S.]& (12) 


The aerodynamic moment about the three-quarter 
chord point when 1/k = O is 


Ma = —(°/s)mpb* (13) 
We have the dynamical equation 
Mr + M, 4 Ma 


For harmonic motion this equation becomes 


1 2 
(2a — 1)bSa + = 


which is identical with Eq. (10), the condition for zero- 
speed flutter. 

We see then that the condition for zero air-speed 
flutter is that the nodal line is at the three-quarter 
chord. 

The importance of the nodal-line location in flutter 
is then apparent. One should therefore investigate the 
variation of flutter speed with nodal-line location as one 
of the parameters. This has been done for various 
combinations of parameter values. When the struc- 
tural damping g is not zero, the flutter speed does not 
vanish for a node location at the three-quarter chord; 
however, it reaches a minimum in the vicinity of that 
location. The variation of flutter speed with nodal- 
line location is illustrated in Fig. 4 for two values of the 
structural damping. 


VORTEX-FREE FLUTTER 


There is a case of occurrence of zero air-speed flutter 


in a two-dimensional airfoil for which flutter occurs. 


also at all other speeds. This type of flutter is as- 


Fic. 3. Mechanical and geometric parameters. 
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Fic. 4. Flutter speed versus nodal line location (node location 
in chord lengths). 


sociated with the vanishing of the circulation about the 
airfoil. The aerodynamic force and moment experi- 
enced by an oscillating airfoil is made up of two parts— 
a noncirculatory and a circulatory. If the latter, 
which also is associated with the shedding of vorticity 
from the trailing edge of the airfoil, is assumed to 
vanish, the equations of motion become 

ale + Sah + [- + 
k\k wb 


aM + + Saw + =0 (15) 


For harmonic motion, in nondimensional form, 
Eqs. (14) and (15) become 


I, 1 


(17) 


Since there is no circulation about the airfoil, the 
Kutta condition becomes 


Vat ht — ala = 0 (18) 


In nondimensional form 


From Eggs. (16) and (17), 


+ (X + (1/k)] [Z + _ 
| 


and from Eqs. (17) and (19) 


(1/2 — a) — i(1/k) 1 
[Z — i(1/k)] (21) 


where 
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X = (Iq/mpb*) [1 — (wa?/w?)] 
Y = (M/xpb?) [1 — (w,?/w*)] 
Z= Sa/ 


Eq. (20) becomes 


Y[X +:(1/k?)] — [Z? + (1/k*)] = 0 = (22) 
while Eq. (21) becomes: 
(/2-a)Y-Z= of 
Eqs. (23) are identically satisfied for all values of 1 /k if 
Y=1 
wad 
Substituting Eq. (24) into the determinant, Eq. (20), 
X = (1/2 — a)? (25) 
as the condition for flutter. Since Eq. (25) is inde- 
pendent of 1/k, it is seen that this type of instability 


will persist at all air speeds—i.e., flutter will occur at 
any speed from 6 to ©. 


From Eq. (19) it is seen that, when 1/k = 0(V = 0), 
h/b + ('/2 -a)a = 0 (26) 


which is the condition that the nodal point be at the 
three-quarter chord. A further examination of the 
Kutta condition, Eq. (18), reveals that the motion of the 
airfoil during this vortexless flutter is one of sinusoidal 
gliding at the three-quarter chord point as portrayed in 
Fig. 5. 


We have then the following conditions for vortexless 


flutter: 

( twa? (: : 

2 
y= =1 (ii) 

apb? w? 
Se 1 
Z xpb? a (iii) 
Node at three-quarter chord (iv) 


However, conditions (iv) and (iii) imply Eqs. (i) and 
(ii); therefore, the necessary and sufficient conditions 
for vortex-free flutter are, from Eqs. (iii) and (iv), that 


Fic. 5. Motion of the wing in vortex-free flutter, 
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the node is at the three-quarter chord point and that 
the static moment about the elastic axis is a small value 
determined by 


Sa = (27) 


The frequency of.flutter w is independent of the veloc- - 


ity, and the frequency ratio w,/w, is near unity. There- 
fore, when, in addition to a three-quarter chord loca- 
tion of the nodal line, the frequency ratio is near unity, 
the flutter speed will be extremely sensitive to changes 
_ in the frequency ratio. 


THREE-DIMENSIONAL CONSIDERATIONS 


We must realize that a wing may be assumed to be 
composed of a continuous distribution of two-dimen- 
sional sections of the type we have been describing. 
Each of these sections would have its own flutter char- 
acteristics if it were separated from the rest of the wing 
and placed in a two-dimensional stream. For example, 
section A of Fig. 1, whose nodal line is at 40 per cent 
of the chord, would have a high flutter speed, but sec- 
tion B, whose nodal line is close to the three-quarter 
chord, would have a low flutterspeed. Ata given speed, 
section A would be dissipating energy to the stream, 
while section B would be absorbing energy from the 
stream. At the air speed at which the amount of en- 
ergy lost by such sections as A equals that absorbed by 
those such as B, flutter occurs. The obvious guide is 
to design the wing so that the nodal line is as far from 
the three-quarter chord as is practicable at all spanwise 
stations of the wing. 

During the early stages of design, the nodal line and 
the effect of various design changes on its location can 
be determined analytically. 

Another important consequence of the present point 
of view is the tendency of vortex-free flutter to be inde- 


pendent of aspect ratio. It is to be expected that cer- 
tain types of flutter occur at or near the vortex-free 
régime. In this case, correction of three-dimensional 
aerodynamic effects will have little influence on the 
flutter speed. 

It should also be noted that vortex-free flutter in- 
volves only perfect fluid theory and does not introduce 
the action of viscosity since no circulation is produced. 
This is probably a contributing factor to the remark- 
able experimental verification of the theory which has 
been frequently observed. It may well be that the 
substitution of experimental aerodynamic coefficients 
measured on oscillating models can lead to less reliable 
predictions of flutter in certain cases. 


CONCLUSIONS 


(a) It is recommended that a study of locations of 
nodal lines be employed to ensure against the occurrence 
of flutter. 

(b) The designer should attempt to locate the nodal 
line as far from the three-quarter chord as is practic- 
able. 

(c) A limiting case of flutter without vorticity may 
occur. This may account for the small effect that as- 
pect ratio corrections have had on flutter in many cases 
studied. 

(d) The importance of the theoretical aerodynamic 
coefficients should be emphasized, since it is entirely 
possible that in certain types of flutter they may lead 
to more reliable predictions of the flutter speed than 
the experimental coefficients. 
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Letter to the Editor 


Dear Sir: 


I wish to make two corrections in the article written by myself 
and Mr. C. D. D. Pancu in the November, 1947, issue of the 
JOURNAL OF THE AERONAUTICAL SCIENCES, entitled ‘‘Use of 
Vectors in Vibration Measurement and Analysis.” 

First (through no fault of the publisher), the heading is mis- 
leading, since it gives the impression that I am at Consolidated 


Vultee Aircraft Corporation. 


I am not connected with Con- 


solidated Vultee and have not been connected with any other 
aircraft concern since I left the Curtiss-Wright Research Labo- 


ratory in 1944. 


Second, the titles of Figs. 3 and 6 should be interchanged. 
Fig. 3 should read ‘“‘Lissajous ellipse,’ and Fig. 6, ‘‘Force equi- 


librium.”’ 


CHARLES KENNEDY 
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The Analysis of a Circular Ring with 
Propped Floor Beam 


J. S. TAYLOR* ann S. S. GILL* 
4. V. Roe and Company Lid. 


ABSTRACT 


The paper sets out a method of analyzing a circular fuselage 
former with a floor crossbeam supported by four vertical struts. 
The structure has eight redundancies, but from symmetry of 
loading and geometry there are only five unknowns to be evalu- 
ated by strain-energy methods. 

The general equation for the bending moment and the partial 
differential coefficients with respect to the unknowns are stated 
in tabular form to assist computation. 


INTRODUCTION 


¥ thos TREND OF DESIGN of cabin aircraft has more or 
less standardized a type of fuselage former which 
consists of a ring and transverse floor beam. The load- 
ing of the fuselage consists in part of the loads applied 
to each floor beam by the internal equipment, etc. 
These loads are carried to the fuselage skin through the 
floor beam and ring and are counterbalanced by shear 
loads from the skin to the ring. 

In the case of noncircular rings, a graphical analysis 
is often necessary or convenient, but the advent of 
pressurized fuselages has required the use of circular 
rings of which there are many analytical solutions.'~* 

An example of an economical design of a former for 
an aircraft with circular fuselage is shown in Fig. 1. 
In this particular case, the floor beam is propped by 
four supports. It is assumed that these four vertical 
supports are pin-jointed at their ends and that the 
floor crossbeam is pin-jointed to the ring. However, 
note that the ring is considered continuous for bending 
at this point. This is consistent with practical 
design. The main ring would be a continuous ring, 
and the crossbeam and vertical struts would be riveted 
to it. 

This structure has eight redundancies, and, at first 
sight, its analysis by strain-energy methods along the 
lines of existing analyses**® presents formidable arith- 
metical difficulties because*of the work entailed. This 
paper presents a method of tabular layout of equations 
and computations which simplifies the otherwise tedious 
work. Combined with the fact that symmetry of load- 
ing and geometry reduces the unknowns to five, the 
solution by this method becomes tractable. 

The loading will be assumed to be a uniformly dis- 
tributed load over the floor beam, but the use of a 
similar analysis for other symmetric and antisymmetric 

Received August 27, 1947, 
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loadings would make it possible to solve a former of 
this type for any loading.® 

In this analysis the strain energy of the structure 
due to end loads, shear forces, and curvature of the 
ring is neglected. This is customary practice and is 
justified for large rings of small cross section.* 

The ring is cut at the top centerline, and the re- 
storing end load, X;, and end moment, M,, are applied. 
The vertical struts have loads W; and W, in them, and 
the floor beam acts as a tie as well as a beam. The end 
load in the beam, Xo, is the fifth unknown to be evalu- 
ated. . 

The total down load on the structure is 2wc. This 
is resisted by shear loads from the aircraft skin. It will 
be assumed that the distribution of this shear load (act- 
ing tangentially to the ring) is sinusoidal, representing 
a reasonable approximation. In actual fact, its dis- 
tribution will depend on the relative stiffnesses of the 
ring and the aircraft skin and their attachments, but 
widely. different assumptions of this distribution do not 
appreciably affect the final result. The moment 
due to the offset of the aircraft skin from the 
centroidal axis of the former is ignored because this is 
small. 


NOTATION 


x = horizontal distance of point on the former from point 
A (top centerline) 
y = vertical distance of point on former from point A (top 
centerline) 
restoring couple at cut at A 


M; 1 = 

X, = restoring load at cutatA 

X2 = end load in crossbeam 

W, = load in inner floor support strut 

We: = load in outer floor support strut 

w = uniformly distributed floor load per in. run 

I, = moment of inertia of former cross section 

Iz = moment of inertia of floor beam cross section 

= 

R_ = radius of neutral axis of former 

@ = angle between point on former under consideration and 
vertical centerline 

M, = bending moment at any point of the former caused by 


shear load from aircraft skin applied to former be- 
tween A and the point under consideration. 

M = bending-moment at any point caused by all loads ap- 
plied to the former 

U = total strain energy 


Other symbols used are obvious from the sketches of the former 
or from their context. 
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TABLE 1 
Bending Moment Due to ~ oM OM OM 
Part M, Xi 2 1 W, w OM, OX, ow, ow, 
AB 1 y M, 1 y we 
BC 1 y y—d x—c c(¢c — x) M, 1 y y—d 
CD 1 y y—d x—c c(c — x) M, 1 y y—d x—-c 
DE | y y—d a-c¢ b-—c c(c — x) M, 1 y y—d a-c b—¢ 
FG a-c¢ b-c (c? — x*)/2 a-—-c b—¢ 
GH x—c b-c (c? — x?)/2 x—-c b—¢ 
HB x—c x—c (c? — x?)/2 
ae TABLE 2 
Part M, Xi Xo W, w M, Xi Xo W; w 
4 AB 1 y y y? ae My 
BC 1 y y—-d x-—-c x-—c cc—x) M, y yy—d) yx—c) — x) My 
CD 1 youd? e-c b—c M, y WHy—d) wWb—c) My 
| HB 
TABLE 3 
(M/I)(0M/0X2) 
4 Part M, Wi 2 w s 
‘ BC y-d yy — d) (y — d)? (y — d)(x — c) (y — d)(x — ¢) c(c — x)(y — d) M,(y — d) 
ae CD y—d yy — d) (y — dj? (y — d)(x — c) (y — d)(b — c) c(c — x)(y — d) M,(y — d) 
— d) (y — d)? (y — d)(a — c) (y — d)(b — c(c — x)(y — d) — d) 
HB 
(M/1)(0M/oW,) 
Part M, Xi Xe 2 w 
; BC x—-c¢ y(x — c) (y — d)(x — c) (x — c)? (x — c)? —c(c — x)? M,(x — c) 
»| CD x—C y(x — c) (y — d)(x — c) (x — c)? (x — c)(b — c) —c(c — x)? M,(x — c) 
DE a-c y(a — c) (y — d)(a — c) (a — c)? (a — c)(b — c) x)(a — M,(a — c) 
FG k(a — c)? k(a — c)(b — c) ae Ke 
| HB k(x — c)? k(x — c)? 
TABLE 5 
Part. (M/1)(0M/oW2) 
M, Xi Xe W, 2 w 
My BC x-c¢ y(x — c) (y — d)(x — c) (x — c)? (x — c)? —c(c — x)? M,(x — c) 
i CD b-c¢ y(b — c) (y — d)(b — c) (b — c)(x — c) (b — c)? c(c — x)(b — c) M,(b — c) 
! DE y(b — c) (y — d)(b — c) (6 — c)(a — c) (b — c)? — — M,(b — c) 
GH k(b — c)(x — 0) k(b — c)? a0. =. 
HB k(x — c)? k(x — c)? 
DERIVATION OF EQUATIONS M,+ Xyt+ + 
The equations for the bending moment at any 


point in the ring and floor beam are given in tabular 
form in Table 1. 

: The bending moment (B.M.) at any point on the 
” ring between C and D, for example, equals 


and the bending moment at any point of the floor beam 
between G and H, for example, equals 


(x — c)W, + (6 — + [w(c? — x*)/2] 


For 1 
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CIRCULAR RING WITH 


The coefficient of each applied load or redundancy is 
tabulated in the appropriate column for each part of the 
ring or beam. The tables of 0M/0M,, 0M/0X1,0M + 
0X2, OM O0M/OW, for each part of the structure 
in the same table are self-explanatory. 

The next step is to tabulate M(OM/0M)), etc., for 
each part of the former. 

The coefficients of the unknowns and applied loads 
for each part of the ring are multiplied by the appro- 
priate differential coefficient (see Table 1), giving the 
Tables 2 to 5. 

Each ‘‘block” of coefficients gives one of the five 
equations 


ET0U/0M, = 0, EI,0U/o0X, = 0, 
EI,0U/0X2 = 0, EI,0U/0W, = 0, 
EI,0U/0W2 = 0 


As an example from Table 2 


EI,0U/OM, = M, ds + X1 yds + X2 X 
(y — d)ds + Wil fp? (x — dds + (a — ods] + 
(x — olds + (b — c)ds] + we X 

(c —x)ds + Mids 


Each coefficient is multiplied by ds and is integrated 
over the appropriate range. Note that integrals 
for the beam are multiplied by k, the ratio of the mo- 
ments of inertia of the ring and the beam. 

Examination will show that all these coefficients 
are multiples, or the addifion of multiples, of the fol- 
lowing standard “‘geometric’”’ and ‘‘loading”’ integrals, 
which must be evaluated for each part of the ring or 
beam: 


For the ring: 


S ds, xds, yds, y*ds, xyds, x*ds, 
S Mads, M.yds, M,xds 


For the beam: 
S ds, xds, x*ds, J 


Integrals for the ring are required over the ranges 
AB, BC, CD, DE. Integrals for the beam are re- 
quired over the ranges FG, GH, HB. All other ranges 
can be found from these basic ones by addition. 

Note that in practice it is simpler and more con- 
venient to write the Tables 1 to 5 side by side on a long 
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SHEAR LOAD 
FROM SKIN. 


Winch 
X, AW, I, 
W,- 
2 
Cc 
E|~ oO 
NUMERICAL DATA: 
R=SO 
b=25 
d=80 Tel, 
Fic. 1. 
sheet of paper. The tables are split up for convenience 


in printing. 

The various integrals are evaluated for a ring whose 
dimensions, etc., are given in Fig. 1. First they are 
evaluated analytically for any range @, to ®. The 
general mathematical analysis is given below. On 
Tables 6 and 7, 6; and 4. and x; and x2 are taken at their 
appropriate values for the different ranges, and the 
trigonometric functions and integrals for each range 
are tabulated. 

Note that these tables are considerably abbreviated. 
In practice it is easier to tabulate a large number of 
small items, as is done in the first eleven lines, thus 
keeping all the calculations in tabular form and avoiding 
separate calculations. 


ANALYTICAL EVALUATION OF THE INTEGRALS 


For the Ring 


Se ds = — 6) 
So” yds = R*[(6. — 6;) — (sin 6. — sin 6,)] 
Je xds = —R*(cos 6. — cos 6;) 
Je." xyds = R*[1/2(cos? 6. — cos? 6:) — (cos 6 — cos 6;)] 
= R°[3/2(0. 2(sin 6 — sin + 1/4(sin 26, — sin 26;) 
Se” x*ds = R*[1/2(6. — 0;) — 1/4(sin 26. — sin 26,)] 


OM 
oW, A 
b-—¢ 
b-—¢ 
\ 
: 
My 
My 
eee . . + 
— d) 
— d) 
— d) 
—c) 
c) 
— c) 
c) 
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TABLE 6 
Range —— 
d 3 , Item AB BC CD DE 
2.2143 2.6180 2.9403 3.1416 
6; 0 2.2143 2.6180 2.9403 
02 — & 2.2143 0.4037 0.3223 0.2013 
6,2 4.9031 6. 8539 8.6454 9.8697 
: 0,2 0 4.9031 6.8539 8.6454 
022 — 0,7 4.9031 1.9508 1.7915 1.2243 
4.4286 5.2360 5.8806 6. 2832 
26; 0 4.4286 5.2360 5.8806 
sin 262 — 0.9600 —0.8660 —0.3918 0 
sin 20, 0 —0.9600 —0.8660 —0.3918 
: sin 26. — sin 26, —0.9600 0.0940 0.4742 0.3918 
sin 6. — sin 0.8000 —0.3000 —0.3000 —0.2000 
cos 6. — cos — 1.6000 —0.2660 —0.1138 —0.0202 
Zz. cos? @ — cos? 6; —0.6400 0.3899 0.2101 0.0400 
. (8. — 0) — (sin 6 — sin 6) 1.4143 0.7037 0.6223 0.4013 
, 1/,(cos? 6. — cos? 6) —(cos 6 — cos 4) 1.2800 0.4609 0.2188 0.0402 
3/2(0. — 0) — 2(sin — sin + 
pate 1/,(sin 26. — sin 20,) 1.4815 1.2291 1.2020 0.79995 
— 0)*—1/4(sin 26. — sin 20;) 1.34715 0.17835 0.04265 0.00270 
02 COS 02 — cos — 1.3286 —0.9368 —0.6137 —0.2607 
sin 20. — sin 20, —2.1257 —0.1415 1.1152 1.152 
ae ; 62 COS 202 — & cos 26; —0.6202 1.9211 1.3961 0.4365 
Total: A = — 0) 2.2143 0.4037 0.3223 0.2013 
; ‘ —3/,(sin 6. — sin 0) —1.2000 0.4500 0.4500 0.3000 
cos 02 — 0; cos — 0.6643 —0.4693 — 0.30685 —0. 13035 
: ’ A 0.3500 0.3844 0.46545 0.37095 
Total: B = —2(6.2 — 6,7) —9.8062 —3.9016 —3.5830 —2.4486 
—16(cos 6. — cos 6;) 25.6000 4.2560 1.8208 0.3232 
J +2(6. sin 20. — sin 26,) —4,2514 — 0.2830 2.2304 2.3040 
Dy +5(cos 26, — cos 26;) —6.4005 3.9005 2.1000 0.4000 
4 B 5.1419 3.9719 2.5682 0.5786 
Total: C + 24(@. — 53.1430 9.6889 7.7352 4.8312 
“t —40(sin 6. — sin 4) — 32.0000 12.0000 12.0000 8.0000 
8(62 cos 02 — cos — 10.6288 —7.5088 —4.9096 —2.0856 
: 5(sin 20. — sin 26,) —4.8000 0.4700 2.3710 1.9590 
+ cos 262 — cos 26;) 1.2404 —3.8422 —2.7922 — 0.8730 
: "34 6.9546 10.8079 14.4044 11.8320 
JS ds 110.715 20.185 16.115 10.065 
i T yds 3,535.75 1,759.25 1,555.25 1,003 . 25 
JS xds 4,000.0 665.0 284.5 50.5 
: JS xyds 160,000 .0 57,612.5 27,350.0 5,025.0 
yids 185, 187.5 153,637.5 150,250.0 99,993.75 
S x*ds 168,393.75 22,293.75 5,331.25 337.5 
f x A —22,282.0 —24,472.0 —29,631.0 —23,615.0 
f Meds = x B —1,022,950.0 —790,186.0 —510,928.0 —115,109.0 
f = x —1,383,576.0 —2,150,167.0 —2,865,669.0 —2,353,905.0 
A For the Beam 
A St ds = x2 — X1, xds = (x2? — x7) /2 
" INTEGRALS DUE TO SHEAR LOAD 
? 
Distribution of Shear Load 
; Assume the distribution varies as sin ¢ (see Fig. 2) equals f sin ¢. Then vertical component equals f sin’ ¢. 
i Therefore, total vertical force = fo" f sin? ¢ Rd ¢ = f(Rr)/2 = we 
Therefore, 
f = 2we/Re 
Qwe 
The bending moment at 0 due to the shear force at ¢ equals — Re sin ¢ {R[1 — cos(6 — ¢)]}Rd¢ 
Therefore, 
6 
2we 2wcR 6 sin 
= 2% sin cos (@ = 1— — cos 6 
Rr Jo 2 
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Hence, we have 


Se" = M,Rdo 
where M, = M,atangle@. Therefore, 


2wcR? 6 2wcR? 3 
Mds = (1 cos 0) dé = — 6) — —(sin 6 — sin + 
T a 2 2 


1 
5 (02 cos 6, — 6; cos a) | 


Similarly, 


62 

2weR® 

4 f M.xds = = [—2(6,2 — 6,2) + 5(cos 20 — cos 20;) + 2(62 sin 20 — 8; sin 20,) — 16(cos 6 — cos 61)] 
iT 

and 


02 | 
2wecR* + ‘0 
M,yds = [24(0. — 0) — 40(sin 6 — sin 6) + 5 (sin 26, — sin 26) cos 62 — 6; cos 
2(02 cos 20, — 6, cos 26;)] 


TABLE 7 

Item FG GH HB 

n 0 10 25 
Xe 10 25 40 
S ds = 10 1 15 
= — 50 262.5 487.5 

3 

— 333.33 4,875 16,125 
Sxids = — 2,500 95,156.25 542,344 


BENDING 
MOMENT 
DIAGRAM 


LINEAR SCALE: 


10 
B.M. SCALE: 
—ve 
SHEAR LOAD. 
Fic. 2. \ 
\ 
From these basic integrals we can evaluate numeri- ¥\ 
cally the coefficients of the partial differential coeffi- \ 


cient equations as tabulated in Tables 2 to 5. 

The layout is identical, each square representing 
the coefficient. For example, coefficient of W; in 
EI,(0U/0W,) = Ois: 


(x — + SC? (x — + 


Fie. 3. 


(a — +k (a — + 
k (x — oc)? ds +k (x — o)*ds 


This is the sum of multiples of /x°ds, {xds, and 
JS ds over ranges BC, CD, DE, FG, GH, HB. 


& 
DE 
3.1416 
2.9403 
0.2013 
9.8697 
8. 6454 
1.2243 
6. 2832 
5.8806 
0 
-0.3918 
0.3918 
-0. 2000 
-0. 0202 
0.0400 
0.4013 
0.0402 
). 79995 
). 00270 
-0. 2607 
1.152 
0.4365 
0.2013 
2.4486 | 
0.3232 ; 
2.3040 
0.4000 
0.5786 
4.8312 
8.0000 | 
2.0856 
1.9590 
0.8730 
1.8320 
10.065 | 
003.25 
50.5 4. 
1,025.0 
993.75 
337.5 
615.0 
109.0 
YS 
500 w.lb.ins. 
+ve 
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TABLE 8 
(M/I)(0M/0X;) 
» Part M, 2 w 
AB 3,535.75 185,187.5 — 1,383,576 
BC 1,759.25 153,637 .5 + 153,637 .5 +57,612.5 +57,612.5 +2,214,800 — 2,150,167 
— 140,740.0 —70,370.0 —70,370.0 — 2,304,500 
CD 1,555.75 150,250.0 +150,250.0 +27,350.0 — 23,336 .0 +2,489,200 — 2,865,669 
— 124,460.0 —62,230.0 — 1,094,000 
DE 1,003.25 99,993.75 +99,993.75 —30,097.5 — 15,048.0 +1,605,200 — 2,353,905 
— 80,260.0 — 201,000 
TOTAL 7,854.0 589,068.75 + 50,421.0 —77,735.0 — 51,142 +3,309,700 — 8,753,317 
TABLE 9 
_ (M/I)(0M/0W,) 
Part M, Xi X2 W, i w 
BC +665.0 +57,612.5 +64,592.0 +22,293.75 +22,293.75 — 1,291,840 — 790,186 
— 53,200.0 — 53,200.0 — 53,200.0 +2,128,000 +978,880 
— 807.4 —70,370.0 +57,612.5 +32,296.0 +32,296.0 — 891,750 
—70,370.0 
CD +254.5 +27,350.0 +5,331.25 — 4,267.5 — 1,031,360 — 510,928 
— 22,760. 
— 645.6 — 62,230.0 +27,350.0 — 22,760.0 +9,669.0 +910,400 +1,185,240 
— 62,230.0 +25,784.0 — 213,250 
DE —301.95 —30,097.5 —30,097.5 +9,058.5 +4,529.25 +60,600 +708,450 
+24,155.0 — 483,120 
FG +9,000k +4,500k — 240,000k 
+5,000k 
GH +4,875k —3,937k —47,578k 
—21,000k +9,000k +210,000k 
— 24,000k —480,000k 
+97,500k 
HB +16,125k * +16,125k —271,172k 
—39,000k —39,000k +390,000k 
+24,000k +24,000k —480,000k 
+322,500k 
TOTAL — 805.45 —77,735 — 13,379 + 18,803 .5 +11,320.5 — 812,320 + 1,571,456 
+18,000k +10,688k — 493,750k 


coefficient = 18,803.5 + 18,000 k 


(See Table 9.) 
The other coefficients are tabulated and evaluated in 


the same manner. 

Note that Maxwell’s reciprocal theorem is a useful 
check on this work—e.g., coefficient of W; in (M/I) X 
(OM/0W:) = coefficient of in (M/l)(0M/0W,) 
and so on. 

Table 8 shows the computations for the equation 
0U/0X,; = 0. Table 9 shows the computations for the 
equation 0U/OW, = 0. Similar tables are required for 
0U/0M,, 0U/0X2, OU/OW2, but these are omitted to 


save space. 
The five simultaneous equations are: 


157.08M, + 7,854X, + 609.05X, — 805.45W, — 
535.10W2 + 34,184w — 100,000w = 0 
7,854M, + 589,068.75X1 + 58,421X, — 77,735W, — 
51,142W.2 + 3,309,700w — 8,753,317w = 0 
609.05M, + 58,421X, + 9,472.25X_ — 13,379W, — 
8,305W2 + 573,800w — 1,152,301w = 0 
—805.45M, — 77,735X, — 13,379X2 + (18,803.5 + 
18,000k) + (11,320 + 10,688k) W2 — (812,320 + 
493,750k)w + 1,571,456w = 0 
—535.10M, — 51,142X, — 8,305X2 + (11,320 + 
10,688k) W; + (7,280 + 6,750k) We — (482,910 + 
299,610k)w + 987,384w = 0 


These five equations may be solved simultaneously 


for M,, X1, X2, Wi, W2 for any value of k. If k = 1, 
M, = —95.62w, X, = +9.44w, X, = +50.05w 
W, = +22.54w, We. = +10.52w 


The bending moment diagram for this case is shown 
in Fig. 3. Positive bending moment tends to decrease 
the curvature of the ring and, for the beam, tends to 
make the beam concave upward. 
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High-Speed Flow Through Cambered 
Rotating Grids 


FRANK L. WATTENDORF* 
Air Materiel Command, Wright Field 


SUMMARY 


In a previous paper,' an analysis was given for the ideal per- 
formance of rotating blades in curved lattice arrangement. The 
lattice type of blading represents an approach to higher stage 
pressures by utilizing blades that are more closely spaced and 
highly cambered than airfoil theory permits. The blades exert 

a predominant guiding action on the flow, and the direction of 
flow relative to the trailing edge of a blade is esseptially deter- 
mined by the angular setting of the trailing edge itself. 

For incompressible flow, ideal pressure coefficients based on 
blade tip speed can be computed as a function of the flow coef- 
ficient. Experimental results of several investigators show good 
comparison with the theory in the neighborhood of the design 
point. 

For compressor application, the stage temperature ratio or 
pressure ratio replaces the pressure coefficient in significance. 
The temperature ratio can also be expressed in terms of the veloc- 
ity diagram and Euler’s equation. 

The present paper shows various configurations of flow dia- 
gram which result in considerable increases in the ideal stage 
pressure ratio, and the following four cases are considered: 

(1) Subsonic flow throughout. 

(2) Subsonic axial flow, subsonic flow in the rotor, supersonic 
flow in the downstream stator. 

(3) Subsonic axial flow, supersonic flow in rotor and stator. 

(4) Supersonic flow throughout. 

In all cases considerations are limited to ideal adiabatic flow 
through curved lattice types of blading. 

It is the purpose of this paper to show the extremely large 
potential increases in ideal stage pressure ratio theoretically 
possible in axial compressors, with a view toward stimulating 
experimental investigation of the losses involved and of possible 
methods of decreasing these losses. Application would be espe- 
cially promising for cases where simplicity and compactness out- 
weigh efficiency in importance. 


INTRODUCTION 


COMPRESSORS are becoming increasingly im- 
portant as components of gas-turbine propulsion 
systems and other applications because of their small 
frontal area in relation to their output. Their chief dis- 
advantage, as compared for instance with the centrifu- 
gal compressor, lies in the many stages and resulting 
length of the unit required to obtain a given pressure 
ratio. A subsonic axial-flow compressor designed by 
airfoil theory seldom has a stage pressure ratio of more 
than 1.25; therefore, such a compressor would need 
many stages to equal a centrifugal compressor having 
a stage pressure ratio of 3 or 4 to 1. 
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It is the purpose of this paper to show how, disregard- 
ing for the moment the question of efficiency, the stage 
pressure ratio of an axial compressor cam be raised to 
compete with the centrifugal compressor. This particu- 
lar discussion is based on the utilization of closely 
spaced, highly cambered blades, and the four cases 
mentioned in the Summary are considered. 

Ideal adiabatic flow will be assumed in this study, and 
none of the cases considered should be assumed competi- 
tive with other methods from the point of view of 
efficiency. 


AIRFOIL THEORY 


Axial compressors are usually designed by the applica- 
tion of single-airfoil theory, whereby each blade element 
is treated as an airfoil and the pressure produced is ob- 
tained by suitable resolution and integration of the 
component forces acting on the blade element. A 
typical vector diagram showing the forces acting on a 
blade element is shown in Fig. 1, where W, and W, are 
the relative inlet and the relative outlet velocities of 
the rotating blades, Wy x. is the absolute outlet velocity 
of the upstream contravane or stator, and W ia is the 
absolute inlet velocity of the downstream contravane. 


Upstream Stator Gutflow. 


{Downstream Stator Inflow. 
Z ‘ 
> 
6 Bm Be 
Ur U, 
Rotation 
Fic. 1. Velocity diagram for a single-stage compressor (airfoil 


theory). 
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Fic. 2. Velocity diagram for a single-stage compressor (lattice 
theory). 


U, is the circumferential velocity of the rotating blade 
element, and W,, is the mean relative velocity through 
the rotating blades. Since the pressure rise is essentially 
proportional to the lift force acting on the airfoil, it 
follows that increased stage pressure should be obtained 
in two principal ways—namely, increase of the relative 
velocity, Wm, or increase of blade area. 

The relative velocity is usually limited by approach 
to the speed of sound in normal compressor design. 
Utilization of supersonic velocity in the airfoil type of 
compressor will not be discussed in this paper. 

Increase of blade area involves increase of the ratio 
of blade chord to spacing. When this ratio essentially 
exceeds unity, mutual interference between the blades 
seriously impedes the applicability of single-airfoil the- 
ory, and a different approach becomes necessary, as 
follows. 


CURVED LATTICE THEORY 


The curved lattice theory, as discussed in reference 1, 
represents an approach to high stage pressures for any 
given relative velocity by utilizing higher values of the 
ratio of chord to spacing than airfoil theory permits. 
This is done by using a system of cambered, closely 
spaced blades, the passages between which represent 
curved expanding channels and are analyzed as such. 
The blades exert a predominant guiding action on the 
flow, and the direction of flow relative to the trailing 
edge of a blade is essentially determined by the angular 
setting of the trailing edge itself, as discussed, for in- 
stance, by Constant.? 

A typical vector flow diagram for a lattice system is 


- shown in Fig. 2, where the notations for the velocities are 


the same as that of Fig. 1. Here, Up is the change in the 
circumferential component of the velocity relative to 


the rotating blades. In reference 1 the pressure coef- 
ficient was expressed in terms of the flow coefficient and 
was integrated over the effective disc area. In this way, 
ideal characteristic curves could be calculated for a 
series of curved lattice systems. Experimental results 
showed good comparison with the theory in the neigh- 
borhood of the design point. 


STAGE TEMPERATURE RATIO 


For compressor application, the stage pressure ratio 
or temperature ratio replaces the incompressible pres- 
sure coefficient in significance. For compressible flow, 
the adiabatic temperature head developed by a com- 
pressor stage can be expressed in terms of Euler’s equa- 
tion, which ties directly into the velocity diagram, 
shown in Fig. 2. The force acting against the direction 
of motion f the rotating blades per unit mass flow per 
second is 1- Uy. The work done per unit mass flow per 
second by the rotating blades is then 1-Uy-U,. If the 
absolute inlet velocity to the rotating blades is equal to 
the absolute outlet velocity from the stationary blades, 
there is no change in the kinetic energy of the gas stream 
in passing through the stage. According to the energy 
equation for gas flow, this work done, UrU,, must ap- 
pear as increase in enthalpy of the gas. Thus, 


JgC, AT = UrU, 


where J is the mechanical equivalent of heat, g is the 
gravitational constant, C, is the specific heat of the gas 
at constant pressure, and AT is the stage temperature 
rise. The stage temperature ratio then becomes 


Tit aT _ |, Url, 
Ti JgC,T; 


Since isentropic flow is assumed, the pressure ratio 
p2/p; across the stage must be given by 


= [1 + (UrU,/JgC,T;) 
where K is the adiabatic exponent. 


LIMITS FOR SUBSONIC FLOW THROUGHOUT 


It is of interest to consider a special case that will 
indicate limiting ideal conditions for subsonic flow. The 
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Fic. 3. Velocity diagram where velocity relative to rotor inlet 


equals velocity relative to inlet of downstream stator (Wi = 
ide) 
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following assumptions are’ made: (1) constant axial 
velocity (i.e., by suitable area change); (2) Mach Num- 
ber relative to rotor inlet approaches unity; and (3) 
Mach Number relative to downstream stator inlet ap- 
proaches unity. 

In Fig. 3 is shown the vector diagram whereby the 
relative velocities approaching the rotor blades and 
downstream stator vanes are equal. This condition 
can be satisfied by a wide combination of rotor speed 
and rotational velocity imparted to the air. 

Thus, 


= 1+ (UpU,/JgC,T1) 
(Ur + U,)/2 = W; cos B; 
Ur = cos Bi U, 
UrU; = 2W, cos B,U, — U,? = 
2W.U,[V 1 — (V/W,)*] — 


By introducing the speed of sound, a, and the gas 
constant, Ro, it follows that 


T, = a*/KgRo 

T2 | M,\? 

14.(K —1)| - (8) - ue] 
T, ( ) i 
where 

M; = W,/a, M, = ,/a, M, = V/a 


The combination for maximum ideal temperature rise 
can be found by differentiating. Thus: 


T2/T, = 1 + (UrU,/JgC,T;) 
U;,U, = 2W; cos B,U,-— U,? 


For max. 72/T,, UrU, = max. 


0(UzU,)/0U, = 2W; cos B; = 
U, = W, cos B; = Ur 
T2/T; = 1 + (W? cos? B;/JgC,T;) 
T, = a?/KgRo 


Ts KRo 2 2 
TC, cos? = 1 +5 M? = 


+ (K 1)M,? 


where M, = W,/a and M, = U,/a. It is seen that the 
theoretical maximum occurs when the absolute inflow 
into the rotor is axial (Fi ig. 4) and the circumferential 
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Fic. 4. Velocity diagram for maximum stage pressure for case 
when W; = Wize. 
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Fic. 5. Velocity diagram for impulse-type rotor blading and 
supersonic velocity entering downstream stator. 
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velocity of the rotor blade element is equal to the rota- 
tional velocity imparted to the air by the rotor blade 
element. 

The general case represented by Fig. 3 can be illus- 
trated by plotting ideal temperature ratio against rota- 
tional tip speed, as shown in Fig. 6. As an illustrative 
case, specific values were assigned to the relative veloc- 
ity entering the rotor blades and downstream stator 
vanes. 


SUPERSONIC FLOW IN STATOR VANES 


In the previous section, special cases of limiting ideal 
stage pressure ratios were considered, whereby veloci- 
ties relative to both rotor and stator were limited by 
the speed of sound. Further increase of stage pressure 
would involve lengthening the horizontal vector, which 
would give rise to supersonic velocities entering either 
the rotor or stator. The special case that will now be 
considered is the case of limiting subsonic flow in the 
rotor and supersonic flow entering the downstream 
stator.* The corresponding vector diagram is shown in 
Fig. 5, where is taken greater than 90° as an example. 
For this more general case, 


T2/T, = 1 + (U,U,/jgC,T1) 


U, = W, cos 8; — (W; sin B,/tan Bo) 
= WM1 — (V/W,)? — (V/tan bo) 
T; = a?/KgRo 
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- 5 M. The rotor blade is a pure impulse type, whereby the 
TT, 1+ (K — 1) | vir ieee Se | M, entrance angle of the rotor blade referred to the plane 
of rotation is 29° and the trailing edge, 151°. The Mach 


Since the tip speed of a rotor is usually chosen as 
high as the strength limitations allow, it is of interest 
to plot curves showing temperature or pressure ratio 
as a function of air outlet angle for constant values of 
circumferential speed. This is shown in Figs. 7 and 8 
for two values of circumferential Mach Number of the 
blade element. 

In each case, the limiting deflection angle is taken as 
equal and opposite to the inlet angle, and the relative 
exit velocity is equal to the limiting relative inlet veloc- 
ity; in other words, the compressor would be analogous 
to the pure impulse turbine. 


SUPERSONIC FLow THROUGH ROTOR AND STATOR: 
Supsonic AXIAL VELOCITY 


If the restriction of subsonic flow is removed from 
both rotor and stator, still further increase of ideal stage 
pressure is theoretically possible. Limitations then are 
associated chiefly with the ‘axial velocity and rotor tip 
speed. A family of curves is plotted in Fig. 9 for sub- 
sonic axial velocity and a blade element circumferential 
speed of 1,600 ft. per sec. 

This general case was proposed by Tsien* who worked 
out a specific illustrative example as follows. 


Number of the axial approach flow is 0.92. The Mach 
Number of the relative flow entering the rotor blade is 
1.62 and that of the flow leaving the blade is 1.55. The 
flow entering the downstream stator vanes-has a Mach 
Number of 3 at an entrance angle of 14.5°. Dr. Tsien 
estimated the pressure rise for the entire stage by first 
computing the pressure drop in accelerating the axial 
flow from atmospheric stagnation pressure to the section 
immediately preceding the rotor; then assuming no 
pressure change in the rotor blades; and, finally, esti- 
mating the pressure recovery in a special design of 
downstream stator vanes utilizing multiple oblique 
shock waves. Dr. Tsien’s calculations made in this 
manner show a possible compression ratio of 15 in a 
single stage with an estimated adiabatic compression 
efficiency of 65 per cent, which would not seem unfavor- 
able in the light of the extremely high compression ratio 
involved. 


COMPLETELY SUPERSONIC FLOW 


In this case the axial flow entering the compressor is 
assumed supersonic, as well as the flow relative to both 
rotor and stator. The only limitations on the stage 
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pressure or temperature ratio, will be the structural 
limitation on the rotor tip speed. Again taking an 
assumed constant value of 1,600 ft. per sec. for the blade 
element circumferential speed, several calculated ideal 
curves are shown in Fig. 10. 


CONCLUSION 


In conclusion, it is emphasized that the present dis- 
cussion intends only to point out that theoretically 
high ideal stage pressure ratios are attainable with the 
use of highly cambered blades of the curved lattice type. 
It is realized that highly cambered blades are not widely 
used in practice because of the tendency toward lower 
efficiency. However, it seemed worth while to point 
out the large potential increase in stage pressure ratio 
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theoretically possible with a view toward stimulating 
experimental investigations of the losses involved in 
‘such a system and of possible methods of decreasing 
these losses. Application would be especially promising 
for cases where simplicity and compactness outweigh 
efficiency in importance. 
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An incorrect graph (Fig. 9, page 63) appeared 
in the paper “The Kinematics of Turbulence” 
by F. N. Frenkiel in the JOURNAL OF THE AERO- 
NAUTICAL SCIENCES, Vol. 15, No. 1, pp. 57-64, 
January, 1948. The correct figure is repro- 
duced here. 


The following errors also appeared in the paper: 

(1) Table 1, column 3, “‘instant’’ should be 
“instant ¢.”’ 

(2) Table 1, column 4, “R,(x) = R,(x)” 
should be ‘‘R,(x) = R,"(x).”’ 

(3) Page 59, in Eq. (6), “®,”’ should be 

(4) Page 61,line 9,“ should be 

(5) Page 62, column 1, line 5, “‘Agm? +1 + 
m? + 1”’ should be ‘‘(Aogm? + 1)/(m? + 1).” 

(6) Page 63, line 3, “8 exp” should be ““B 
exp.” 
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Fic. 9. Spectrum of turbulence measured by Dryden." 
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Services 
of the Libraries of the 


Institute of the Aeronautical Sciences 


The services of the Libraries are available to all mem- 
bers of the Institute, to Corporate Members, to advertisers 
in the ENGINEERING Revigw and Agro- 
NAUTICAL ENGINEERING CaTALoa, and, under usual library 
limitations, to the public. Four specialized services are 
available. 


The Paul Kollsman Lending Library 


This lending library service makes available, without 
charge, the latest and more important aeronautical books. 

Members may request the loan of any acronautical or 
technical book they wish to borrow. Through an exchange 
agreement with the Engineering Societies Library, any 
book on general engineering may be borrowed from its great 
collection of over 160,000 volumes. 

A photostating service is available at usual library rates. 

Applications for membership in the library and further 
information will be sent on request. 


The W. A. M. Burden Reference Library 


This reference library contains over 12,000 aeronautical 
books, magazines, pamphlets, and reports gathered from 
world-wide sources and is one of the most complete aero- 
nautical libraries in the world. Material from this library 
is not available for loan but may be used for reference pur- 


poses. 


The Pacific Aeronautical Library 


6715 Hollywood Boulevard 
Los Angeles 28, California 


Established in cooperation with the aircraft companies 
the library serves. The leading aircraft companies in or 
near Los Angeles participate in its support and operation. 

This service library for aeronautical research is available 
to the public for reading privileges. Source material in- 
cludes aerodynamic and structural research reports, as well 


as books on drafting, production methods, history, and al- 
lied sciences. It furnishes books, periodicals, and pamphlet 
material to the participating aircraft companies to supple- 
ment their engineering libraries. 


Technical Information Service 


* This service has experienced personnel under the super- 
vision of trained aeronautical engineers to compile any in- 
formation desired. The services range from listing special- 
ized reference books to the preparation of exhaustive 
bibliographies, digesting of reports, and general surveys of 
any acronautical subject. Some of the available services 
are: 


Bibliographies on any aeronautical subject. 

Reports on any aeronautical subject. 

Digests of aeronautical books, papers, periodicals, and refer- 
ences. 

Translations. 

Engineering investigations of special aeronautical subjects. 

Biographies of individuals engaged in aeronautics. 

Photostats of any acronautical or general engineering 
material. 

Microfilms made on special order. 

Photographs made from the Institute’s photographic collec- 
tion. 

Drawings and tracings made. 


In addition to the services mentioned any commission 
which comes within the scope of the Service will be ac- 
cepted. Special arrangements may be made for work re- 
quiring several weeks or months. 

Translators are available for accurate transcriptions of all 
foreign language data. Translations are carefully edited 
by trained engineers. 

Reproductions of any material in the collections of 
the Institute may be ordered at standard photostat rates. 
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